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Bestimmung von ebeuien sowie von kreis- und 
kugelsymmetrischen Raumladungsfeldern mit Hilfe einfacher 
Widerstandsketten mit zusitzlichen Stromquellen 


Von GREGOR CREMOSNIK und Max J. O. Strutt, Ziirich!) 


1. Einleitung 


Bei den Elektronenréhren treffen wir oft eindimensionale Felder an, wie ie 
das bei der ebenen Diode, der Zylinderdiode und bei kugelsymmetrischen 
Elektrodensystemen der Fall ist. Analoge Messungen kann man auch an zwei- 
dimensionalen Netzen (Widerstandsnetze fiir ebene und rotationssymmetrische ak 
Felder) ausfiihren [4, 5, 6]*), doch ware es von praktischen Gesichtspunkten 
aus eine bedeutende Aufwand- und Zeitersparnis, wenn man die Messungen an 
gewohnlichen Widerstandsketten ausfiihren kénnte. Dies ist durchaus méglich, By 
wie folgende theoretische Betrachtungen sowie praktische Messungen an ~ es 
solchen Widerstandsketten zeigen sollen. ae 
Die Poissonsche Gleichung fiir allgemeine Bacneion ae Felder kann man’ 

in der Form oes ane 
| AV = Pte (1) eg 


Or? Y Or € 


cre. 


schreiben. Hierin bezeichnet V das Potential, 7 die eindimensionale Koordinate, 
o die Raumladungsdichte und ¢ die Dielektrizitatskonstante, die im Fall der 
Hochvakuum-Elektronenréhren gleich ¢, ist. Auf der rechten Seite haben wir — 
das positive Vorzeichen, da wir die negative Elektronenraumladungsdichte 
—o beriicksichtigt haben. Nun kann man leicht einsehen: 
a) fiir »y = 0 erhalten wir die Raumladungsgleichung fiir das ebene eindimen- 
sionale Feld; . 
b) fiir y = 1 folgt die Raumladungsgleichung fiir das kreissymmetrische Feld;. 
c) fiir vy = 2 folgt die Raumladungsgleichung fiir das kugelsymmetrische Feld. 


2. Die Widerstandskette als Analogiegerat Be 

Rs 

Eine Widerstandskette als Analogiegerat fiir die Lésung der eindimensio- os 
ee 


_ nalen Raumladungsgleichung ist in Figur 1 skizziert worden. In der Reihen- 
 folge der Widerstande sind jene mit ungeraden Indizes verwendet worden, in : 


1) Institut fiir hhere Elektrotechnik der ETH. Pai eee 
2) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 359. vy 


a 
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Analogie zu der zweidimensionalen Anordnung, wo die zweite Richtung gerade 
Indizes erhalt [4,5]. Der Einspeisestrom 7 soll die vorhandene Raumladung 
nachbilden. Die Widerstandskette stellt eine Analogie zur Differenzenform der 
Poissonschen Gleichung (1) dar. 


ee OP See es 
_— SP ; << 
4 Ig '9 
Figur 1 


Widerstandskette zur Lésung von eindimensionalen Feldproblemen mit Nachbildung der Raum- 
ladung durch zusdtzliche Speisestrome. ; 
Die Differenzenform der zweiten Ableitung ist nach Figur 1 


8a Vee ea 
Ove ad : (2) 


wo d der gegenseitige Abstand benachbarter Maschenpunkte (zum Beisp:el — 


mit den Potentialen V, und V,) ist. Fiir das Glied erster Ordnung der Glei- 
chung (1) ist die Differenzenform 


» OV y 


ater aires Hl 


Vz— V1) » (3) 


da man das betrachtete Intervall tiber zwei Maschenweiten (2d) erstrecken 


muss, analog wie die zweite Ableitung. Somit kann man die Differentialglei- 
chung (1) in der Differenzenform schreiben, indem man die Ausdriicke aus den 
Gleichungen (2) und (3) beriicksichtigt: 


AV = = [N+ %-2%+ 32 (%-V)] =. (4) 


v 
27 € 


Schreiben wir noch die Differenzengleichung (4) in der Form, die fiir die 
Analogiebestimmung vorteilhaft ist: 


v% (1-32) +%, (14+ 34) -2y%,-£a=0. (5) 


Um den Analogieausdruck zur Gleichung (5) zu erhalten, wenden wir den 


Kirchhoffschen Satz fiir die Stréme im Knotenpunkt V, (Figur 1) an, indem- 


wir die Widerstande R, und R, zu 
R 


R 
Re as Soak und 135 a 6 Bi (6) 7 


1 a eae 
ZY wares 


wahlen, wobei R ein noch zu bestimmender Widerstandswert ist. Nach Figur 1 
folgt dann d 


(1-35) +%(1+ 44) -2%4+8R=0. (7 


—— 
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Aus dem Vergleich der Gleichungen (5) und (7) erkennt man, dass sie voll- 
kommen identisch sind, wenn der Einspeisestrom i, zu 


<> eo -@ 

gare (8) 
gewahlt wird. Das negative Vorzeichen zeigt, dass der positive Strom aus der 
Kette herausfliessen muss. R ist ein Widerstand, der im allgemeinen von der 
Lage des Punktes im Feld abhangt. 

Nehmen wir nun nach Figur 1 an, dass der Punkt V, ein allgemeiner 
Punkt m im Feld ist, so wird der Punkt V, zum (m + 1)-ten Punkt. Fiir den 
m-ten Punkt gilt nach Gleichung (6) mit der Widerstandsbezeichnung R = 1/A,, 

1 v 1 v 
R740 (1= 55) md 9 = 4n(14+4h), (9) 


2m 


wobei das Verhaltnis r/d mit m bezeichnet wurde. In Gleichung (9) ist R, als 
oberer Widerstand fiir den Punkt m angenommen. Nun sollte sich der gleiche 
Widerstandswert fiir R, ergeben, wenn wir ihn als unteren Widerstand fiir den 
Punkt (m+ 1) annehmen: 


& Age (t= tera): (10) 


Da aus den Gleichungen (9) und (10) der gleiche Wert fiir R, folgen sollte, muss 
die allgemeine Beziehung fiir alle Felder gelten: 


g ES (2m+2—v)m 
Ama  (2m+ rv) (m+ 1)" 


(11) 


Mit Hilfe dieser wichtigen Bezichung kann man dann die Werte der Wider- 
stande sowie auch die Einspeisestréme 7, fiir verschiedene Widerstandsketten 
bestimmen: 

a) Fiir das ebene, eindimensionale Feld ist » = 0. Wenn wir y = 0 in die 
Gleichung (11) einsetzen, erhalten wir 


ee eee x oder Ay = Anis = ==) (12) 
0 


An4a 
In einer Widerstandskette fiir das ebene Feld sind alle Widerstande gleich, 
was schon von der Konstruktion des zweidimensionalen Widerstandsnetzes 
her [4,5] bekannt ist. Mit der Widerstandsbezeichnung K = 1/A,, folgt aus 
Gleichung (8) fiir den Einspeisestrom der Wert 


iy = = d? A, . (13) 
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Fiir das ebene Feld ist A,, = 1/R,, wobei Ry der Wert der Widerstande in der 
Kette ist [vgl. Gleichung (12)]. Somit wird nach Gleichung (13) der Einspeise- 
strom 2 

to =—_— @Ro OQ 'e (14) 
Die Raumladung o kann man durch die bekannte Stromdichte 7 und durch die 
Elektronengeschwindigkeit v ausdriicken: 


o=1. YS a5) 


Aus dem Energiesatz 

1 
— 2S ey. 
zy my 


(e Betrag der Elektronenladung, m Elektronenmasse) erhalten wir die Ge- 
schwindigkeit wag Ts 
g=\2 Pe (16) 


wobei e/m = e, geschrieben wurde. Unter Beriicksichtigung der Gleichungen (15) 
und (16) erhalten wir fiir den Einspeisestrom 7) nach Gleichung (14) 


wee a? ee ae SN 7 
ip a@erk, pve es (17) 


b) Im Falle des kreissymmetrischen Feldes ist y = 1. Wenn wir v = 1 in die 
Gleichung (11) einsetzen, erhalten wir die Beziehung 


A m 


Agi? tee 


oder 
Aa oes Mens ae ; (18) 


m m+1 


Giinstig ist, alle Widerstande in der Kette auf den ersten Widerstand zu be- 
ziehen. In diesem Fall ist das der Widerstand zwischen dem Punkt an der 
Rotationsaxe und dem nachsten benachbarten Punkt. Aus der Gleichung (9) 
fiir den Widerstand R, (m= 1, y= 1) folgt unter Beriicksichtigung der Glei- 
chung (18) ; 


a is 


R,=7=R. (19) 


Somit folgt aus der Gleichung (9) fiir den Widerstand R, 


) R 
BL RE eee (20) 


Fir m = 2 erhalten wir 
me 
. 5) 
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und somit werden die Widerstande der Reihe nach 


ee 
Ros a es sumte " (21) 


Fir den Einspeisestrom folgt unter Beriicksichtigung der Gleichungen (18) 
(19) aus der Gleichung (13) der Wert 


2m a? 


ip — Qe So. (22) 


é€ 


Hierbei ist m= vy/d und R, als der gewahlte erste Widerstand in der Kette 
bezeichnet. Fiir die Stromdichte 7 gilt im kreissymmetrischen Fall die Beziehung 


1%= Jun » (23) 


wobei 7;, die bekannte Stromdichte an der Stelle k ist (zum Beispiel an ‘der 
Kathode). Wenn wir die Gleichung (15) fiir die Raumladungsdichte @ in die 
Gleichung (22) einsetzen, erhalten wir unter Beriicksichtigung der Gleichung 
(16) fiir den Einspeisestrom Z) den Wert 


pte: 2a? CA etl 

: sVee, BIS es or 
c) Eine analoge Betrachtung fiir das kugelsymmetrische Feld (vy = 2) ergibt 

nach Gleichung (11) die Beziehung 


(24) 


oder 


ge a 7 ae aa 

laa capamane ati 3 (25) 
Nach den Gleichungen (9) und (25) folgt fiir den ersten Widerstand R, der 
Kette der Wert 


i age A a te By (26) 


Da der erste Punkt ausserhalb des Kugelzentrums m= 1 ist, wird R, = oo. 
Das Feld (Potential und Feldstarke) im Zentrum der Kugel nimmt auch 
unendliche Werte an. Da wir in Wirklichkeit bei einer Elektrode im Zentrum 
immer mit einem endlichen Halbmesser zu rechnen haben, stért das die Feld- 
messung nicht. Der erste Punkt der Kette liegt dann auf der zentralen Elek- 
trode. Fiir den Widerstand R, erhalten wir [nach Gleichung (9)] 

i dL 


Rate ete Guay cmatet 8s aa 
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Hier ist es giinstig, alle Widerstande auf den Widerstand R; zu beziehen. Fur 
m = 2 folgt R; = R,/3, und fiir héhere m werden die Widerstande der Reihe 
nach Plena eee 

CO, Ry; ek Son 10’ ed ag (28) 
Bestimmen wir noch die Einspeisestroéme als Nachbildung der Raumladung. 
Fiir das kugelsymmetrische Feld erhalten wir analog zur Gleichung (22) 


’ m2 ah 
19 = a ag (29) 


Fiir die Stromdichte gilt jetzt die Beziehung 
j=, Gh (30) 


Aus der Gleichung (29) folgt unter Beriicksichtigung der Gleichungen (15), (16) 

und (30). phe aD ( Vip y Tie ; (31) 
2 02, Ry Wed VV 

wobei 7, wieder die bekannte Stromdichte bei 7 = 7; ist. 

Hiermit haben wir die praktisch wichtigen Ausdriicke fiir die Anwendung 
der Widerstandsketten bei den Feldmessungen erhalten. Der grosse Vorteil der 
Widerstandsketten liegt in ihrer einfachen Konstruktion sowie in der leichten | 
und raschen Durchfiihrbarkeit der Messungen. 

Bei der Bestimmung des Raumladungsfeldes treten oft Fehler vor den 
Elektroden auf. Wenn man das ganze Feld bei Nachbildung der Raumladung 
in Differentialflachenelemente aufteilt, sieht man, dass die Raumladung vor 
den Elektroden mit den Speisestrémen nicht berticksichtigt wird. Man muss — 
die Speisestréme in den Punkten vor den Elektroden korrigieren oder die 
entsprechende Maschenweite vor den Elektroden zu d/2 wahlen (naheres in 
[15]). Dieses Ergebnis stammt von Herrn P. TscHopr in diesem Institut. 


3. Ebene Diode 


Als praktische Anwendung untersuchen wir zuerst den Fall der ebenen 
Diode. Bei ebenen Dioden kann man das Feld genau berechnen. So erhalten’ 
wir eine Méglichkeit, die praktischen Messergebnisse mit den genauen berech- 
neten Werten zu vergleichen. 

Bei der Bestimmung des Raumladungsfeldes ist das Verfahren der schritt- 
weisen Naherung anzuwenden. Fiir die Bestimmung der Einspeisestrome als 
Nachbildung der Raumladung sollte die Stromdichte 7 sowie auch die Poten- 
tialverteilung V bekannt sein. Die Stromdichte 7 kann man aus der bekannten 
Anodenspannung U, und aus der Geometrie der Anordnung berechnen. Die 
Potentialverteilung unter Beriicksichtigung der Raumladung ist zundachst 
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nicht bekannt. Ausgehend von dem gemessenen statischen Feld Vy (ohne 
Raumladung) bestimmt man die entsprechende Raumladung nach Gleichungen 
(15) und (16) bzw. den Einspeisestrom nach Gleichung (17). Nach der Nach- 
bildung der Raumladung mit diesen Einspeisestrémen Andert sich die Poten- 
tialverteilung. Mit den neuen Potentialwerten bestimmt man wieder die Raum- 
ladung bzw. die notwendigen Einspeisestrome usw. Endlich ist die Anderung 
des Potentialverlaufs bzw. der einzuspeisenden Stréme so klein, dass sich die 
Potentialverteilung nicht mehr andert. Das Verfahren nahert sich im allge- 
meinen sehr rasch einem Endwert. Bei kleiner Raumladung sind nur wenige 
Schritte notwendig, bei starkerer Raumladung jedoch mehrere. Modelle mit 
Widerstandsketten haben nur wenige Messpunkte im Feld, so dass die Durch- 
fiihrung der Messungen wenig Zeit in Anspruch nimmt. 


3.1. Austritisgeschwindigheit der Elektronen aus der Kathode gleich Null 


Zunachst erlautern wir kurz den Einfluss der Vergrésserung im Modell bei 
der Bestimmung der Einspeisestréme. Aus der Differentialgleichung (1) geht 
folgendes hervor: Wenn wir die linearen Dimensionen u-mal vergréssern, so 
miissen wir auf der rechten Seite der Gleichung (1) durch m? dividieren. Das 
bedeutet, dass man die Raumladung ?-mal kleiner machen muss. Nur dann ist 
die Poissonsche Gleichung am Modell erfiillt. Der Einfluss von Fehlern in den 
Widerstanden des Modells auf die Messergebnisse sowie der Einfluss der end- 
lichen Maschenweite wurden an anderer Stelle erértert [4,5]. Die verwendeten 
Widerstande des Modells wurden alle innerhalb + 0,2% auf die errechneten 
Nominalwerte eingestellt. 

Der einfachste theoretische Fall einer ebenen Diode wird nur kurz behan- 


delt. Die Stromdichte ist gegeben durch 
3/2 


j = 2,34- 10-8, (32) 


wobei U, die Anodenspannung gegen die Kathode und } den Abstand zwischen 
den Elektroden bezeichnen. Der Potentialverlauf ist theoretisch durch 


u=u (=) (33) 


gegeben. Als Beispiel nehmen wir eine ebene Diode mit 
b=100um=10-?cm und U,=30V. 

Am Modell wurde die gleiche Spannung und ein Abstand von 10 Maschenweiten 

zwischen den Elektroden angenommen. Die Vergrésserung wird: 


20 
Oe 


Y= =P Sees 
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Aus der Gleichung (32) folgt fiir die Stromdichte 
j = 3,845 A/cm?. 


Diese Werte konnen bei praktischen Dioden nicht auftreten, die Stromdichte 
j ist stets viel kleiner, ebenso die Anodenspannung U,. Durch obige Annahmen 
erhalten wir jedoch eine gréssere Raumladung, was fiir die Untersuchung und 
Priifung des Verfahrens giinstig ist. Aus der Gleichung (17) kénnen wir die er- 
forderlichen Einspeisestréme bestimmen. Mit d = 2cm und R, = 3,5 kQ wird 
c = 217,768 [Gleichung (17)] und somit unter Beriicksichtigung der Vergrésse- 
rung ” 

1g = 54,442 . ——— w —— A. (34) 


Die Resultate der schrittweisen Naherungen, ausgehend von dem Potentialver- 
lauf (0) ohne Raumladung sind aus der Figur 2 ersichtlich. Fir verschiedene 


0 
Fal aes 40 
b 
Figur 2 
Schrittweise Naherungen bei den Messungen des Potentialverlaufs U einer ebenen Diode mit Aus- 
trittsgeschwindigkeit der Elektronen aus der Kathode gleich Null. Kurve 0: Homogenes statisches 


Feld ohne Raumladung; I: erster Schritt mit Nachbildung der Raumladung; Br: berechneter — 


Potentialverlauf und zugleich dritter Naherungsschritt am Modell. 


Abstande x/6 kénnen wir die theoretisch berechneten Potentialwerte Br nach - 


Gleichung (33) sowie die nach der dritten Naherung erhaltenen Messwerte U- 
in der folgenden Tabelle vergleichen: 


x[b | 0,1 | 0,2 | 0,3 | 04 | 05 | 0,6 | 07 | 08 | 0,9 
at Ne OE SE a eS a 


Br[V]| 1,392 | 3,51 | 6,03 | 8,85 [11,91 | 15,18 118,65 
| ULV) 1,41 | 3,54 | 6,08 | 8,85 | 11,85 | 15,00 | 18,50 


26 07 
26,00 | 


222i 
fafa AO) 


ON i ol ol Tg i en ’ Ad 
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In der Figur 2 ist nur die erste Naherung I und die letzte dritte Naherung, die 
mit dem berechneten Potentialverlauf By gut lbereinstimmt, gezeichnet sui 
den. Die zweite Naherung II weicht so wenig von der dritten Naherung bzw. von 
dem berechneten Verlauf ab, dass man sie schwer einzeichnen kann. Die pro- 
zentualen Fehler im Potentialverlauf fiir verschiedene Naherungen sind in der 
Figur 3 angegeben worden. Fiir die Abweichung des statischen Feldes A, gilt 


aS x ] F 


> 


Aq” Statische's Felli 


Grenze 
Instrum. 
Goerz 


Figur 3 
Prozentualer, relativer Fehler A [°] der schrittweisen Naherungen. Fiir Jy (statisches homogenes 
Feld) gilt der linke vertikale MaBstab, fiir 4, (nach dem ersten Schritt) gilt der rechte MaBstab, 
ebenso fiir A, (nach dem zweiten Schritt) und A, nach dem dritten Schritt. Eingezeichnet: die 
Grenze der Messgenauigkeit des verwendeten Universalinstrumentes. 


der zehnfache vertikale MaBstab (linker MaBstab in Figur 3). Aus dieser Figur 
sieht man deutlich, wie die Fehler bei den schrittweisen Naherungen jeweils 
herabgesetzt werden. Die Fehler der dritten Naherung A, fallen schon in den 
Bereich der Messgenauigkeit des verwendeten Universalinstrumentes, dessen 
Grenze in Figur 3 angedeutet ist. 

Wenn wir die prozentualen relativen Fehler in Abhangigkeit von der Zahl 
der ausgefiihrten schrittweisen Naherungen aufzeichnen, erhalten wir den 
Verlauf nach Figur 4. Aus der Figur 4 geht deutlich hervor, dass die Abnahme 
des Fehlers A nach den einzelnen Schritten ungefahr einer Exponentialkurve 
entspricht. Somit kénnen wir mit guter Naherung diesen Verlauf durch die 
Funktion 

A = K exp(—« S) (35) 
darstellen, wobei S die Nummer der Naherung (S = 0, 1, ... ) bezeichnet und 
der Exponent « aus den Kurven hervorgeht. Fur schwache Raumladung ergab 
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sich, dass der Exponent « ungefahr den mittleren Wert 2 hat (1 < « < 3)i Bei 
starkerer Raumladung sinkt « und erreicht Werte, die kleiner als 1 sind, wie 
wir im Abschnitt 3.2 sehen werden. 


eS p= 


aj} 
: 


A{%] 


}— 
10 %=01 
)2 
5 
(0) 
19 le | 
(6) Ue I ill 


a Ss 
Figur 4 
Prozentualer relativer Fehler A [%] der schrittweisen Naherungen als Funktion der Nummer S des 


Schrittes: O, I, II und III fiir die Maschenpunkte vor der Kathode mit x/b = 0,1, 0,2... (Lage 
des Punktes als Parameter). 


3.2. Konstante Austrittsgeschwindigkeit der Elektronen aus der Kathode 


Betrachten wir den Fall, wo sich zwischen Kathode und Anode infolge der 
vorhandenen Raumladung ein Potentialminimum ausbildet, so dass der Po- 
tentialwert im Minimum U,, > 0 ist und die Anodenstromdichte j, gleich der 
Kathodenstromdichte 7 (keine riickkehrenden Elektronen). 


Fur die theoretische Betrachtung ist es vorteilhaft, eine normierte Strom- 
dichte 7* einzufiihren : 


jot, (36) 
3 
wobei die Stromdichte j, nach der Gleichung (32) : 
jet let (37) 
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ist, 7 die vorhandene Stromdichte in der Diode, unter Beriicksichtigung der 
konstanten Austrittsgeschwindigkeit der Elektronen aus der Kathode, und U, 
die Anodenspannung. Damit der obige Fall mit dem Potentialminimum 
zwischen Kathode und Anode vorliegt, darf die Stromdichte 7 nicht beliebig 
gesteigert werden. Die bekannte Theorie [1, 2, 3] ergibt, dass die maximale 
Stromdichte nach Gleichung (36) 


ite= (14/2) (38) 


betragt. Dabei bedeutet U, die Spannung an der Kathode, die der konstanten 


_Austrittsgeschwindigkeit entspricht, und U, die Anodenspannung im Falle, dass 


U, > U, ist. Fiir die Spannung des Potentialminimums folgt dann 


U,/U, __ ; 


sen = U, ———— (39) 
* (1+ VU/U,)? 
und fiir die Lage 
Say 2 ES od Ag (40) 
: (1 + VU,/U;) 


Somit ist es méglich, nach den Gleichungen (38), (39) und (40) aus den vor- 
gegebenen Elektrodenspannungen U, und U, die noch zulassige Stromdichte, 
die Spannung des Potentialminimums sowie seine Lage zu berechnen. 

Als Beispiel nehmen wir eine ebene Diode mit dem Elektrodenabstand 
b= 100 wm, mit einer Anodenspannung U,= 25 V und mit der Kathoden- 


mere U, = 0,01 U, = 0,25'V. 
Aus der Gleichung (37) folgt die Stromdichte 
§g td, 04-2 10—* +125 = 2,925 Ajem*. 


Wenn wir den Fall der maximal zuldssigen Stromdichte nehmen, so erhalten 
wir aus der Gleichung (38) 


Gt = (2 POA* = 1,331 
und aus der Gleichung (36) die vorhandene Stromdichte 7 
j = 2,925 - 1,331 =.3,8932 A/cm?. 
Fiir die Spannung des Potentialminimums erhalten wir nach der Gleichung (39) 


Um, = 0,2066 V 


und fiir die Lage in bezug auf den Elektrodenabstand 6 nach Gleichung (40) 


oa ee oe. 


Xm — 
a ed 0,0233 . 
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Da die Anodenspannung in bezug auf die Kathodenspannung sehr gross ist, 
tritt ein Spannungsminimum in der Nahe der Kathode auf. Wir haben einen 
Fall, wo die Raumladung nicht gross ist. 

Bei der Berechnung des Potentialverlaufs integrieren wir die Poissonsche 
Gleichung 


CUS .9¢. 2 3 Jae (41) 
dx* £0 V2 & & VU 


Unter Beriicksichtigung der Integrationskonstanten ergibt sich fiir den ge- 
suchten Potentialverlauf die Formel [2,3] 


+0 b (0 Hg)® E+ [le — ae? E (oR + Ee — a) FS (42) 


ae Va + (x — aq)? = [w= x)* (US + Ee aibas 
wobei 


a= eg V2 6, = 2,34 - 10-8 A/VEP 


ist. Die praktische Berechnung des Raumladungsfeldes ist somit schwierig und 
verlangt sehr genaue numerische Rechnungen wegen des Differenzgliedes in 
der Gleichung (42). Im Vergleich hierzu ist die Ermittlung des Potentialverlaufs 
mit Modellmessungen an der Widerstandskette mit den zusatzlichen er 
strémen viel einfacher und weniger zeitraubend. 

Als Modell nehmen wir. 10 Maschenweiten zwischen den Hiekadee mit 
einer Vergrésserung » = 2+ 10%. Nach der Gleichung (17) unter Beriicksichti- 
gung der Vergrésserung folgt dann fiir die Einspeisestréme 


ty = 54,442 -—_ » —= pA. (43) 


Nun wurde das beschriebene Verfahren der schrittweisen Naherungen ange- * 
wandt. Die Messergebnisse sind in der Figur 5 gegeben worden. Der fiinfte 
Schritt bzw. die fiinfte Korrektur der Raumladung ergibt einen Potentialver- 
lauf, der sehr gut mit dem berechneten nach Gleichung (42) iibereinstimmt. In 
der nachfolgenden Tabelle sind die gemessenen Werte U nach der fiinften 
schrittweisen Naherung sowie die berechneten Werte U, nach Gleichung (42) 
zum Vergleich angegeben worden. Zuletzt ist noch der prozentuale Fehler A 


angegeben [Abweichung der gemessenen Werte von den genau berechneten 
Werten nach Gleichung (42)]: 
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Potentialverlauf U [V] bei schrittweisen Messungen im Modell einer ebenen Diode mit konstanter 

Austrittsgeschwindigkeit der Elektronen aus der Kathode U, = 0,01 Uy. O: homogenes Feld; 

I: erste Naherung bei Nachbildung der Raumladung; II: zweite Naherung ... usw. Br: berechneter 
Verlauf nach Gleichung (42) und zugleich fiinfte Naherung. 


x|b ! 0,1 | 0,2 | 0,3 | 0,4 | 0,5 | 0,6 | 07 | 08 | 0,9 


U [Vv] | 0,64 | 2,03 | 3,93 | 6,29 | 8,85 | 11,70 | 14,90 | 18,25 | 21,45 
U, [V]| 0,63 | 2,01 | 3,90 | 6,25 | 8,81 | 11,66 | 14,95 | 18,32 | 21,50 
A[%]| 1,5 | 1 9 aan A | 


/ ‘ 
Der grésste Fehler tritt fiir die Punkte vor der Kathode auf. Wenn eine sechste 
Korrektur der Raumladung bzw. des Einspeisestromes in den Punkten vor der 
_ Kathode ausgefiihrt wiirde, so ware der Fehler iiberall unter 1%. 

Untersuchen wir noch ein Beispiel mit grésserer Austrittsgeschwindigkeit - 
der Elektronen und somit auch mit grdésserer Raumladung. Die Anodenspan- 
nung soll unverandert bleiben: U, = 25 V. Die Kathodenspannung, die der 
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Austrittsgeschwindigkeit der Elektronen aus der Kathode entspricht, wird 


jetzt 
U, = 0,20, = 5 V: 


Praktisch kénnte dieses Beispiel in Mehrgitterrdhren vorkommen, wo die 
Schirmgitterebene als neue Emissionskathode angenommen werden kann und 
die Elektronen dann grosse Emissionsgeschwindigkeiten haben. Die Strom- 
dichte, welche als Normierungsgrésse dient, bleibt gegen das vorige Beispiel 
unverandert 7, = 2,925 A/cm?, da die Anodenspannung U, die gleiche ist. Die 
maximal zulissige Stromdichte ergibt sich aus der Gleichung (38) 


jee, = (1+V0,2) = 3,03, 
und somit aus Gleichung (36) die vorhandene Stromdichte in der Réhre 
4a = 2,925 3,03 = 8,863 Ajem*. 


Aus der Gleichung (39) erhalten wir fiir die Spannung im Potentialminimum 
den Wert U,, = 2,4 V und fiir die Lage des Minimums in bezug auf den Elek- 
trodenabstand 


ay 0,227 


Das Spannungsminimum liegt jetzt weiter von der Kathode entfernt, und 
gegen die Kathode weist dieser Punkt eine negative Spannung auf im Betrage 
von 


3— 2,4 = 2,6 V. 


Die Vergrésserung des Elektrodenabstandes 6 = 10-2cm am Modell mit 
10 Maschenweiten betragt m = 2-108. Nach der Gleichung (17) folgt somit 
fiir die Berechnung der Speisestréme als Nachbildung der vorhandenen Raum- 
ladung in der Réhre: 


482,5 


j Sea ole (44) 


et Oh 
Den Potentialverlauf kénnen wir mit Hilfe der Gleichung (42) berechnen. Die- 
Messergebnisse der schrittweisen Naherungen sind in der Figur 6 gezeichnet. 
Aus dem Verlauf sieht man deutlich, wie die Naherungen nach 7 Schritten zu 
dem genauen Verlauf fiihren. Bei den ersten zwei Schritten mussten alle Speise- 
stréme korrigiert werden, beim dritten Schritt war bei den letzten zwei Punkten _ 
vor der Anode eine Korrektur nicht mehr notwendig. Bei dem fiinften Schritt 
war dies schon bei den letzten vier Punkten der Fall, und beim letzten Schritt 


‘waren nur in der Halfte vor der Kathode, wo das Potentialminimum liegt, 
noch Korrekturen der Einspeisestrome erforderlich. 


m=. 


Vol. VIII, 1957 Bestimmung von Raumladungsfeldern 343 


25 \ a? sae e En a ee gh 


U [V] 


20 


tio) 


at, 
0 Q2 04 06 08 jo 


Figur 6 


Potentialverlauf U [V] bei schrittweisen Naherungen in einer ebenen Diode mit konstanter Aus- 
trittsgeschwindigkeit der Elektronen aus der Kathode U, = 0,2 U,. Gleiche Bezeichnungen wie 
in Figur 5. 


In der nachfolgenden Tabelle sind die nach Gleichung (42) berechneten 
Potentialwerte sowie die Messwerte U des letzten, siebenten Schrittes gegeben 
worden. In der letzten Zeile ist noch der prozentuale Fehler A berechnet 
worden. 


x|b | 01 | 0,2 | 0,3 | 0,4 


U,[V]| 3,24 | 2,40 | 2,66 | 3,94 | 6,15 | 8,74 | 12,25 | 16,00 | 20,25 
Uv] | 3,26 | 2,42 | 2,68 | 3,96 | 6,10 | 8,71 | 12,25 | 16,05 | 20,25 
FO O68 108 VS OB te O4n | weed cae cies | 


Der Fehler ist somit kleiner als + 1°%, was eine gute und praktisch hinrei- 


_ chende Genauigkeit bedeutet. Dabei ist noch zu sagen, dass der Fehler nicht 
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auf die gesamte Spannung zwischen den Elektroden bezogen wurde, wie es oft 
geschieht, sondern in bezug auf den genau berechneten Wert in dem entspre- 
chenden Punkt. Der obige Fehler ware sonst noch viel kleiner. Aus den Mess- 
ergebnissen folgt deutlich, dass man mit dem beschriebenen Verfahren der 
schrittweisen Naherungen den Betrag und die Stelle des auftretenden Poten- 


75 


50 


25 


Figur 7 
Prozentualer Fehler A [%] der schrittweisen Naherungen der Figur 6 als Funktion des Abstandes 
von der Kathode: x/b. Bezeichnungen der Kurven 4 1: nach der ersten Naherung, usw. 


tialminimums sowie den gesamten Potentialverlauf sehr gut und rasch bestim- 
men kann. Die Berechnung des Potentialverlaufs ist kompliziert, die Resultate - 
sind praktisch auch mit Fehlern gleicher Gréssenordnung behaftet, wie sie sich 
beim letzten obigen Schritt ergeben. Die prozentualen Fehler A fiir verschiedene 
Schritte sind in Figur 7 veranschaulicht worden, A, bedeutet den Fehler nach 
dem ersten Schritt usw. Aus dem Kurvenverlauf geht die rasche Verringerung 
des Fehlers deutlich hervor. Der maximale Fehler tritt in der Nahe des 
Potentialminimums auf. Wenn wir den prozentualen Fehler als Funktion der 
Nummer S der schrittweisen Naherung aufzeichnen, erhalten wir die Figur 8. 
Da die Raumladung grésser ist als im vorherigen Fall, findet die Abnahme des 
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Fehlers langsamer statt als dort, was zu einem Exponent « nach Gleichung (35) 
fiihrt, der kleiner als 1 ist. Fiir die Punkte direkt vor der Kathode ist die 
relative Abnahme am geringsten. Hier ist « am kleinsten (x ~ 0,7). Fiir die 
Punkte vor der Anode, die in der Figur 8 dargestellt sind (x/b = 0,7 ,... , 0,9) 
ist Sera starker, so dass in diesem Fall der Exponent « fast den Wert 1 
erreicht. 


Figur 8 


Prozentualer Fehler A [%] der schrittweisen Naherungen als Funktion der Zahl S der Schritte 
(S = 0, 1, 2,..., 7), fiir die Maschenpunkte vor der Anode (#/b = 0,5, ... , 0,9). 


3.3. Maxwellsche Geschwindigheitsverteilung der emittierten Elektronen 


Fiir die Nachbildung des Einflusses der Maxwellschen Geschwindigkeits- 
verteilung der emittierten Elektronen aus der Kathode an Widerstandsketten- 
modellen wurde ein approximatives Verfahren entwickelt : 

1. Man bestimmt die Temperaturspannung U; der Kathode, die dem qua- 
dratischen Mittelwert v; der thermischen Geschwindigkeit entspricht, aus den 
Energiegleichungen, wobei die Kathodentemperatur T als bekannt voraus- 
gesetzt wurde: 


mM v2 
kT =¢,Up=——. (45) 


Hierin ist & die Boltzmannsche Konstante k = 1,38 -10-?8 J/Grad abs. 
2. Nun verwendet man das entwickelte Verfahren fiir konstante Aus- 
_trittsgeschwindigkeit der Elektronen, indem man im Modell U, = Ur und 
U, = Ur + U, nachbildet, wie dies im Abschnitt 3.2 durchgefithrt wurde. 


346 GREGOR CREMOSNIK und Max J. O. STRUTT ZAMP 


3. Aus dem Endpotentialverlauf (nach der letzten schrittweisen Naherung), 
welcher mit der vorausgesetzten konstanten Geschwindigkeit Uy erhalten 
wurde, bestimmt man die Korrektur der Speisestroéme vor der Kathode bis zum 
Potentialminimum nach der vorgegebenen Maxwellschen Geschwindigkeits- 
verteilung, indem man den riickkehrenden Elektronen Rechnung tragt. 

4. Da somit nicht alle Elektronen zur Anode gelangen, wie das bei kon- 
stanter Austrittsgeschwindigkeit der Fall ist, muss man die Speisestrome ver- 
ringern nach der Gleichung 

lg = tg — lw, (46) 
wobei 7, der Speisestrom vor der Korrektur ist, der aus dem bekannten Poten- 
tialverlauf bei konstanter Austrittsgeschwindigkeit der Elektronen folgt, und 
F ein Faktor, den man bestimmen muss aus dem am Modell beriicksichtigten 
zuriickkehrenden Strom. 

Man integriert die Maxwellsche Verteilungsfunktion jeweils von einem 
betrachteten Punkt zwischen Kathode und Spannungsminimum bis zum Span- 
nungsminimum. Die zuriickkehrenden Elektronen tragen in diesen Punkten 
zweimal zur Raumladung bei. Das Verfahren der schrittweisen Naherung 
wiederholt man, bis die Integration das richtige vorgegebene Verhaltnis zwi- 
schen dem Anodenstrom und dem Sattigungsstrom ergibt. Das Verfahren soll 
an einem praktischen Beispiel erlautert werden. 

Im Gebiet hoher Temperaturen der Kathode (iiber 1000° K) kann man mit 
guter Naherung annehmen, dass die emittierten Elektronen eine Maxwellsche 
Geschwindigkeitsverteilung aufweisen. Die Zahl der durch eine Flacheneinheit 
pro Zeiteinheit hindurchtretenden Elektronen sei Ny. Dann ist die Anzahl dN 
dieser Elektronen, deren Geschwindigkeit zwischen v und v + dv liegt, 


dN =N => exp(-"7) dv , 


oder, unter Beriicksichtigung der Gleichung (45), 


: - 
aN = My exp v;) dU. (47): 


Nehmen wir an, dass die Zahl N, der Elektronen der Sattigungsstromdichte ie 


entspricht. Die Zahl der zuriickkehrenden Elektronen erhalten wir durch 
Integration 
Um/Ur 


N= N, [ exp(-V) AVN E Bs exp(— "| ; (48) 


6 i 


wobei Um die Spannung im Potentialminimum ist und V = U/U; im Integral 
substituiert wurde. Somit ist dann die Zahl der Elektronen, die zur Anode 
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gelangen, 


U, 
N = N,exp i: gai \ : (49) 
/ 
Da die Zahl N, der Elektronen der Sattigungsstromdichte entspricht, erhalten 
wir fiir die Anodenstromdichte 
N 


ie= Ny Ts (50) 
und fiir die Stromdichte der zuriickkehrenden Elektronen 

Whe 

In = Fyn de (51) 

. . 0 

Somit gilt 
und analog 

Is= Jat Ir + (52) 


Als Beispiel betrachten wir eine ebene Raumladungsdiode mit der Satti- 
gungsstromdichte7, = 100 mA/cm?undeiner Anodenstromdichte 7, =40mA/cm?. 
Diese Zahlen haben keine praktische Bedeutung in Réhren, sondern sollen nur 
als Rechenbeispiel dienen. Die Temperatur der Pastekathode sei 900°C oder 
1173°K mit einer Kathodenoberflache von 1,5 cm?. Der Abstand zwischen 
Anode und Kathode sei 6 = 50 um = 5- 10-3 cm. 

Die Lésung der Poissonschen Gleichung (41) mit Beriicksichtigung der 
Maxwellschen Geschwindigkeitsverteilung der emittierten Elektronen ist nur 
unter numerischer Auswertung der letzten Integration médglich. Die Resultate 
stammen von I. LANGMuIR [7], der die Integrale in Tabellen ausgewertet hat. 
Diese Tabellen wurden von P. H. J. A. KLEYNEN [10] vervollkommnet. Eine 
Erweiterung der Ausfiihrungen von I. LancmuIR hat A. VAN DER ZIEL [11] 
gebracht. Ohne auf die Resultate der Ausfiihrungen naher einzugehen, geben 
wir nur die berechneten Werte fiir ein angenommenes Beispiel an. Die Tempe- 
raturspannung der Kathode erhalten wir aus 


U, = 086 = 10-8 Tae 0,10 Vi. 
Fiir die Spannung im Potentialminimum gilt 
Divan In = 0,092 V 
und fiir die Anodenspannung 
Us = te Ned aoe = 9,064.V. 


Yq ist ein aus den genannten Tabellen entnommener Wert und 7, = U,»/Ur. 
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Fiir die Lage des Potentialminimums erhalten wir mit Hilfe der Tabellen 


x 


Pm is 
b 


= § > = 0,338. 


Die berechneten Spannungen sind klein. Deshalb wurden im Modell alle Span- 
nungen mit dem Faktor 50 multipliziert. Somit ergibt sich: 


Up Ven U, = 3.2: V vund? 0, = tee 


Nach dem beschriebenen Verfahren betrachten wir zuerst den Fall konstanter 
Austrittsgeschwindigkeit mit: U, =5V und U, = Ur + U,= 8,2 V (vgl. Ab- 
schnitt 3.2). Fiir die Stromdichte erhalten wir nach den Gleichungen (36), (37) 
und (38) 

] = To is i 12,4 A/cm?. 


Die Rechnung ergibt dann U,, = 1,58 V oder die Spannung gegen die Kathode 
Um = —5 + 1,58 = —3,42V und fir die Lage des Potentialminimums 
%mlb = 0,404. Die obige Spannung des Potentialminimums unterscheidet sich 
von der Spannung, die aus der Rechnung mit Maxwellschen Geschwindigkeits- 
verteilung der emittierten Elektronen folgt, um mehr als 25%. Die mittleré 
thermische Geschwindigkeit mtisste um einen Faktor vergréssert werden, 
damit die Ubereinstimmung besser wiirde. Diese notwendige Vergrésserung 
kann man nachher bei der Korrektur infolge der Maxwellschen Verteilung 
beriicksichtigen, wie in dem nachfolgenden Beispiel gezeigt wird, indem man 
die Stromdichte vergréssert. 

Im Modell wurden 10 Maschenweiten fiir den Abstand zwischen den Elek- 
troden gewadhlt. Die Vergrésserung betragt m = 20/(5 - 10-3) = 4- 10%. Aus der 
Gleichung (17) folgt unter Beriicksichtigung der Vergrésserung fiir den Ein- 


speisestrom 
217,768. 4 170 


4 . St , == 
0 16 yu jU 


pA. (53) 


Die Messresultate der schrittweisen Naherungen sind in der Figur 9 gegeben 
worden. Nach dem fiinften Schritt kommen wir zur Endkurve, die sehr gut 


mit dem gerechneten Verlauf U,, fiir die genannte konstante Austrittsgeschwin- ~ 


digkeit tibereinstimmt. Zum Vergleich ist in Figur 9 auch die gerechnete 
Kurve Br bei Beriicksichtigung der Maxwellschen Geschwindigkeitsverteilung 
eingezeichnet. Man erkennt gleich die grosse Abweichung der beiden Kurven 
U, und Br voneinander. 


Nun wurde die M-Verteilung (Maxwellsche Geschwindigkeitsverteilung) — 


der emittierten Elektronen berticksichtigt. Der erste Schritt wurde so durch- 
gefihrt, dass Js = Jq Angenommen wurde. Fiir die Punkte vor der Kathode bis 
zum Potentialminimum wurde der Speisestrom 7, gewahlt, der aus der Span- 


i ie 
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Mh 


Figur 9 


Potentialverlauf U [V] bei schrittweisen Naherungen im Modell einer ebenen Diode mit konstanter 
Austrittsgeschwindigkeit der Elektronen aus der Kathode U, = Uf. 1: erste Naherung... usw. 


U;, letzte, fiinfte Naherung, die mit dem berechneten Verlauf gut Ubereinstimmt. Kurve Br: 
berechneter Potentialverlauf bei Maxwellscher Geschwindigkeitsverteilung (M-Verteilung) der 
emittierten Elektronen. 


nung des Minimums nach Gleichung (53) folgt. Bei der Beriicksichtigung der 
zuriickkehrenden Elektronen, die zweimal zu der Raumladung beitragen, wur- 
den die Stromdichte bzw. die Speisestréme vergréssert. Wenn wir die Abkiir- 
zung V = U/U; einfiihren, ergibt sich nach Gleichung (48) die Zahl der zuriick- 
kehrenden Elektronen zwischen zwei Punkten V, und V, 


Vs 
N, = Ny | exp(—V) aV = Ng [exp (—V;) — exp (—Val] (54) 
V; 


[ol-st 


| 


Figur 10 
Schrittweise Naherungen bei Beriicksichtigung der M -Verteilung. U;,: letzte Naherung nach Figur 9 
mit konstanter Austrittsgeschwindigkeit der emittierten Elektronen; Uj: letzte Naherung mit 
Beriicksichtigung der M-Verteilung; Rr: theoretisch berechneter Potentialverlauf. 
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Die Korrektur wurde erhalten, indem der Speisestrom vor dem Potential- 
minimum um den Faktor vergréssert wurde, der aus der Gleichung (54) her- 
vorgeht. Fiir einen bestimmten Punkt V, muss man die Integration bis U,, 
bzw. V,, ausfiihren. Dann ersetzt V,, die obere Grenze V2 in der Gleichung (54). 
Die schrittweisen Naherungen und die gemessenen Potentialkurven bei Beriick- 
sichtigung der M-Verteilung sind in Figur 10 gezeichnet worden. Die Potential- 
kurve U;,, ist der erzielte Potentialverlauf bei Annahme einer konstanten Aus- 
trittsgeschwindigkeit Uy, und Uy ist die letzte (dritte) Naherung bei Beriick- 
sichtigung der M-Verteilung. Die Methode soll jetzt praktisch bei der dritten 
Naherung erlautert werden. Als Grundlagen der Korrektur dienen die gemes- 
senen Potentialwerte der vorletzten zweiten Naherung. In der Tabelle sind nur 
die Werte von der Kathode bis zum Potentialminimum angegeben: 


xb 0,1 0,2 0,3 0,4 0,5 
U [V] 2148 3,48 4,10 4,12 3,68 
V 0,436 0,696 0,820 0,824 
exp (—V) 0,647 0,499 0,440 0,439 
exp (—V) — exp (—V») 0,208 0,060 0,001 = 
ig (uA) 255 202 184 180 


Aus der Gleichung (50) folgt unter Beriicksichtigung der Gleichung (49) 
1pe= 0 A394, 


Da die vorgegebenen Werte 7, = 0,47, ergeben, ist somit der Anodenstrom 
noch zu gross. Ein weiterer Schritt ist notwendig, um das vorgegebene Ver- 
haltnis 0,4 zu erhalten. 

Aus der vorletzten Zeile der obigen Tabelle fiir den Punkt x/b = 0,1 ersieht 
man, dass nur 0,208 7, als zuriickkehrender Strom erfasst wird. Somit wird die 
Korrektur des Anodenstroms bzw. der zugehérigen Speisestrome nach Glei- 
chung (46) mit F = 0,208 

tq = 15 — 0,208 4, , 


wobei 7) aus dem Potentialverlauf nach der Gleichung (53) berechnet wird. Fiir ° 


das Potentialminimum U,, bei x/b = 0,4 erhalten wir als notwendigen Speise- 
strom 


to = 180 pA. 


Diesen Speisestrom muss man nach Gleichung (54) korrigieren, indem die zu- 
riickkehrenden Elektronen beriicksichtigt werden. Der Betrag dieser Korrektur 
folgt aus der vorletzten Zeile der Tabelle, indem man den betreffenden Faktor 


verdoppelt. Hierdurch ergeben sich die Speisestréme, die in der Tabelle in der 


a on 
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letzten Zeile angegeben worden sind. Der nichste Schritt ergibt nun 
a = 9,409 7, , 


was mit dem vorgegebenen Wert schon gut iibereinstimmt. Auf diese Weise 
wurde der Potentialverlauf Uy erhalten. Die Lage sowie die Spannung des 
Minimums kann man somit sehr gut am Modell bestimmen. Die aus der Modell- 
messung erhaltenen Werte U,, und die berechneten Werte U, des Potential- 
verlaufs sind in der folgenden Tabelle zum Vergleich angegeben worden. 


x|b | 0,1 | 0,2 0,3 | 0,4 | 0,5 | 0,6 | 0,7 | 0,8 | 0,9 | 
(aes — 2,655 | — 3,975 | —4,545 | "4,49 | — 3,9775 | — 3,09 — 1,895 | — 0,43 -f- 1,27 

Uy |—2,32 |—3,72 |\—4,42 |—4,47 —401 —3,13}—1,93 !—0,45| + 1,28 

A% |—12,5 |—6,4 —2,75 |—0,445|]+ 0,8 +13 |+2 + 4 +1 


In der letzten Zeile ist noch der prozentuale Fehler 


ae 100% 

FQ) U, 0 
angegeben worden. Die Abweichung der aus der Modellmessung bestimmten 
Werte von den berechneten Werten ist grésser als in den bisher beschriebenen 
Fallen. Die grésste Abweichung tritt fiir die Punkte vor der Kathode auf. Das 
ist auch physikalisch leicht verstandlich, da wir die Raumladung zwischen der 
Kathode und dem ersten Messpunkt im Modell nicht beriicksichtigt haben. Da 
nach Messungen von A. Hutson [14] die Geschwindigkeitsverteilung nur ange- 
nahert der Maxwellschen entspricht, darf die erreichte Genauigkeit fiir prak- 

tische Zwecke als befriedigend bezeichnet werden. 


4. Kreiszylinderdiode 


Eine Kreiszylinderdiode ist eine kreissymmetrische Anordnung, das Feld 
hangt nur von einer Koordinate ab. In der Widerstandskette werden nicht alle 
Widerstande den gleichen Wert haben, sondern die Werte sind der Reihe nach 
durch die Gleichung (21) gegeben. Damit der Strom durch die Kette klein 
gehalten wird, wurde der erste Widerstand R, = Ry zu 9 kQ gewahlt. Alle 
- Widerstainde wurden auf ihren Nominalwert mit der Genauigkeit ++ 0,3% ein- 
gestellt. Die notwendigen Speisestréme fiir die Nachbildung der Raumladung 
folgen dann aus der Gleichung (24). 


4.1. Austrittsgeschwindigkeit der Elektronen aus der Kathode gleich Null 


Das statische Feld in einer Zylinderdiode ohne Raumladung kann man 
leicht durch Integration der Laplaceschen Gleichung bestimmen. Als Loésung 
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der Gleichung (1) fiir y = 1 und die rechte Seite gleich Null fiir die vorgege- 
benen Randbedingungen U = 0 an der Kathode (ry = 7;,) und U = U,\anider 
Anode (7 = 7,) erhalten wir: 
U, Y 

Als praktisches Beispiel nehmen wir eine Zylinderdiode mit 7, = 0,4 mm und 
v, = 2,4mm sowie mit der Anodenspannung U, = 10 V. Das Verhiltnis 7,/7, 
ist somit gleich 6. Am Modell war 7;, gleich 2 Maschenweiten (2 E) und 7, gleich 
12 E. Da die Maschenweite am Modell zu E = 2 cm gewahlt wurde, ist die Ver- 
grosserung 2 == 100. Die berechneten Werte U, des statischen Potentialverlaufes 


ohne Raumladung nach Gleichung (55) sowie die am Modell gemessenen Werte 


U sind in der nachfolgenden Tabelle zusammengestellt worden : 


i] 


ae Se ee Eee 2 
U, | 2,265| 3,87 | 5,115| 6,13 | 7,0 | 7,73 | 8,40] 8,98 | 9,52 
U. .| 2,28 | 3,85) 5,40. | 603 | 70 | 7,75-) 6 AOumeg ate, ee 


Die Ubereinstimmung ist gut (prozentualer Fehler kleiner als 1%). Das ist fir 
die Bestimmung des Raumladungsfeldes sehr wichtig, da man von dem stati- 
schen Potentialverlauf ausgeht, wie das im Abschnitt 3 erértert wurde. 


Eine exakte Lésung der Gleichung (1) unter Beriicksichtigung der Raum- _ 


ladung fiir den Fall der Zylinderdiode [vy = 1 in Gleichung (1)] und fiir die vor- 
gegebenen Randbedingungen U = 0 beiv=7, und U = U, beir = 7% ist nicht 
bekannt. Eine Naherungslésung wurde von I. LANGMuIR und K. B. BLODGET 
[8] angegeben. Den Potentialverlauf erhalten wir unter Beriicksichtigung der 
Raumladung aus : 


Un ry ie Ue (56) 


Yq Bi 


wobei der Korrekturfaktor 6? eine Funktion 7/7, ist mit dem Anfangswert 
b> =0 bei y= 7, (Erfiillung der vorgegebenen Randbedingung). f? ist der 
Wert des Korrekturfaktors 6? fiir y= 7,. Die fiir die Berechnung wichtigen 


Werte des Korrekturfaktors 6? sind in Figur 11 dargestellt worden. Daraus ent- 


nehmen wir {7 = 0,836. Die Grésse £? nimmt somit in unserem Beispiel. 


Werte zwischen 0 und 0,836 an. Bei den Modellmessungen wird die Raumladung 
durch die Einspeisestréme nachgebildet. Die Gleichung (24) fiir die Bestimmung 
der Nachbildungsstréme kann man in der Form 


1g = Cy 


tie 
Vo (57) 


schreiben, wobei c, eine Kettenkonstante ist, die fiir die gebauten Widerstands- — 
ketten den Wert c, = 3,39 - 102 (V"? cm2) hat. Bei Beriicksichtigung der Ver- 


rats 
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Fiir die Berechnung des Potentialverlaufs in einer Kreiszylinderdiode gilt der Korrekturfaktor 6? 


als Funktion des Verhaltnisses r/r,. 


grosserung ~ im Modell kann man die Raumladungsdichte bzw. die Strom- 
dichte 7;, durch m? dividieren oder aber die Stromdichte j, direkt fiir die Modell- 
abmessungen berechnen. Die Stromdichte 7, erhalten wir aus der Gleichung 


3/2 


jn = 2,34 - 10-8 
Valk 


(58) 


_ Fiir die Stromdichte 7; in der wirklichen Zylinderdiode bei U, = 10 V erhalten 


wir 7, = 9,225 mA/cm? oder fiir die Modellstromdichte 


Tp 12250. A jem 


Fiir die Bestimmung der Speisestréme als Nachbildung der Raumladung er- 


halten wir aus der Gleichung (57) 


(59) 


Nun wurde am Modell das gleiche Verfahren der schrittweisen Naherungen an- 
-gewandt wie bei einer ebenen Diode (Abschnitt 3). Die berechneten Potential- 
_werte U, nach Gleichung (56) sowie die am Modell gemessenen Werte U sind 


in folgender Tabelle zusammengestellt worden: 


-ZAMP VIII/23 


— v. 
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a ee el RRS E i: | 58 | 
U.[v] {1,085 | 2,31 | 3.48] 4.57 | 5,59| 6,56 | 7,48 | 8,37 | 9,20 | 
iv} |1,090 | 2,32 | 3,48| 458 | 5,60 | 6,55] 7,50 | 8,35 | 9,20 | 


# 


Aus den angegebenen Werten U, und U sieht man deutlich, dass die Ab- 
weichung der gemessenen Werte U von den berechneten Werten U, sehr klein 
ist. Der prozentuale Fehler, bezogen auf den genauen Wert im entsprechenden 
Punkt, ist kleiner als 19%. Obwohl die Messungen mit geeichten Universal- 
instrumenten durchgefiihrt wurden, liegt die Ablesegenauigkeit in den erreich- 


1 


—=— S| 


05 


Figur 12 


Potentialverlauf bei schrittweisen Messungen am Modell einer Kreiszylinderdiode r,/7;, = 6 mit 

einer Austrittsgeschwindigkeit der Elektronen aus der Kathode gleich Null. Kurve O: statisches 

Feld ohne Raumladung; I: erster Schritt mit Nachbildung der Raumladung; II: zweiter Schritt; 
Br: berechneter Potentialverlauf nach Gleichung (56). 


ten Fehlergrenzen. Die berechneten bzw. die gemessenen Werte der Potential- 
verteilung im statischen Feld (0) ohne Raumladung sowie die schrittweisen 
Naherungen bei der Nachbildung der Raumladung sind aus der Figur 12 er- 
sichtlich. Nur zwei Schritte waren notwendig, ahnlich wie bei dem ebenen Fall. 


Die zweite Naherung II stimmt schon gut mit dem berechneten Verlauf 
tiberein. . 


4.2. Konstante Austrittsgeschwindigkeit der Elektronen aus der Kathode 


Die Poissonsche Gleichung (1) ist fiir den Fall der Tylinderdiods (» a Ng 


mit konstanter Austrittsgeschwindigkeit der Elektronen aus der Kathode 3 
U, = const analytisch nicht lésbar. Eine Naherungsmethode, die sehr zeit- 


raubend ist, wurde von L. Pace und N. J. Apams [13] angegeben. Der Rech- 
nungsaufwand ist so gross, dass sie praktisch kaum anwendbar ist. Einfacher 


ist der Fall, in welchem eine virtuelle Kathode mit der Spannung U = 0 auf- | 
tritt. Die entsprechende Lésung wurde zuerst von A. VAN DER ZIEL [12] ange- | 
geben. Um so grésser ist hier die Bedeutung des einfachen Messverfahrens an _ 


der Widerstandskette, das in schrittweisen Naherungen leicht durchfithrbar ist. 


a 
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Wir nehmen wieder die Zylinderdiode mit dea geometrischen Abmessungen 
7, = 0,4mm und 7, = 2,4 mm an, wie bei dem Fall der Austrittsgeschwindig- 
keit der Elektronen aus der Kathode gleich Null. Die Anodenspannung wahlen 
wir gleich 10 V und die Kathodenspannung, die der konstanten Austrittsge- 
schwindigkeit der Elektronen aus der Kathode entspricht, gleich U; = 2 V. 
Wenn wir die Stromdichte an der Kathode nach Gleichung (58) mit 744 be- 
zeichnen, so wird die vorhandene Stromdichte j,, unter Beriicksichtigung der 
konstanten Austrittsgeschwindigkeit grésser, wie auch im ebenen Fall (vel. 
Abschnitt 3.2.). Wir kénnen dann schreiben: 


Tk = % Iho» (60) 


wobei x ein Faktor grésser als 1 ist. Nun kann x nicht beliebige Werte anneh- 
men, denn physikalisch ist der Faktor begrenzt, ahnlich wie bei der ebenen 
Diode. Schon bei x = 2,5 wurde der Fall einer virtuellen Kathode mit der 
Spannung Null festgestellt, wobei man die Stromteilung bzw. auch die riick- 
kehrenden Elektronen beriicksichtigen muss. Bei x = 2 tritt aber ein Potential- 
minimum U,, = 1,73 V bei 7/7; = 1,5 auf, wie die nachfolgenden Messungen 
zeigen werden. Mit Beriicksichtigung der Gleichungen (59) und (60) mit x = 2 
werden die notwendigen Einspeisestréme fiir die Nachbildung der Raumladung 


durch 
626 


to= VU MA (61) 
gegeben. Mit den eingestellten Randbedingungen am Modell (U, — U; = 8 V) 
wurde zuerst das statische Feld ohne Raumladung gemessen. Die gemessenen 
Werte wurden mit den nach Gleichung (55) berechneten Werten verglichen. 
Die Abweichungen sind kleiner als 1%, was wir als befriedigend annehmen 
konnen. Ausgehend von dem statischen Feld ohne Raumladung wurden die 
notwendigen Einspeisestréme als Nachbildung der Raumladung nach Glei- 
chung (61) bestimmt und die Einspeisestréme im Modell eingestellt. Es waren 
vier Schritte notwendig. Nach dem vierten Schritt war die Anderung im Po- 
tentialverlauf so klein, dass praktisch keine Korrektur der Einspeisestréme 
mehr méglich war. Die gemessenen Werte U des Potentialverlaufs nach dem 
letzten, vierten Schritt sind in der nachfolgenden Tabelle zusammengestellt 


worden: 


rite | 15 | 2 | 25 | 3 [ aat | 4 [eb 5 J 55 | 
Ti | 1,73 | 2,12 | 2,85 | 3,76 | 4,75 | 5,80 | 6,83 | 7,95 | 8,95 


Die Ergebnisse der schrittweisen Naherungen bei der Bestimmung des Poten- 
tialverlaufs sind in der Figur 13 gezeichnet worden. Mit (0) wurde der statische 
Potentialverlauf ohne Raumladung bezeichnet und mit I bis IV der Potential- 
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verlaut nach schrittweisen Naherungen. Aus dem gezeichneten Verlauf kann 
man leicht die Stelle sowie auch den Spannungsbetrag im Potentialminimum 
ermitteln. Da die theoretische Berechnung kaum ausfiihrbar ist, bilden die 
einfach durchfiihrbaren Modellmessungen bis jetzt die einzige praktische 
Methode zur Lésung dieses Problems mit einer konstanten vorgegebenen Aus- 
trittsgeschwindigkeit der Elektronen aus der Kathode bei Kreiszylinderdioden. 


Figur 13 
Potentialverlauf bei schrittweisen Naherungen in einer Kreiszylinderdiode r,/7,, = 6 mit konstanter 
Austrittsgeschwindigkeit der Elektronen aus der Kathode U;, = 0,2 U,. Kurve O: statisches Feld 
ohne Raumladung; I: erste Naherung... usw. IV: letzte (vierte) Naheruns bei Nachbildung der J 
Raumladung. 


5. Kugelsymmetrisches Elektrodensystem 


Praktisch haben die kugelsymmetrischen Elektrodensysteme eine viel 
kleinere Bedeutung als die kreissymmetrischen (Zylinderdiode). Da sie doch 
auftreten kénnen und auch die theoretische Naherungslésung nach I. LANG- 
MUIR und K. B. BLopcEt [9] bekannt ist, werden wir kurz ein solches Elektro- 
densystem behandeln. Die Messungen sind einfach und schnell ausfiihrbar. All- 
gemein kénnte man mit dem beschriebenen Messverfahren analoge Probleme 
fiir den -dimensionalen Raum lésen. Da sie aber praktisch keine Bedeutung 
haben, beschranken wir uns auf den kugelsymmetrischen Fall » = 2 in Glei- 
chung (1). 

Als Beispiel nehmen wir ein Elektrodensystem an mit Vn = 0,2 mm und 
7,=1,6mm. Die Anodenspannung gegen die Kathode (U,=0V) betragt 
U,=10V. Nach den Ausfithrungen im Abschnitt 2 hat der sae Widerstand 
in der Kette den Wert R, = co. Da wir am Modell erst den nachsten Punkt 
vom Kugelzentrum aus als den Kathodenpunkt wéhlen, stért das die Mes- 
sungen nicht. Damit der Strom in der Kette klein bleibt, wurde der Wider- 
stand R; = Ry = 9kQ gewahlt. Die weiteren Widerstandswerte in der Kette 
folgen aus der Gleichung (28). Alle Widerstande wurden mit der Genauigkeit 
(Toleranz) + 0,3°% auf den nach Gleichung (28) bestimmten Wert eingestellt. 


ie 
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Hierdurch wird die erzielte Genauigkeit bei den Messungen des Potentialver- 
laufs befriedigend. Die fiir die Nachbildung der Raumladung notwendigen 
Einspeisestr6me werden nach Gleichung (31) ermittelt. Im Modell wurde an 
der Widerstandskette r, mit einer Maschenweite nachgebildet. Da das Verhilt- 
nis 7,/7;, = 8 ist, waren fiir die Nachbildung der Anodenhalbmesser 7, acht 
Maschenweiten (8 £) notwendig. 

Die Lésung der Laplaceschen Gleichung [rechte Seite der Gleichung (1) 
gleich Null] fiir den kugelsymmetrischen Fall y = 2 ergibt 


eg eid ae, Oke a ( 
IGA Me |. (62) 
Dabei wurden die Randbedingungen U;,, = 0 bei r= 7, und U = U, bei r=% 
beriicksichtigt. Die nach Gleichung (62) berechneten Potentialwerte U, und 
die am Modell gemessenen Werte U sind zum Vergleich in folgender Tabelle 
zusammengestellt worden: 


a RE Og ee a a en ee ae 
U, 5,715 | 7,62 8,575 | 9,15 9,52 9,80 
U 5,7 7,60 8,55 9,15 9,50 9,80 


Ein Vergleich der berechneten Werte U, und der gemessenen Werte U des 
statischen Potentialverlaufs ohne Raumladung zeigt, dass die Ubereinstim- 
mung gut ist. Der auftretende relative Fehler ist kleiner als 0,3%. 
Die Gleichung (31) fiir die Einspeisestréme als Nachbildung der Raumla- 
dung kann man analog zur Gleichung (57) in der Form schreiben 
to os Cy ve , (63) 
wobei die Kettenkonstante c, in unserem Fall den Wert c, = 42,35 (V‘? cm?) 
hat. Die Stromdichte 7, wird jetzt 


U2? 
jx = 2,34 - 10-86 *,, (64) 


2 
a’, 


-wobei «2 der Wert des Korrekturfaktors «? [vgl. Gleichung (66)] bei r= 7, 
ist. Bei der Bestimmung der Einspeisestréme nach Gleichung (63) kann man 
die Vergrésserung n = 100 im Modell beriicksichtigen oder die Kathodenstrom- 
dichte direkt fiir die Modellabmessungen bestimmen. Fiir die Stromdichte 4; 17 
. am Modell erhalten wir fe bigts gah lone 


‘und nach Gleichung (63) fiir die Einspeisestrome 


- — 9500 65 
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Betrachten wir nur den Fall der Austrittsgeschwindigkeit der Elektronen aus 
der Kathode gleich Null. Die Lésung der Poissonschen Gleichung (1) fiir den 
kugelsymmetrischen Fall y= 2 in Gleichung (1) nach I. Lancmurr und 
K. C. BLopcGEtt [9] wird 
2 

oe (Se) rig (66) 
Der Korrekturfaktor «2 ist eine Funktion von 7/7,, ahnlich wie bei der Zylinder- 
diode f2. Er erfiillt die Randbedingungen an der Kathode und an der Anode. 
Mit «2 ist der Wert des Korrekturfaktors a? bei 7 = 7% bezeichnet worden. Die 
nach Gleichung (66) berechneten Werte U, des Potentialverlaufs mit Beriick- 
sichtigung der Raumladung sowie die am Modell nach dem letzten Schritt mit 
schrittweisen Naherungen gemessenen Werte U sind in nachfolgender Tabelle » 
angegeben worden: 


rl, | 2 | 3 | 4 | 5 | 6 | 7 | 
Cavite 3,50 5,67 7,0 8,07 8,85 9,50 | 
U [Vv] 3,60 5,65 7,0 8,05 8,85 9,50 


Die Ubereinstimmung ist befriedigend, die relative Abweichung ist kleiner als 
1%. Der statische Potentialverlauf (0) ohne Raumladung sowie die schritt- 
weisen Naherungen bei Beriicksichtigung der Raumladung sind in der Figur 14 


ue 
Ug 
05 


Ef 


Figur 14 


Potentialverlauf bei schrittweisen Messungen in einem kugelsymmetrischen Elektrodensystem 
Y al" = 8 mit einer Austrittsgeschwindigkeit der Elektronen aus der Kathode gleich Null. Kurve 0: 
statischer Verlauf ohne Raumladung; I: erste Naherung bei Nachbildung der Raumladung ... usw. 

III: letzter, dritter Schritt, der mit dem berechneten Verlauf Br gut tibereinstimmt. . 


' 


veranschaulicht worden. Hier ist der Einfluss der Raumladung stark. Es 
waren drei Schritte notwendig. Erst nach dem dritten Schritt war die Anderung 
im Potentialverlauf vernachlassigbar klein. Die dritte Naherung III stimmt gut 
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mit dem berechneten Verlauf By iiberein. Ganz analog wie im Beispiel der 
Zylinderdiode kénnte man auch hier den Fall der konstanten Austrittsge- 
schwindigkeit der Elektronen aus der Kathode behandeln. Da er aber praktisch 
nicht wichtig ist und auch bei der Anwendung der Modellmessungen mit 


Widerstandsketten nichts Neues ergibt, werden wir diesen Fall hier nicht weiter 
betrachten. 
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Summary 


Resistance chains, used as analogy computers, allow the determination of 


be 


static potentials without and with space charge of high vacuum diodes, the fields 
of which may be considered as onedimensional. In the case with space charge, 


A viii a 4 fs 2 at 
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additional currents, fed into the resistance chain, are essential for the represen- 
tation of the space charge. Such resistance chains with additional current sources 
afford solutions of Poisson’s equation 

OFWs 2 gpd Ves ted @ 


Or? y Ov ie 


The first problem under consideration is the plane diode, for which »y = 0. Here 
the solution of Porsson’s equation is obtained by iteration: 

(a) with zero emission velocity at the cathode; 

(b) with finite emission velocity ; 

(c) with Maxwellian emission velocity distribution at the cathode. 

The solutions take much less time using the resistance chain than using the 
well-known formulas for these problems. The second problem concerns the cir- 
cular cylindrical diode (vy = 1) with zero and with finite emission velocity at the 
cathode. In the latter case, no straight-forward mathematical solution is known. 
Again, practical solutions are easily obtained by iteration using the resistance 
chain. Finally, a spherical diode (vy = 2) is considered at zero emission velocity. 
Here also, the chain affords easy solutions. The computed chain solutions were 
compared with known formulas in several cases, showing discrepancies not over 
1% im most Cases. 


(Eingegangen: 19. Februar 1957.) 


Plastic Yielding of Notched Strips Under Conditions 
of Plane Stress 


By Hueu Forp and GEorGE Lianis, London?) 


Introduction 


The problem of thin notched strips subjected to tension has been studied by 
HILt [2]?) for the case where the notches are symmetrically placed on each side 
of the strip. The deformation corresponds to plane stress and, using the theory of 
characteristics, HILL finds the upper bound solutions for the two cases of 
V-notched and circularly notched strips. BisHop [1] has shown that HILt’s 
upper bound is the true solution for the V-notched case, but gives no proof 
for the circularly notched case. 1 

In this paper, a possible statically admissible stress field (lower bound) is 
found for these two cases, but the main purpose is to examine single notched 
strips subjected to bending in the plane of the sheet. Upper bounds are found for 


single V- and circularly notched strips in bending and the general lower bound 
solution is applied to these two cases. 


. 
1) Imperial College of Science and Technology (University of London) 
Mechanical Engineering. 


2) Numbers in brackets refer to References, page 382. 


Department of 
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Statically Admissible Stress Field for Plane Stress — Simple Tension 


The symmetrical field shown in Figure 1 is a statically admissible field if it 
satisfies the following conditions: 

(1) The boundary stress conditions. This is satisfied for the directions of 
the stresses shown in the figure. 

(2) The equilibrium equations. Only discontinuities of the direct stresses 
normal to the surfaces AB, CD, AO, BO, CO, and DO, are allowed. Since these 
conditions are satisfied automatically across AB and CD, only the discontinui- 


Figure 1 


Stress trapezoid. 


ties AO and CO need to be examined. These provide four equations of equi- 
librium. i 
(3) The yield criterion must nowhere be violated. This condition provides 
three inequalities for the regions AOB, COD, and AOC. is 
The above four equations and three inequalities are sufficient to evaluate 
the seven unknowns 04, 0», /, $1, Sz, y, and 6, in terms of the angle f. In addition, 
the condition that the load acting on AB must be a maximum provides another 


equation to evaluate f. 


>. 


le Ee ee hd > ee ir 
aNd “a. 


From condition 2 


From the continuity of the normal (c) and shear (t) stress across AO it 


follows: 
o = 0,sin*y + o, cos*y = fsin®(B + y) , (1) 


1 = (0, — o) siny cosy = f sin(B + y) cos(B +7) - (2) 
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Therefore sin(B + y) cosB 
=} siny ; 8) 

sin(B + y) sin B 
o,=} cosy, 4) 


Also from the continuity across CO it follows: 


o’ = s, sin?6 + s,cos*6 = fsin?(d — 8), (5) 
1’ = (s; —S) sind cosé = f sin (6 — B) cos(6 — B) . (6) 
Therefore 
a ey a 
Spee sin ( met cos B (7) 


sin(é — 6) sin B 
cos 0 ‘ 


Sg=f 


From condition 3 


The von-Mises yield criterion for plane stress in terms of the principal 
stresses 1s: 
2 2, 2 
0% + OF — 9, 6, & be 


where Y is the yield stress in simple tension. Therefore: 
For the region AOB: 


of + 02 — 016, & Y?. (9} 
For the region AOC: 
ORS (10) 
For the region COD: 


From the equations (3), (4), (7), and (8) and the inequalities (9), (10), and (11), 
the seven unknowns can be evaluated. © 

Substituting o,, og, from equations (3) and (4) into (9) and s,, ss, from 
equations (7), (8), into (11) it follows: 


Ses rie al cos? B sin? B sin? B We q 
Calta ah a cos?y a = pe’ (12) 
aaa Garces cos? B sin? B sin? B YE. | 
eae re cos? 6 aad = hoe (13) 


From (12) and (13) y and 6 can be found in terms of Y/f.. 
The load acting across AB is a maximum if o, is a maximum. If the region 
AOB is fully plastic from von Mises’ ellipse (Figure 2) it follows that maximum 
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Lf) 


Figure 2 
Yield locus for plane stress. 


o, corresponds to the point P: 


2 ae 
— ye ) a ee 
O71 3 O2 V3 (14) 
From (3), (4) and (14), it follows 
Cot Pp cote 2, (15) 
sin(P + y) cosB 2 We 
siny in V3 , , (16) 
or 
TPSy, 
2 = Ve 1. SS a eee 1 
cost = 2 ( 2 r) (17) 
From (10) and (17) it follows 
cos?B 2 (1-7) or B= 23°. 


For B = 23°, f = Y and the region AOC is fully plastic. 

For £ = 23° inequality (13) is satisfied for 6 S 61° 15’. For 6 = 61° 15’, we 
find s, = 0-650 Y, sz = — 0-505 Y. The four regions of the trapezoid are ina 
fully plastic state if 


2 1 J reed 
eh ee ee ae oint P in Figure 2) | 
O71 V3 2 V3 (p § ) | 
5,=0650Y, s,=—0-505Y, | (18) 
B= 23°, y= 49°40’, 6 = 61°15’. 
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V-Notched Strips Loaded in Tension 


For « <2a/2— B < 67° the fully plastic trapezoid created from equation 
(18) can be inscribed into the strip without cutting the boundaries (Figure 3). 


Figure 3 


V-notched strip in tension. 
In this case the constraint factor is constant. Ly = 2IV3 = 1-155 for0 =a 6/e 
This is a lower bound for L. HiLv’s upper bound for the constraint factor is — 


Ly= Eos == 15155) dori" 102 | 


) (19) 
ae V(1 + 3 sin? @) + 3 cosy 


23 
where y is found from the equation 2 tang + tan2 (« + 2 p’ — gy) = 0 and 
B’ = 54° 44". Therefore the upper and lower bound coincide for 0 < a < 67°” 
which gives the collapse load of the V-notched strip. 


For a 2 67° it is taken 6 = 2/2 — «, in other words the sides of the notch 
and the trapezoid coincide. In this case: | 


for 70° 3273 = 00" | 


Bee cos (y = sin Y, (a) ce cos (y a cos « Y, (b) 
cos(é + «) sing cos(d + a) cosa (20) 
3 as Yo we) So. = ¥ sod) 


cos 6 
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For the constraint factor to be maximum it follows from equation (20a) 
that y must be the minimum possible and it is- defined from inequality (13). 
The latter for 8 = m/2 — « and f = Y becomes 


(21) 


a SIN? oe COS* a sin? « 
wee &) sin?» “ty cos*y sin? y =it 

The following table gives the optimum value for y and the lower bound for 
the constraint factor: 


Table 1 
ee ee ee ee ee 
lo | 70° | Too | 80° | 85° | 90° 


50° 30’ SS FO W522. 30. 54° = 
20° a5? 10° an Oe 


47° as 
1-060 | 1-000 


Om if e7°a0".|- 38°30 


¥ 

B 

B. 1-152 1-132 1-103 1-058 1-000 
Ne 

L 1-155 1-141 1-105 


For the region 67° < « < 90° the two bounds almost coincide (Figure 4). 
This agrees with BisHop’s proof [1] that the upper bound is the true solution. 
BisHoP, however, showed this to be so for values of « between 0° and 90°; the 
above limit analysis gives the true value of the constraint factor of 1-155 for 
0 <a < 67°, but a value slightly below the true value from 67 < a < 90°. 


HILL'S upper bound | 


=e 


1-15 


lower bound 


— 
eck 
(Ter 


Constraint factor —— 
3 
—+— 


l= 


e Cems cm me oCume0cr = h00 es WANl SSNs 2025 10° 0° 
o=Notch angle —=— 
Figure 4 
Double V-notched strip under tension. 
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Circularly Notched Strips in Tension 


For double-circular notches, radius 9, with the stress trapezoid drawn 
tangential to the notch surface (Figure 5), 


Figure 5 
Circularly notched strip in tension. 


oa1- 48 (==). (22) 


a a cos B 


Therefore, the lower bound for the constraint factor is given by 


Se; 20 (1—cosB 
te Y [2 a: ( cos B )| : ae 
The angle f is defined so as to maximize qs Reo | 
0(0,/Y) 2 05/1 = cosh Ope MAN ih bos 4 
op [1 a ( cos B )| Via s-coste os et) 
where 
0(0,/Y) Oy 
AE and y 


d 
are defined from equation (3) with f = Y for 0° < 6 < 23° and from equation 


(18) for 6 = 23°. Solving numerically equation (23) and (24) the values ae 
in Table 2 are found. 
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The upper bound solution found by Hit [2] is given in the following 
approximate formulae: 


2 


4 


Pomel ot 205 tora nbeee eee Le. 0S "= 0-664. 
a-+ a ra) 


a+ 


L,=—-— 0080-2 for0<%<0-465, ie, 0664 —*  <10. 
V3 a a a+ o- 

We can see from Table 2 and Figure 6 that the mean value differs by less 

than -+ 0.9% from the bounds. BisHop did not prove that HILv’s upper bound 


is a true solution. The present limit analysis gives a lower limit close to the 
upper bound and indicates that HILv’s solution is the true solution. 


Table 2 


L, | 1-155 | 1-133 | 1-118 | 1-102 | 1-087 | 1-072 | 1-057 1040 | 1.025 1-012 
1-155 | 1-146 | 1-135 | 1-120 | 1-096 | 1-075 | 1-057 i ieee 1-012 


Constraint factor —— 


tes 


Figure 6 
Circularly notched strip under tension. 
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Statically Admissible Stress Field for Single-Notched Strips Under- 
going Plastic Bending Under Conditions of Plane Stress 


Let a horizontal band stressed uniformly by the stress s be added between B 
and D (Figure 7). The state of stress in the regions I, II, and III, is shown in 
Figures 8a, b, and c. The von-Mises criterion: 


O20 6 6 31 ee 
applied in the three regions gives 
—Y<s<=+/Y (from Figure 8a) , 
(a) | v (=> + 3 cos?B — 1) <0 (from Figure 8b) , 
| s <0 [since sz < 0, equation (8)] (from Figure 8c) . 


The above inequalities are satisfied for 


RNS 


(02s 23 cos? =1)) == for S51) ae 
/3 
or 
Dy pied ees Vee a for 23° < Bp < 35°10’ 
or 


Sa SY op sin (Frey) 008 et for 0 = 5 S223"3 
siny 


(c) If we keep therefore 6B = 23° and s=—Y then o,=const = 2 Y/V3, 
whilst for 0° < B S 23°, s = —Y, 9, is given by equation (3). 


The dimensions of the field are found from the geometrical data of the 
strip and the notch, the requirement that the oblique side of the trapezoid 
must be tangential to the notch, and the equilibrium equation along the 


minimum cross section of the strip. The latter gives (Figure 7) d = c o,/Y and- | 


therefore 


2 
ae oor 423° = 6 355ee (25) 


_ , Sin(B + y) cosB ° ° 
dee a for O° S08 = 23>, 
where y as a function of 8 is given by (12) with the equality sign and f= Y, 
and o, from equation (3) or (18). 
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| 


| Notch surface _. 


V 


— a——> 


<— § — > 


Figure 7 
Stress field for single-notched strips undergoing plastic bending 


ZY, ~ 
: 


Ysin?a 


ee ae Ta —>[¥cos?B +5] 


YsinB cosB Y 
ly 


S2 


a 


52 
Figure 8 
Elements showing states of stress in regions I, IT, and III in Figure 7. 
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The constraint factor is 
Lie aii , (26) 
and from Figure 7 it follows 
pata a 
Therefore 
1,-42 (2 +94). (27) 


Case 1: Plastic Bending of Single V-Notched Strips 
(A) Upper Bound 
Solution 1 


The slip-line field is shown in Figure 9. In OAB the field consists of 
straight lines inclined to the stress free surface at an angle 


p= tales V2 = 54° 44’ 


(Hix [2], p. 28). The state of stress is uniaxial tension Y = V3k parallel to 
the edge. H1Lt [2] has shown that the stress field around the singularity is: 


6,=kcosp, o,=2kcosy, T,=ksing, (28) 


Figure 9 
Characteristics for V-notched strips in bending. Solution 1. 
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where the polar angle is measured from the base line which is inclined at the 
angle 2 8 to the edge AB. The one family of the characteristics are straight 
lines while the other family, such as BC, have polar equation: 


vy sing = const . (29) 


The field to the left of C consists of parallel straight lines and the state of 
stress is biaxial tension parallel to the axes. In order that there should be no— 
shear stress along ON the angle q, is defined from the equation [HILt [2], 
equation (14)] 

2 tang, + tan2 (2 + 28 — g,) =0. (30) 


The angle CNH must be greater than the angle defined by a radial field 
along the end radii of which the states of stress are those in the stress discon- 
tinuity point N in the regions I and IT. Otherwise the yield criterion is violated 
in some parts of the rigid material near the corner N. The analogy of HENcKY’s 
first equation for plane stress is [HILL [3], equation (44), p. 302]: 


d(A — w) = 0 along «-lines, (a) a 
d(A + w) = 0 along f-lines, (b) | 
where 
A= tant (2tanP)—+-, a) | 
(32) 


vy Es We 
O= T+ t+ zat e% 2° (b) 


Pa Pg are anticlockwise angles of the characteristics from a fixed direction x, 
and the angle V’ is defined from the ratio of principal stresses 0, 0: 


fy | On KOs 33) 
Ww = sin eae | (33) 
From equations (31a) and (32b) it follows that the angle of rotation of NC to 
meet NH for a radial field consistent with the state of stresses in I and II is 
equal to 


Py Py 
Ay — Ag + LG 
Therefore e zt 
I II 


_ From Figure 12 it follows 


ot Uv 
Ag=a+3B-m- >> Y= 2 2B — Pe) — 3.2 Pay ae rs 
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Substituting in (34) and by means of equations (30) and (32a) we find: 


[tan-1(cot m,) — tan! (2 cot2 f)] Sa +3 B— %— 7 
or 


[F — gy, — tan-1(2 cot2 p)| <a4+3B—%- ; 


or 
a 2a—3 8 -—tan-1(2 cot2 f). 


But since tan B = V2, tan-1(2 cot2 £) =2/2— 8, 
O23. ie Ai On glae S1e. (36) 


By means of the angles f, yp) and equation (29) the field is fixed if c is known. 
c can be found from the equation of equilibrium. The stress normal to ON is 
[Hit [3], equation (41)] 

k (2 sin ¥, ab 1) 


“Vit 3sin’y, 


and because of equation (34) 


Y [3 cos2 («0 + 2 B— mo) + 1] 


m . ey 
V3 V1 + 3 cos?2 (x + 2 B — @) Se, 
Therefore 
COm = Y (a—c) (38) — 
and 
es 30 
a ny [3 cos2 (a + 2 B — mo) + 1] : (39) 
V/3/1 + 3 cos®2 (a + 2 B — qm) 
where @, is found from equation (30). 
The yield moment is 
WG, 6 ee (ae)? 
Hs Zt By 
and because of (38) 
oe Cin 
M, == o»,(1 “ 4 
Therefore the constraint factor is (Figure 12) 
vr 2 Ory, 
4 
_ WY ye 
ee (40) 
(1+) 


This solution however holds as long as Yo = 9 or from equation (30) for 
ee 10° 32h 


Vol. VIII, 1957 Plastic Yielding of Notched Strips Under Conditions of Plane Stress S73 


Solution 2 


For a= 70° 32! the field inside OA BCN (Figure 9) shrinks to a coincident 
pair of characteristics along the transverse axis (Figure 10) and it retains its 
shape for 

ee es fen 
ON is therefore a local neck. 


Figure 10 


Characteristics for V-notched strips in bending. Solution 2. 


The deformation mode along the neck is shown in Figure 10. There is a 
discontinuity of the normal component of the velocity along the neck of 


amount 
An] = 2% 


Along the neck Y = z/2 and the principal stresses and strain rates are 


I 


ey . (normal to the neck) , 
é,=0 (along the neck) , 


é,= = (normal to the strip plane) , 
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where b is the local breadth of the neck. 

" 2 fe SA Cc 1 
thee a a ee ea a ey 
V3 Valea wpe em 


V3 


in which C = ON. The constraint factor for 0 S « < 70° 32’ has the constant 


? 


value (Figure 12) 
2 2 2 
[en ee ate ig (1 i ca 21-072. 
(: ae reah V3 V3 
y3 


The field remains unaltered for this region. In Table 3 the sizes of the field, 
the maximum stress, and the upper bound for the constraint factor are given. 


(B) Lower Bound 


Figure 11 shows the statically admissible stress field discussed previously 
fitted to V-notched strips. 


(a) For 0 Sa < 67°, B = 23° = const: 


Figure 11 


Lower bound for V-notched strips in pending. 


vt, Be a Te Wo ee 
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and therefore 1 


; (1 + i i= = 1-072 = const. 


375 


We see that for 0° < « < 67° both upper and lower bounds coincide the value 


1-072 being therefore the true constraint factor. 


Table 3 
a 90° | 85° | go° | 15° | 70° 67° OF ess OF 
Yo 54° 447 47° 38°:30/'| 27° 30’ 0° 0° | 
pe Be 1-0 1-053 1-095 1-124 1-155 1155 1-155 
cla 0-50 0-486 0-477 0-470 0.464 0-464 0-464 
nae 1-00 1-025 1-045 1-060 1-070 1-072 1-072 
B 0° 5° 10° 15° 20° 25° 23° 
ayy 1-00 1-053 1-102 1-132 1-150 1-155 IS) 
L, 1-00 1-025 1-045 1-060 1-070 1-072 1-072 


K upper bound 


1:06 TH lower bound 
pal = 
-> 


a 
Ke 
tt ese | 


L 


Constraint factor —— 
ea ies 
Mm Ww 
Bees 


Figure 12 
Single V-notched strip in bending. 
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(b) -Bor67° Sa = 90": 
We take B = 2/2 — «. The constraint factor is given by 


i an rae ae (41) 

(1+ -¥) 
Using equations (12) and (21a) o,/Y in terms of B can be found. Substituting 
this value in equation (41) the lower bound for the constraint factor is calcu- 
lated (Table 3 and Figure 12). We can see that the difference is almost zero, 


and the two bounds practically coincide for this region as well. The maximum 
stress along the notch axis also is almost the same in both bounds. 


Case 2: Plastic Bending of Circularly Notched Strips 
(A) Upper Bound 
Solution 1 


The field of characteristics in the upper region is radially symmetrical, as 
in a hole expanded by internal pressure (Figure 13). It is defined by the para- 


Figure 13 
Characteristics for circularly notched strip in bending. Solution 1. 


_—— 
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metric equations: 


oS 2h sin (6 2: =) alee as ¥. sin (0 - 2) CGV eV8(0—a10) 


where 6 = tan-1(/3 sin). 

The field is drawn as follows. Let us take the characteristic AS: 

sinceo, S o,,@4=« and sinceo,=0, 6,=2/6, VY, = 19° 28’ = 0-340, 
Ag (a/6) = 0-447 (Hix [3], p. 350, Table 2). Also wy = « — gp. Along the line 
AS it is: w + A=const. Therefore 


Aq(0) +a — p = 0-447 + & 


or 


Q= sin-*( 7 sin 8) — = sin-t fa tan a) — 0-447. (43) 


Equations (42) and (43) are the equations of the characteristic AS in para- 
metric form. 


The upper and lower field extend until they meet in N. 
From the equilibrium equation it follows: 


o+c o+c 


2 . d | V3 ee 
(a — c) ¥ = [oy dr = 2h @ [ sin(9 + a) oy | seco | do, 
: Q @ : 
and therefore 
i 7 2 1+c/o 
ide = ya [sin(0+ <) a(=). (44) 


i 


By means of equation (44), c/a in terms of a/(a+ @) is calculated and 
plotted in Figure 14. 

As was sought previously, the angle of rotation Ap of the characteristic ZN 
(in order that HN can meet ZN) must be greater than 


(Ay): — (Ay) + if. 7 ee) ' 


But from Figure 13 it follows 


‘ Therefore 
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Coe 


Figure 14 
Curve giving the solution of equation (44). 


But it is always A < /4 so this condition is fulfilled until a neck is formed. 
The neck is formed when the characteristics coincide on the axis. In this case: 


Cc 


el % _ 4.265, ie = 07s 


ies pan? eA) om il 
2° 5 ates 2 ae 


Substituting in equation (44) we find 


a 


(*) eae BE 


This solution is therefore valid for 


ee ean. 
0S so 50002, (45). 


CONSTRAINT FACTOR 


Taking the moments of the normal stresses across the axis we find 


6 a 


i 


[ae 4c ff sin (0 + 2/6) cos(6 — 2/6) 2V3(0—m/6) 40 


cos? 


6 aN He 1+ c/a 
i (1 4 ie assy ). 
The integral is evaluated numerically. For a certain value of al(a + 0), cla 


and therefore c/o is found from Figure 14. Substituting in equation (45) we 
calculate L, (Table 4 and Figure 17) 
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Solution 2 


For 0-692 < a/(a+ 0) <1 the field of characteristics extends over an 
angular span of 38° 56’. A neck is formed joining the corners of the upper and 
lower domains (Figure 15). Normal to the neck a constant tensile stress is 


Figure 15 


Characteristics for circularly notched strips in bending. Solution 2. 


applied equal to 2V3 Y. The rigid parts rotate around N. Again along the 
neck there is a discontinuity of the normal velocity of amount: 


|Av| = 2. 


From the equilibrium condition it follows 


2-0710 P 
ames faye Ae amc ae 
V3 | 


- But from equation (42) and Figure 15 it follows that 


2-071 


[ovdr=1172 V0. 
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Therefore 
fa AAO Olt 0-464 :04 oe (47) 


Having c, and c, the sizes of the field are fixed. 


CONSTRAINT FACTOR 


The point of application of the resultant force along ON has a distance 
from K equal to 


2.0710 
[6 vy dv 
1:827 02 Y 
° els 
+= “3.0710 PijzZoy ane 
| % dy 


@ 


By taking the moments with respect to N we find 


2 (oe 1 
MM, = 3112.0 ¥ 10, + OSE 0) + /3 — ee Ne 3 (a 1-071 o+-.6.)*- 


By means of equation (47) we find 
2 
L, = 1072 — 0-123 = + 0-022 = 


2” 


0:692 = See = 48 
— a ka Q =—=— 2 ( ) . 
Values of iz for this region have been plotted in Figure Ws 


(B) Lower Bound 


_, The. statically admissible stress field fitted for circular notches is shown | 
in Figure 16. From equilibrium it follows 


(1 — cos B) a Ce eS o/a (1 — cos f)/cos B 


ie aes @ 
y and d+c+ ee ka isiv eu 


Substituting in equation (27) we find the formulae for the constraint factor 7 


a: @ (1—cosf) ]2 0,/Y . 
Ly 2]: eo | earae. (49) 


B is defined in terms of @/a so as to maximize L,, and o,/Y(f) is found from | 
equation (3) or (18). In Table 4 and Figure 17 both bounds are shown. We can . 
see that they coincide in the greater part, their greatest difference from a _ 
mean value being less than + 0:3%. 
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Lower bound for circularly notched strips in bending. 


Figure 16 


{ 


upper bo 


| 


Figure 17 
Circularly notched strip in bending. 
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Table 4 
alia+e] 0 | o-1 | o-2 | 0-3 | o-4 | 0-5 | 0-6 | o-7 | 0-8 | 0-9 | 1-0 
L, | 100 | 1-001 |1-0025 | 1-005 | 1.008 | 1-011 | 1-014 | 1-022 | 1-0424| 1.059 | 1.072 
L, | 1-00 | 1-001 | 1-025] 1.005 | 1-008 | 1-011 | 1-014 | 1-022 |1.038 | 1-053 | 1.072 
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Zusammenfassung 


Das Problem diinner gekerbter Flachstabe, die unter Zug stehen, ist bereits 
von Hitt [2] fiir den Fall, wo die Kerben an jeder Seite des Flachstabes sym- 
metrisch angeordnet sind, untersucht worden. Die Verformung kommt dem 
ebenen Spannungszustand gleich, und HILL vermag mit Hilfe der Charakteristi- 
kentheorie eine obere Grenzlésung fiir die zwei Falle von V-formig und kreis- 
formig gekerbten Flachstaben zu geben. Spater konnte BiscHop [1] zeigen, dass 
Hitts obere Grenze die wahre Lésung fiir den V-f6rmig gekerbten Fall darstellt, 
bleibt aber fiir den kreisf6rmig gekerbten Fall den Beweis dafiir schuldig. 

In der vorliegenden Arbeit wird fiir beide Falle ein statisch zulassiges mdg- 
liches Spannungsfeld (untere Grenze) gefunden; dennoch hegt die Hauptaufgabe 
in der Behandlung einfach gekerbter Flachstabe bei Biegung in der Flachebene. — 
Obere Grenzen konnen bei Biegung sowohl fiir einfach V-f6rmig wie fiir kreis- — 
formig gekerbte Flachstabe angegeben werden, doch ist es méglich, auch eine all- 
gemeine untere Grenzlésung fiir beide Falle zur Darstellung zu bringen. 


(Received: March 22, 1957.) 


Eine Stibchenmethode zur DK-Messung, 
Anwendung auf (NH;),H3JOs¢ bei 3 cm Wellenlinge 


Von HEINI GRANICHER?!) und WERNER SCHURTER?), Ziirich 


1. Einleitung 


Die Untersuchung der Frequenzabhangigkeit der dielektrischen Eigen- | 
schaften von Isolatoren ist nicht nur technisch, sondern auch wissenschaftlich : 
von grossem Interesse. Die elektrische Polarisation, die in einem Ionenkristall | 


beim Anlegen eines E-Feldes aufgebaut wird, setzt sich aus einem Elektronen- 
‘beitrag (Verschiebung der Elektronenhiille relativ zum Atomkern) und einem 


1) Physilcalisches Institut der ETH. 
2) Physikalisches Institut der ETH; jetzt: Imperial Chemical Industries Ltd., Billingham 


England. 
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Ionenbeitrag (Verschiebung der verschieden geladenen Ionen relativ zueinan- 
der) zusammen. Die beiden Anteile unterscheiden sich in ihrer Frequenzabhan- 
gigkeit: Die Elektronenpolarisation kann Wechselfeldern mit Frequenzen bis 
etwa 101° Hz folgen. Der Ionenanteil-dagegen ist nur bis zu den Gitterfrequen- 
zen im infraroten Spektrum (etwa 1012 Hz) wirksam. Molekiilkristalle, deren 
Molekiile ein permanentes Dipolmoment besitzen, weisen einen zusatzlichen, 
in der Regel grossen Polarisationsbeitrag auf, der von der Ausrichtung der 


BaTiO, [ferroe/ektrisch] 


1000 


SrTi0; [ paraelektrisch] 
-—_NHa)2H3J0g CaTi03[parae/ektrisch] 


om, 


° 


100 “Ti09 (Rutil) 


10° 10° (ap se 


Figur 1 

Dispersion der DK von Titanaten bei 25° C [2] und von (NH,),H,JOg: 
o—o €, bei 20°C [12], 

o--e, bei 20°C | 


, [13] und vorliegende Arbeit. 
e  «, bei —70°C | 


Dipole im angelegten Feld herrithrt. Da in diesem Fall ganze Molekile oder 
Molekiilteile umorientiert werden miissen, zeigt dieser Anteil eine Dispersion 
schon im Frequenzgebiet zwischen 10% und 1011 Hz. Die dielektrischen Unter- 
suchungen polarer Stoffe haben sich seit den grundlegenden Arbeiten von 
DEBYE zu einem Ausserst wichtigen Hilfsmittel des Physikochemikers ent- 
wickelt, da ein enger Zusammenhang zwischen der Molekiilstruktur und den 
dielektrischen Eigenschaften besteht [1]?). 

Ferroelektrische und antiferroelektrische Kristalle weisen in der Regel sehr 
hohe Werte der Dielektrizitatskonstanten (DK) auf, besonders in der Nahe 
ihrer Umwandlungstemperaturen. In diesen Fallen sind zwat keine permanen- 
ten Dipole beteiligt. Hingegen stehen die durch die thermischen Gitterschwin- 
gungen erzeugten induzierten Dipolmomente untereinander in derart starker 


3) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeivhnis, Seite 399. 
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Wechselwirkung, dass unterhalb der Umwandlungstemperatur spontan eine 
wesentliche Polarisation von paralleler bzw. antiparalleler Orientierung auf- 
tritt. Es ist daher naheliegend, wie an den Stoffen mit polaren Molekiilen auch 
hier das Dispersionsverhalten zu untersuchen, um Auskunft tiber den Polarisa- 
tionsmechanismus zu erhalten. Diese Méglichkeit ist jedoch bis heute noch 
kaum ausgentitzt worden. Im Bereich der Zentimeter- und Millimeterwellen, in 
den die Dispersion fallt, ist es schwierig, DK-Werte von 100 — 10000 genau zu 
messen, und man benétigt dazu fiir jede Frequenzoktave eine verschiedene 
Apparatur. Aufschlussreiche Resultate bedingen Messungen an Einkristallen 
verschiedener Orientierung, die oft nicht in geniigender Grdésse geztichtet wer- 
den kénnen. Ausserdem ist es nétig, das Dispersionsverhalten als Funktion der 
Temperatur zu untersuchen, insbesondere oberhalb und unterhalb des Um- 
wandlungspunktes. 

Einigermassen vollstandige Angaben bis zu Frequenzen von 3 - 101° Hz 
(1cm Wellenlinge) liegen fiir die Erdalkalititanate in polykristalliner, kera- 
mischer Form vor. Diese Materialien haben technisch als Dielektrikum fiir 
raumsparende Kondensatoren grosse Bedeutung erlangt. Die in Figur 1 zusam- 
mengestellten Daten [2] wie auch die Messungen an (Ba—Sr)TiO;-Mischkristal- 
len [3-5] zeigen, dass in der paraelektrischen Phase*) bis tiber 101° Hz hinaus 
keine Dispersion auftritt. Der Abfall der DK auf den Wert nu? ~ 6 erfolgt im 
Gebiet von 1012 Hz, wo an BaTiO,- und SrTiO,-Kristallen verschiedene Ab- — 
sorptionsbanden beobachtet wurden [6]. In der ferroelektrischen Phase von 
BaTiO, oder von Mischkristallen dagegen beginnt die DK bereits oberhalb 
3+ 108 Hz abzusinken und strebt oberhalb 3 - 101° Hz einem Wert e’ & 100 zu. 
Dieser Abfall im Zentimetergebiet wurde von KitTTEt [7] mit der Tragheit der 
Bewegung der Domanenwande vom BaTiO, gedeutet. 

In einer neueren Arbeit haben AKAO und Sasaxr[8] die Dispersion von 
Seignettesalz im Temperaturintervall — 30 + + 35°C untersucht. Sie fanden 
eine Dispersionsfrequenz®) von der Gréssenordnung 3 - 10° Hz, die nur schwach 
temperaturabhangig ist. Insbesondere hat das Auftreten der spontanen Polari- 
sation selbst keinen Einfluss auf die Relaxationszeit. Das Dispersionsverhalten 
entspricht nicht einer einfachen Debye-Dispersion. Es tritt bei allen Tempera- 
turen eine Relaxationszeiten-Verteilung auf, die mit den bestehenden Theorien 
der Ferroelektrizitat des Seignettesalzes nicht gedeutet werden kann. 

Beim Antiferroelektrikum ist die Bildung einer Domanen- oder Zwillings- ) 
struktur elektrostatisch nicht erforderlich, da makroskopisch keine spontane 
Polarisation auftritt. Bei den bis jetzt bekannten Antiferroelektrika, deren 
Symmetrie in der Regel beim Ubergang in die antiferroelektrische Phase er- 
niedrigt wird, bildet sich dennoch eine Zwillingsstruktur, welche ermoglicht, 


*) Paraelektrisch nennt man die Phase oberhalb der Umwandlungstemperatur, wo fiir die 
DK ein Curie-Weiss-Gesetz e’ = C/(T — 0) gilt. 
4) Vp = 1/27, wot die Relaxationszeit bedeutet (vgl. Seite 398). 


Vol. VIII, 1957 Hine Stabchenmethode zur DK-Messung 385 


_ dass der Kristall als ganzes den hochsymmetrisclien Habitus der Hochtempe- 
raturphase beibehalt. Eine Ausnahme bildet hierbei (NH,),H; JO, (Diammo- 
nium-orthoperjodat): Seine trigonale Symmetrie bleibt beim Phasentibergang 
erhalten, und entsprechend wurde eine Zwillingsbildung bisher nicht beobachtet 
[9], im Gegensatz zum Silbersalz Ag,H, JO, [10, 11]. Somit kann beim Ammo- 
niumperjodat keine Relaxation von Wandverschiebungen vorhanden sein. Es 
tiberrascht daher, dass die Messungen von ¢’ (e’ parallel zur kristallographi- 
schen c-Achse), welche ABOAV [12] mit einem Q-Meter (Marconi Type TF 329 G), 
ausfiihrte, bereits im Frequenzgebiet von 107 + 5-107 Hz einen merklichen 
Abfall der DK ergaben (Figur 1). Es schien uns deshalb besonders interessant, 
gerade bei diesem Antiferroelektrikum durch DK-Messungen im Bereich der 
Zentimeterwellen weitere Auskunft iiber das Dispersionsverhalten und insbe- 
sondere tiber dessen Temperaturabhangigkeit zu erhalten. 
Ammoniumperjodat-Kristalle lassen sich aus wasseriger Lésung mit Ammo- 
niakiiberschuss in hexagonalen Platten senkrecht zur c-Achse ziichten [9, 12, 13]. 
Die Herstellung ist jedoch recht heikel wegen der geringen Léslichkeit des 
Salzes. Die gréssten erhaltenen Kristalle erreichten einen Durchmesser von 
etwa 15 mm und eine Dicke von etwa 2 mm. Infolge der beschrankten Kristall- 
dimensionen kam von den zahlreichen in der Literatur beschriebenen Methoden 
zur DK-Messung bei Zentimeterwellen (zusammengestellt in (2, 14, 15]) nur eine 
der Stabchenmethoden [16-19] in Frage und somit nur die Messung von é, 
(vergleiche Seite 389). Dies entspricht der Richtung, in welcher beim Phaseniiber- 
gang die fiir die Antiferroelektrizitat charakteristische Uberstruktur auftritt [20]. 
Gréssenordnung und Verlauf der zu messenden komplexen DK waren zum 
vornherein nicht bekannt, die Resonanzmethode war deshalb nicht zweckmas- 
sig. Wir beschrankten uns auf eine Variante, die einen minimalen apparativen 
Aufwand erforderte. Der orientierende Zweck unserer Messung gestattete uns, 
dafiir eine geringere Messgenauigkeit in Kauf zu nehmen. Die verwendete 
Methode wurde durch Messung von Priifsubstanzen bekannter DK gepriift. Die 
gesammelten Erfahrungen werden im folgenden ausfiihrlich wiedergegeben. 


2. Messmethode. 


2.1 Prinzp (Figur 2) 


Ein Mikrowellengenerator erzeugt im Innern eines Hohlleiters (HL) fort- 
laufende Wellen, welche vom abschliessenden Hohlleiterkurzschluss reflektiert 
werden. Der Attenuator (Dampfungsglied) verhindert, dass ein zu grosser Teil 
der reflektierten Energie zum Generator zuriickgelangt (Resonanzeffekte!). 
Durch Uberlagerung der einfallenden und der reflektierten Welle entsteht zwi- 
schen Kurzschluss und Attenuator eine stehende Welle, deren Verlauf mit 
Hilfe einer verschiebbaren Sonde (Messleitung) ausgemessen wird. Die an der 
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Sonde auftretende Hochfrequenzspannung wird gleichgerichtet und mit einem 


Galvanometer gemessen. 
Im Idealfall (Kurzschluss und HL-Wande total reflektierend) ist zwischen 
Kurzschluss und Attenuator die elektrische Feldstarke am Ort der Wellen- 


Gleichrichter Galvanometer 


Generator | Messleitung —- 
Attenuator Probe — kurzschlub 
Figur 2 
Prinzipschema. 


knoten null (Figur 3). Wird an einer beliebigen Stelle zwischen Kurzschluss und 
Messleitung eine Probe (Kérper aus dem zu messenden dielektrischen Stoff, 
e’ > 1, tgd6 > 0 §)) eingesetzt, so besitzt die hochfrequente Strahlung im Innern 
dieser Probe eine kiirzere Wellenlange als in Luft. Entsprechend verschieben 
sich Maxima und Minima der stehenden Welle zwischen Probe und Attenuator 
um eine Strecke A (Figur 4). Ausserdem nimmt die Feldamplitude an den 
Knoten zu und an den Bauchen ab. Die Grésse von A wird wesentlich durch é’ 
bestimmt, die Anderung des Amplitudenverhaltnisses wesentlich durch tg6. 
Grundsatzlich kann man aus diesen Veranderungen der stehenden Welle auf 


Ort der Probe kKurzschiuB 
yes ort der Probe —-KurzschluB 


al 
is ! 


i ! 
eX Pa Ay e/ 
i] 2 | ! 
eae ee a 
Figur 3 Figur 4 
Feldverlauf ohne Probe. Feldverlauf mit Probe. 


é’ und tg6 zuriickschliessen, wenn die Geometrie der Probe und die Generator- 
frequenz bekannt sind. In guter Naherung lésbar ist das Problem, wenn das 
Dielektrikum in Form eines Stabchens von kreiszylindrischem Querschnitt 


*) Die komplexe DK ¢ = e’ — ie” wird durch ihren Realteil e’ und den Verlustfaktor tg 0d =e" |e’ 
charakterisiert. 
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nach Figur 5 durch den HL hindurchgesteckt wird. Diese Stabchenmethode fiir 
HL-Wellen wurde in 2 Varianten schon von LE Bor und LE MoNTAGNER [16-18] 
verwendet. Die Impedanzmessung geschieht bei Proben mit kleinem e’ und 
tg am empfindlichsten in einer kurzgeschlossenen Leitung, bei hoheren Werten 
dagegen in einer reflexionsfrei abgeschlossenen Leitung. 


Probe | Hohlleiter 
J 


Figur 5 
Hohlleiter mit Stabchen. 


Fir unsere Zwecke musste die Stabchenmethode aus folgendem Grund 
modifiziert werden: In dem zu untersuchenden Temperaturintervall von 
— 70 — + 20°C war zu erwarten, dass die komplexe DK sowohl kleine als 
auch grosse Werte annimmt. Wahrend einer Messung bei tiefen Temperaturen 
lassen sich aber Kurzschluss und reflexionsfreier Abschluss aus technischen 
Griinden nicht leicht gegenseitig auswechseln. Wir verwendeten ausschliesslich 
die kurzgeschlossene Leitung. Die Ausmessung der stehenden Welle geschah 
mit einer geeichten Messleitung, wahrend der Kurzschluss fiir die Dauer der 
Messung einer bestimmten Probe unverandert blieb. Fiir jede Probe wurde 
durch Verschiebung des Kurzschlusses diejenige Stellung gewahlt, in der die 
-Apparatur am empfindlichsten misst. Das inverse Verfahren (geeichter Kurz- 
schlusskolben, feste Sonde [16—-18]) eignet sich in unserem Fall nicht, weil die 
exakte Messung der Kolbenstellung temperaturempfindlich ware. 


2.2 Der Hohileiter als Leitung [21] 


Im verwendeten HL (siehe unten) tritt nur die H,)-Welle (TE,)-Welle [21)) 
auf, bei welcher das E-Feld iiberall parallel zur kleineren Rechteckseite 6 des 
HL-Querschnittes verlauft. Die Leitungs-Wellenlange A, in Richtung der HL- 
Achse betragt dann fiir ein unmagnetisches, dielektrisch verlustarmes Innen- 
medium (u; = 1, é; beliebig, tg 6; <1): 


(1) 


wobei A, = 2 c/w = Vakuum-Wellenlange im freien Raum; @ = Generator- 
frequenz; Ag = 2 a = Grenzwellenlange; a = grossere Rechteckseite des HL- 
Querschnittes. Ag ist definiert als das A) der tiefsten Frequenz, welche theoretisch 
im gegebenen HL noch eine fortpflanzungsfahige Welle erzeugt. 
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2.3 Dielektrischer Zylinder im Hohlleiter 


Ein dielektrischer Kérper im Innern eines HL wirkt als elektrisches « Hin- 
dernis», das heisst, die elektrischen und magnetischen Feldlinien werden gegen- 
iiber dem Verlauf bei ungehinderter Ausbreitung verzerrt. Daraus folgt, dass 
im Bereich zwischen einem Ort «vor» dem Hindernis und einem Ort «nach» 
dem Hindernis der Feldverlauf als Summe vieler verschiedener Wellentypen 
von nicht vernachlassigbarer Amplitude dargestellt werden muss [21] und da- 
mit der einfachen Leitungsgleichung nicht mehr gehorcht. Dagegen sind in 
geniigender Entfernung vom dielektrischen Kérper die Amplituden der stéren- 
den Wellentypen so weitgehend abgeklungen, dass sie vernachlassigt werden 
kénnen. Man kann deshalb das Hindernis und seine nahere Umgebung stets 
durch einen aquivalenten elektrischen Vierpol ersetzen. Ist die geometrische 
Ausdehnung des dielektrischen Kérpers in Richtung der HL-Achse gentigend 
klein gegen den HL-Querschnitt, so reduziert sich der Vierpol auf eine einfache 
Shuntimpedanz. 

Steht ein kreiszylindrisches Stabchen vom Durchmesser d und aus einem 
homogenen Dielektrikum der DK ¢ = e’ — 7 «” in einem rechteckigen H,, — HL 
nach Art der Figur 5, so betragt nach [21] die 4quivalente Shuntimpedanz: 

2 (4) - +433 
Zig e—1\xd 4 e—1 


2g 
17 


4a 
log x 2 
: Lig 


a 
mit 
il 


a 1 
Pe, Gat eV 357 2. 


Z, ist die charakteristische Impedanz [22] des HL. Die Beziehung gilt auf ein 
paar Prozente genau, falls dja < 1/10 und | ¢| < 20. 

Kennt man umgekehrt die spezifische Admittanz Yj = Yp/Y, = Z.|ZD 
= G' +1 B’") der Probe, so erhalt man daraus ihre DK: 


Rat, 2 B I Y 
oD Riga pes ont ee i ea mn 2) 
mit 7 
_ Bs 4a Z | Ve 
A eee pat, Oe a 
PB jeune an ey en) oe | D=2 (se), 
uw G | 2 Ag 
y=A B® G | A= = 


Unter den Voraussetzungen, unter denen diese Beziehungen gelten, werden beim 


*) Y. heisst die charakteristisc “ 5 . wick 
; c stische Admittanz der Leitung, G’ die ifi 
B’ die spezifische Suszeptanz der Probe. 8; spezifische Konduktanz und 
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Einfitbren der dielektrischen Probe in den HL dic E-Feldlinien in der nachsten 
Umgebung nur wenig aus ihrer Richtung parallel zur Stabchenachse abgelenkt. 
Beriicksichtigt man zudem, dass die Feldlinien beim Eintritt ins Innere des 
Dielektrikums von der Oberflachennormale weg gebrochen werden, so kommt 


man zum Schluss, dass das E-Feld im ganzen Innern praktisch parallel zur 
Zylinderachse steht. Im Falle eines anisotropen Dielektrikums misst man 
somit die DK in dieser Richtung. 


2.4 Impedanzmessung im kurzgeschlossenen Hohlleiter 


Die Shuntimpedanz Zp liegt parallel zur Impedanz Z, des leeren Leitungs- 
stiickes zwischen Probe und Kurzschluss. Da sich parallel liegende Admittanzen 
addieren, erhalt man Yp als Differenz der Admittanzmessung Y,, mit einge- 
setzter Probe und der Messung Y, one Probe: Yj = Y/, — Yj. Nach der all- 
gemeinen Leitungstheorie kann man die Admittanz an einer Leitungsstelle 
berechnen aus dem Verlauf der stehenden Welle « vor» dieser Stelle; sowohl fiir 
die Messung mit Probe als auch fiir die Messung ohne Probe gilt 


a al 
Yas 


y ist der komplexe Reflexionskoeffizient und betragt seinerseits 


sevice y . 

Srsiam Hy EXD 2c)’ - 
Hierbei bedeutet 7 = Emin/E wax das Verhaltnis von minimaler und maximaler 
Amplitude der eleKtrischen Feldstarke der stehenden Welle®), « = ky x1; %, ist 
der Abstand des 1. Feldstarkeminimums vom Ort der Probe und ky = 2a/Ay. 
Damit wird 

Poa cf = G Am tS Cm 1 Genie No t& Xo 
B Nm — 1 tEe&m No — 1 tE%, 


Zerlegt man Yp in Real- und Imaginarteil, so erhalt man 


Gis Nm (1 Be fg? Om) Bs No (A + tg* a6) (3a) 

‘sae. tes sia te? o;, 6 ‘ile tg7a5 , 7 

B’ (= Nin) t& %m ar a No) tE % ‘ (3b) 
Nin + t8* om Ne + tg*x, 


3. Apparatur 
Als Generator wurde ein Reflexklystron [24] vom Typ «2K25» im Frequenz- 


f intervall (8,73 — 8,83) - 10° Hz verwendet, entsprechend A, = 3,44 + 3,40 cm 


= 


it Wes 2 i Dini 


8) Im englischen Sprachgebrauch mit «inverse voltage standing wave ratio» bezeichnet. 
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oder Ay = 5,20 — 5,07 cm. Die Reflektorspannung wurde vor jeder Messung 
(Anderung der Stellung des Kurzschlusskolbens) auf maximale Ausgangslei- 
stung eingestellt, um die Leistungsschwankungen auf ein Minimum zu redu- 
zieren. Jede Anderung der Reflektorspannung bewirkt gleichzeitig eine leichte 
Verschiebung der Frequenz. A, wurde deshalb vor jeder Impedanzmessung neu 
bestimmt. Als Hohileiter wurde ein Al-Rechteckprofil mit polierten Innen- 
flachen und den Innenmassen 22,86 x 10,16 mm verwendet. Als Durchfithrung 
des HL durch die Kryostatwand diente ein 10cm langes Profilstiick aus 
diinnem Neusilberblech. Zur Vermeidung von Warmekonvektion und Kondens- 
wasser enthielt diese Durchfiihrung ein Glimmerfenster von rund 0,01 mm 
Dicke. Im HL-Stiick zwischen Messleitung und Kurzschluss wurden Kontakt- 
kupplungen verwendet. Die tibrigen Kupplungen waren kapazitiv [25]. 

Der Attenuator [14] wurde auf einen Dampfungswert von 12 db eingestellt. 
Die Messleitung («slotted line») war vom Institut fiir Hochfrequenztechnik an 
der ETH gebaut worden und besass eine verschiebbare Sonde von 2,5 mm Ein- 
dringtiefe und 0,5 mm Durchmesser. Wir versahen sie mit Feinantrieb und 
Messuhr, wodurch eine Messgenauigkeit von 0,01 mm médglich wurde. Der 
Gleichrichter (Si-Diode «1N23») war mit der Messleitung zusammengebaut. Er 
besass im verwendeten Bereich in guter Naherung [26] eine quadratische J_/U_- 


Charakteristik, der Richtstrom war also proportional der Intensitat des E-Fel- 
des. Zur Messung des Richtstroms diente ein Lichtzeiger-Galvanometer von 
0,9 wA Vollausschlag und einer Ablesegenauigkeit von 10-® A im empfindlich- 
sten Bereich. 

Der Kurzschlusskolben (Figur 6) war verschiebbar und wurde auf der Ober- 
und Unterseite durch Teflonquader geftihrt. Die vom Generator her von links 


WNe ee es 


PARE 
Messing Teflon(OK=E;) Eisenpulver 
(Absorber) 
. Figur 6 
Kurzschlusskolben. 


einfallende Welle dringt zu einem kleinen Teil in den obern und untern Luft- 
spalt ein, erfahrt an den verschiedenen Luft-Teflon-Grenzflachen Mehrfachre- 
flexion und erreicht teilweise die Absorberschicht (hinter der keine Reflexion — 
mehr auftritt). Die speziell gewahlten Langendimensionen des Kolbens bewir- 
ken, dass sich eindringende und reflektierte Strahlung weitgehend ausléschen. 
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Der Absorber verhindert, dass sich der Grad dieser Ausléschung und damit der 
Reflexionskoeffizient des gesamten Kolbens mit der Kolbenstellung andert. 
Der ganze HL-Abschluss mit Kurzschlusskolben befand sich im Kryostaten. 

Der Kryostat bestand aus einem durch eine 4 cm dicke Schaumplastik- 
schicht isolierten Trog, der ein Alkoholbad enthielt. Die Abkihlung des Bades 
erfolgte durch Trockeneis, die Erwarmung mit Hilfe einer Heizspirale. Die 
Erwarmung des Bades infolge von Warmeleitungsverlusten betrug bei — 70°C 
rund 1/7°C/min, bei héheren Temperaturen entsprechend weniger. Die Proben- 
zylinder waren je auf einen Halter aus Messing gekittet und wurden durch eine 
Bohrung in der Oberseite des HL in diesen eingefiihrt. Die Temperatur wurde 
an der Aussenfliche des HL, unmittelbar unterhalb der Probe, mit einem ge- 
eichten Cu-Konstantan-Thermoelement gemessen. Ein Versuch hatte ergeben, 
dass die Temperatur hier hinreichend genau mit der Lufttemperatur im Innern 
des HL iibereinstimmt. Vor und nach jeder Impedanzmessung (Dauer 1 + 4 
min) wurde je eine Temperaturmessung vorgenommen, wobei der Mittelwert 
verwendet wurde. 


4. Experimentelles 
4.1 Herstellung der Proben 


Von den zu untersuchenden Substanzen wurden zylindrische Stabchen von 
mindestens 13 mm Lange und mit einem Durchmesser zwischen 1 und 5 mm 
hergestellt. Die Zylinder aus Teflon, Polyathylen und Plexiglas wurden auf der 
Drehbank bearbeitet. Das destillierte Wasser wurde in einem Glasréhrchen von 
etwa 0,05 mm Wanddicke gemessen. Nach jeder Messung des gefiillten Rohr- 
chens wurde noch der leere Behalter gemessen und dessen Admittanzwert von 
demjenigen des vollen Réhrchens subtrahiert. 

Die Kristalle wurden zu Stabchen von quadratischem Querschnitt gesagt. 
Zum Rundschleifen wurde das Stabchen auf einen Halter (Figur 7) gekittet, 
welcher um einen Winkel + 45° um die Stabchen- 
achse hin- und hergedreht wurde. Eine horizontal 
gelagerte, glatte Cu-Platte, die mit einer Paste 


ZK 


von Karborundpulver und Petrol bestrichen war, —__ iste | : 2 lee 
iibte von unten her einen leichten Druck auf die 

zu schleifende Kristallkante aus. Analog wurden Figur 7 
nacheinander auch die iibrigen Kanten rund- __ Kristallhalter fiir Rundschleifen. 


geschliffen. Die Ungenauigkeit des Durchmessers : 

der so bearbeiteten Kristalle betrug héchstens + 0,1 mm. Die mechanische 
Beanspruchung der Stabchen beim Rundschleifen war gross; rund die Halfte 
der Kristalle zerbrachen. Bessere Resultate diirften sich bei der Verwendung 
~ eines Pulverstrahl-Apparates erzielen lassen, welcher sich bei der Bearbeitung 
von spréden organischen Finkristallen bewahrte [27]. 
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4.2 Bestimmung von tga, Einstellung des Kurzschlusskolbens 


Die Grosse tga, [Gleichungen (3)] kann man auf die feste Distanz / zwischen 
Probe und Kurzschluss zuriickfiihren. Aus Figur 3 sieht man: x, = ”Aq/2—1 


(n = ganze Zahl), also tea, = —tg (Ry) . (4a) 


Nach Einsetzen der Probe hat man nach Figur 4: *,,, = %1, — 4, also 
tg Om - —tg[ky (J na Ay) fi (4b) 


Der Ort des Minimums der stehenden Welle wurde stets indirekt bestimmt 
als arithmetisches Mittel der zwei x-Werte, aus denen auch 4x (siehe Seite 394) 
bestimmt wurde. Die Annahme, dass das Minimum in der Mitte zwischen den 
Flankenpunkten gleicher Intensitat liegt, ist infolge der Sondenreflexion [14] 
nicht exakt richtig. Der dadurch verursachte Fehler erwies sich als vernach- 
lassigbar. Bei den Proben mit e’ ~ 2 wurde der Ort des Minimums versuchs- 
weise auch direkt bestimmt: Bei d = 1mm wurden auf diese Weise die Resul- 
tate bedeutend schlechter, wahrend sich bei dickeren Proben gleichgute Resul- 
tate ergaben. 

Unsere Erfahrungen zeigten, dass grundsatzlich um so genauere Resultate 
erzielt werden, je groésser tga, und tga,, sind. Bei der Ermittlung der opti- 
malen Stellung des Kurzschlusskolbens hat man daher 3 Falle zu unterscheiden : 


a) Bei kleinen Impedanzen [(«%, — &m) > 2/2] wird man den Kolben so ein- 
stellen, dass a, und «,, symmetrisch zu N a liegen (N = ganze Zahl). Daraus 
folgt als Bedingung: 2 
[oe A 

2 a 

b) Bei matileren Impedanzen (7/12 S («) — &m) << 2/2] sollen «, und «,, sym- 
metrisch zu 2/2 + N aliegen, weil in diesem Fall tga, und tg«,, grésser werden. 
Daraus folgt , 


i= (w+ 5) H 4. 


c) Bei grossen Impedanzen [ (a) — &,) </12] wiirden bei obiger Kolben- 


stellung die tga oo gehen, so dass die Methode ungeeignet wird. Man wahlt- 


in diesem Falle am besten: I 
le (N+ 5] Ft. 
% besitzt dann den giinstigen Wert 


(0 *a = 
Chere os 


und «,, liegt um weniges tiefer, das heisst, tga, wird grosser (ohne aber > co : 


a ae 
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zu werden). Messungen mit 


Im (ve 3) 


bo YO - 5 +Na 
und tga,, gegen kleinere Werte hin verschoben wird, ergaben schlechtere 
Resultate. 

Der voraussichtliche Wert von A muss rechnerisch abgeschiitzt oder durch 
eine Messung in erster Naherung bestimmt werden. N ist so klein wie experi- 
mentell moglich zu halten, damit der Leitungsverlust des HL die Messgenauig- 
keit nicht unnétig vermindert. 

Bei | e | > 20 stellt die Formel fiir Zp auf Seite 388 eine schlechte Approxi- 
mation dar, und bei noch héheren DK-Werten ist die Darstellung des aquiva- 
lenten Vierpols durch eine einfache Shuntimpedanz nicht mehr erlaubt (syste- 
matische Fehler!). Ausserdem kommt dann die Oszillatorfrequenz 101° Hz in 
die Nahe der Resonanzfrequenz des Vierpols. Die Vierpolelemente sind dann 
nicht mehr stark von der DK abhangig (Streuung!). Im Gebiet von e’ ~ 50, 
tgdé ~ 0 liess sich in der Kolbenstellung / = (N + 0,14) Ay/2 reproduzierbar 
messen. 

4.3 Lettungsverlust 


Infolge von Verlusten an den HL-Wanden nimmt die Intensitat der fort- 
laufenden Welle vom Ort der Probe bis zum Ort der Sonde auf einen Bruchteil 


ab: 1 


wobei /; die Lange eines HL-Teilstiickes, #; dessen Leitungsverlust in Dezibel 
pro Langeneinheit. Hierbei ist iiber alle HL-Teilstiicke zwischen Probe und Sonde 
- zusummieren. Der Wert von #, hangt in einfacher Weise von der elektrischen Leit- 
'fahigkeit o, des HL-Metalles ab [28]: #;0¢ o;-1”. Die Reduktion der Leistung 
auf den Bruchteil # hat zur Folge, dass der mit der Messleitung gemessene Wert 
Ewmin|Emax = 7 grésser ist als der gesuchte Wert 74 am Ort der Probe. Die all- 
gemeine Leitungstheorie ergibt 
ye bE ees 1) 
Dl tig a ee ae 


oder aufgelést nach 7 


eas LP) (6) 


4.4 Bestimmung von 1’ 


Messungen von Feldamplituden verschiedener Gréssenordnung werden 
wegen der iiber gréssere Bereiche nicht mehr quadratischen Gleichrichtercha- 
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rakteristik verfalscht. Bei kleinen 4’ misst man deshalb nicht Eyyq, Undrleses 
sondern den Abstand Ax derjenigen zwei Punkte beidseits des Minimums, bei 
welchen die Feldstairke V2 E,,;, betragt, also die Intensitat 2 I j,i, . Nach [2] 


gilt dann 


tee 


Der Korrekturterm 6 beriicksichtigt die Ungenauigkeit der Naherung. Er ist in 
[2] tabelliert und nur fiir sehr kleine 7’ vernachlassigbar. Ax wird fiir sehr kleine 
ny’ oft unmessbar kurz; man kann dann statt dessen die Kurvenweite A, x 
zwischen den zwei Punkten der Intensitat 1 I,,;, messen, wobei nach [14]: 
Ax = A, x|/Vn — 1. Diese Naherungsformel ist mit = 10 bequem. 

Zeitliche Schwankungen der Generatorfrequenz bewirken Schwankungen 
der ortlichen Feldstarke. Zur Bestimmung der Sondenstellungen mit J = 1 J pip 
wurden diese Intensitatsschwankungen wahrend eines langeren Zeitintervalls 
beobachtet und beriicksichtigt (rund 1 min). 


5. Resultate 


5.1 Priifsubstanzen und Phosphate 


Tabelle 1 € 
vena 
Dislekiriicara d 7 ted - 102 ergleichswerte [2] 
mm e tgd - 102 
Teflonsss ale i Ol ale Oana <7 2,08 0,04. | 
Polyathylen . | 1,05 2:0 4 0715 as 2025 0,04 
Plexiglas. . . | 0,99 FMA <8 2,59 0,67 
1,97 2 40410, 1 <2,5 
3,04 2,6) a 001 <i75 
4,00 2i5y e004. <2,5 
5,04 2,62 + 0,06 <3 
NEHPOs | ¢: |) 182 1142 441,5 | 13,7 0,5 
eae 52 39 + 3 1+0,5 | 55,91) | <0,14) 
KH,PO,|¢ ..| 1,42 iy JED Pie St 20:22)" | 20,052) 
a . | 1,95 29 441 24+0,3 | 44,38) | <0,053) 
Wasser . .. | 0,800 65 + 1,5 38 + 0,6 544) 644) 
1) Bei 10¢~ 3-108 Hz. 3) Bei 103 — 108 Hz. 


2 : . 
) Bei 106 =~ 108 Hz, 4) Interpoliert fiir 18°C und 8,8 - 109 Hz. 
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Die Temperatur betrug je nach Messung 18 -~ 20°C, die Frequenz (8,73 
— 8,83) - 10° Hz. Die Fehlerschranken geben den moglichen zufalligen Fehler 
wieder. Die Vergleichswerte [2] verstehen sich fiir 22 ~ 27°C und, wo nichts 
anderes vermerkt, fiir 1019 Hz. 


5.2 Ammoniumperjodat 


Die Temperaturkurve der DK senkrecht zur c-Achse ist in Figur 8 wieder- 
gegeben. Die Messreihe mit fallender Temperatur ist durch volle Kreise dar- 
gestellt, diejenige mit steigender Temperatur durch leere Kreise. Die Generator- 
frequenz betrug 8,74-10° Hz, entsprechend Ay, = 5,18cm. Die Probe hatte 
einen Durchmesser von 1,57 mm. Der maximal mégliche zufallige Fehler 
wurde graphisch wiedergegeben. Die gestrichelte Kurve gibt den Verlauf der 
statischen DK nach BaErtscut [13]. 
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Figur 8 
Temperaturabhangigkeit von e, und tgdg von (NH,4).H3JOg. 
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6. Diskussion der Messmethode 
6.1 Zufillige Fehler 


Der Einfluss einer grossen Zahl méglicher Fehlerquellen wurde abgeschatzt. 
Wichtig sind nur die Schwankungen der Frequenz und die Abweichung der 
Gleichrichtercharakteristik vom quadratischen Verlauf. Bei Einkristallen, 
welche wegen ihrer Sprédigkeit nicht auf der Drehbank bearbeitet werden 
kénnen, spielt auch die Abweichung von der geometrischen Form eines Zylinders 
eine Rolle. 

Die in Figur 8 eingezeichneten Fehlerschranken sind offensichtlich zu gross. 
Das war zu erwarten, da der Einfluss der Frequenzschwankungen betrachtlich 
reduziert werden konnte, indem alle Intensitaétswerte wahrend eines langeren 
Zeitintervalles beobachtet wurden (Seite 394). Die den Resultaten in Tabelle 1 
beigefiigten Fehlerschranken enthalten nur einen Anteil von 1/3 der maximal 
méglichen Frequenzschwankungsfehler. 

Beobachtete Streuung: Die Ungenauigkeit von e’ und tg 6 ist eine Funktion 
der DK und des Stabchendurchmessers. Wir beobachteten fiir tg 6 < 10-? und 
d= 1,5 mm im Mittel folgende Streuungen der Resultate: 


ms Streuung von &’ Streuung von tg dé 
ag (relativ) (absolut) 
25 + 5% + 0,005 
15 + 5% + 0,01 | 
50 + 10% + 0,01 | 


Bei e’ < 15 war die beobachtete Streuung kleiner als die geschatzte. Bei 
é’ » 50 war sie zu gross, was vermutlich auf die verringerte Empfindlichkeit 
der Methode in diesem DK-Bereich (Seite 393) zuriickzufiihren ist. 


Ry AS) ystematische Fehler 


Die Darstellung des elektrischen Hindernisses durch eine Shuntimpedanz 
und die quantitative Berechnung von Zp (Seite 388) beruhen auf einer Nahe- - 
rung, die nur fiir diinne Proben und niedrige DK gilt. Die Verfalschung des 
Messresultates in Abhangigkeit von der Probendicke wurde fiir Plexiglas von 
1 — 5mm Durchmesser untersucht (Tabelle 1). Ein systematischer Gang von 
e’ und tg6 konnte fiir diesen DK-Bereich jedoch nicht festgestellt werden. Die 
Abhangigkeit von der DK kann durch Vergleich mit genaueren Literaturwerten 
beurteilt werden. Im Bereich e’ ~ 15 liegt nur ein Wert fiir e/ von NH,H,PO, 
vor, der bis 101° Hz noch keine Frequenzabhangigkeit aufweist. Fiir ef von 
KH,PO,, das nur bis 108 Hz untersucht ist, darf man dasselbe annehmen. Falls 
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man alle in Tabelle 1 aufgefiihrten Literaturwerte als richtig voraussetzt, er- 
geben sich fiir unsere Messmethode bei tgd < 10-? und d= 1,5 mm folgende 
Verfalschungen der Resultate: 


Systematische Fehler 
von tg6 (absolut) 


Systematische Fehler 
von é’ (relativ) 


+ 0,007 


2,5 ging 5% 
+ 0,03 


—15% 


Bei héheren DK-Werten fehlen zuverlassige Vergleichswerte. Unsere Resultate 
fir e, von NH,H,PO, und KH,PO, liegen um 30 + 35% tiefer als diejenigen 
in [2], welche aber bis heute erst bis 3 - 108 Hz bzw. 108 Hz gemessen sind. Es 
ist denkbar, dass die Differenz mindestens teilweise reell ist, da e/ der Phosphate 
moglicherweise eine niedrigere Dispersionsfrequenz besitzt als ¢’, und dass sich 
die Dispersion bereits bei 10!°Hz in einem DK-Abfall bemerkbar macht. 
Unsere Messung an Wasser hat nur orientierenden Charakter, da man sich bei 
101° Hz mitten im Dispersionsgebiet befindet und da zudem das dielektrische 
Verhalten wesentlich von der Reinheit abhangt. 


6.3 Zusammenfassung 


Das kreiszylindrische Stabchen parallel zum E-Feld ist, abgesehen von der 
Platte senkrecht zur HL-Achse, die einzige dielektrische Diskontinuitat, deren 
aquivalente Impedanz in einer guten Approximation berechnet werden konnte. 
Die hier beschriebene Variante der Stabchenmethode gestattet, mit einer ver- 
haltnismassig einfachen Apparatur eine (isotrope oder anisotrope) DK von 
Stoffen zu messen, die nur in kleinen, langlichen Stiicken zur Verfiigung stehen. 
Die Variante eignet sich speziell fiir orientierende Messungen, bei denen die | 
Gréssenordnung der DK nicht zum voraus bekannt ist oder bei der Messung als 
Funktion der Temperatur wechselt. Direkte Prazisionsmessungen sind an Stof- 
fen mit kleinem e’ und tg6 méglich, wenn man die Messmethode weiter verfei- 
nert: 1. Stabilisierung der Frequenz (Resonator); 2. Entkopplung des Gene- 
rators mit Hilfe eines Gyrators; 3. exakte Eichung des Gleichrichters; 4. Nie- 
derfrequenzmodulation des Generators und Verstaérkung des Diodenstroms, 
womit sich wiederum der Generator besser entkoppeln lasst. Bei Stoffen mit 
hoherer DK ist es notwendig, die Absolutgenauigkeit der Methode durch Eich- 
messungen zu erhohen, wofiir jedoch bei Frequenzen von 2 101° Hz nur wenige 
geeignete Literaturwerte vorliegen. 

In einfachen Fallen kénnen die Nomogramme von [16—18, 29] verwendet 
werden. Bei unserer Messanordnung erwies es sich als notwendig, alle Bezie- 
: hungen und Korrekturen der Gleichungen (1) bis (7) zu verwenden, und eine 
Aufstellung besonderer Nomogramme lohnte sich nicht. 
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7. Diskussion der Messungen an Ammoniumperjodat 


Wie aus Figur 8 ersichtlich, zeigt ¢, bei 3 cm Wellenlange einen analogen 
Verlauf mit der Temperatur wie die statische DK. Man kann mit grosser 
Wahrscheinlichkeit annehmen, dass die Dispersion im Zentimetergebiet als 
Relaxationsprozess und nicht als Resonanzabsorption [30] zu deuten ist. 

Die einfachste Art von Relaxationsdispersion ist die einer Debye-Dispersion, 
gekennzeichnet durch eine einzige Relaxationszeit t (monodispers). Es gelten 
dann fiir Real- und Imaginarteil der DK die. Debye-Gleichungen 

é,, — & 64 — € 
6 (@) = &,, 7 eae e"(o) = een OT 
In praxi findet man jedoch oft, dass der Abfall von e’ im Dispersionsgebiet im 
Vergleich zum Debye-Fall sich iiber ein grésseres Frequenzintervall erstreckt 
und gleichzeitig die Kurve e”(w) verflacht ist. Dieses Verhalten erklart man 
sich dadurch, dass nicht mehr ein einheitliches t vorliegt, sondern eine gewisse 
Relaxationszeiten-Verteilung (polydisperses System). 

Der hochfrequente Grenzwert ¢,, der DK enthalt im Falle des Ammoniumper- 
jodats den Elektronen- und den Infrarotanteil der Polarisation. Der Elektronen- 
anteil wurde von uns bestimmt durch eine Messung des Brechungsindex bei 
sichtbarem Licht. Nach der Methode des Duc DE CHAULNEs wurde fiir die gelhe 
Hg-Linie (578 my) n, = 1,71 + 0,02 gefunden. Da 1, etwa von gleicher Grésse 
sein diirfte, ist e,,,, ~ 3. 


In der Kristallstruktur des Ammoniumperjodats bilden die JO,-Oktaeder — 


Komplexe mit vorwiegend homdopolarer Bindung zwischen J- und O-Atomen. 
Die JO,-Komplexe sind als ganzes fiinffach negativ geladen. Die Hydrogen- 
bindungen, welche benachbarte Oktaeder verbinden, und die Bindungen zu den 
(NH,)*-Ionen sind vorwiegend ionogenen Charakters. Es ist daher mit einem 
verhaltnismassig grossen Infrarotanteil (« lonenpolarisation » [30]) zu rechnen. 

Um zwischen mono- und polydispersem Verhalten entscheiden zu kénnen, 


geniigen die DK-Messung bei einer einzigen Frequenz im Dispersionsgebiet und — 


die Kenntnis von ¢;,,; nicht. Mit Bestimmtheit lassen sich jedoch iiber die DK 


parallel zu a in der Hochtemperaturphase zwei qualitative Aussagen machen, 
auch wenn man beriicksichtigt, dass e/ méglicherweise systematisch zu niedrig - 


gemessen wurde (Seite 397): 


1. Die Dispersionsfrequenz yp = 1/271 ist kleiner als die verwendete — 


Messfrequenz (8,74 - 109 Hz). 

2. Mit sinkender Temperatur nimmt ¢,,, relativ starker zu als ¢,.,. Dies 
bedeutet, dass mit fallender Temperatur yp ansteigt, somit t abnimmt. Darin 
steht Ammoniumperjodat im Gegensatz zu den ferroelektrischen Kristallen 


BaTiO, und Seignettesalz [2, 3-5, 8], bei welchem t am Umwandlungspunkt 
ein Maximum besitzt. 
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Uber die antiferroelektrische Phase lassén sich analoge Aussagen nicht 
machen, obwohl unsere Messungen hier genauer sind. Es sind zusitzliche Re- 
sultate bei andern Frequenzen, insbesondere im Gebiet von 1 cm Wellenlange 
erforderlich. 

Quantitative Ergebnisse erhalt man nur, wenn man rein monodisperses 
Verhalten voraussetzt. Man kann dann t und «4 mit Hilfe der Debye-Glei- 
chungen berechnen. Man erhalt t von der Gréssenordnung 10-!%s, und t(T) 
durchlauft im Umwandlungsgebiete ein Minimum. e’,, wird etwa 46 in der Hoch- 
temperaturphase und etwa 12 in der antiferroelektrischen Phase, die geringe 
Temperaturabhangigkeit ist gleichsinnig wie der Verlauf von ¢4,,,. Man kann 
sich tiberlegen, dass die fiir den monodispersen Fall berechneten ¢/,,-Werte eine 
obere Schranke der méglichen (polydispersen) Werte darstellen. 

Die Verschiebung der Umwandlung zu héheren Temperaturen beim Ersatz 
der H-Atome durch D (um 15° C bei (NH,),H; JO, [31] bzw. 40°C bei Ag,H, JO, 
[10]) zeigt, dass es sich um « order-disorder »-Umwandlungen beziiglich der An- 
ordnung der H-Atome auf den Hydrogenbindungen handelt [9]. Analog wie bei 
KH,PO, sind die Anderungen der H-Anordnung als primare Ursache fiir den 
Konfigurationsanteil der Polarisation aufzufassen. Zur Beurteilung der Polari- 
sationsbeitrage der iibrigen Ionen ist eine genaue Strukturbestimmung im 
Gange. Das von uns beobachtete anomale Verhalten der Relaxationszeit in der 
Hochtemperaturphase steht mit der antiferroelektrischen Wechselwirkung der 
Dipole im Zusammenhang und gibt einen weitern Hinweis auf den Polarisa- 
tionsmechanismus. . 


Wir danken Herrn Prof. Dr. P. SCHERRER fiir sein forderndes Interesse an 
dieser Arbeit und Herrn Prof. Dr. F. TANx fiir die leihweise Uberlassung einer 
Messleitung. Die Herren Dr. W. Nev und K. A. MULLER haben mit wertvollen 
Ratschlagen zum Gelingen der Messungen beigetragen. 
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Summary 


The complex permittivity can be measured in the region of centimeter-waves 
by inserting a cylindrical rod of the specimen across a short-circuited rectangular 
waveguide in the maximum of the E-field. The resulting change in line-impedance 
is measured with a slotted line. This rod method [16-18] was modified and tested 
with materials of known dielectric properties (Teflon, Polythene, Plexiglas and 
water). The conditions for best performance with specimens of high dielectric 
constants were established. Measurements were made at room temperature with 
single-crystals of NH,H,PO, and KH,PO, parallel to the a and c axis. : 
_ The permittivity e, of diammonium-orthoperiodate (NH,),H,; JO, was studied 
in the temperature range of + 20° to — 70°C. In the unpolarized and in the anti- 
ferroelectric phase the real part of the dielectric constant at 8740 Mc/s is less than 
half the static value. The results indicate that in ammonium periodate the disper- 
sion frequency of e4 is in the range of 104 Mc/s and that the relaxation time 
increases with rising temperature above the transition point. This behaviour 
represents the contrary of what is found with ferroelectrics. 


(Eingegangen: 4. April 1957.) 
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Discrete Approximations to the Laplace Transformation’) 


By WoLFrcGanc Wasow, Madison, Wis. 2) 


1. Introduction 


The form and the substance of finite difference methods depend in an 
essential manner on the continuous analogue that motivates them. For electrical 
engineers the Laplace transformation is the most familiar tool for the solution 
of differential and difference-differential equations. It is therefore a natural 
step for them to introduce the sum 


h OH f(a h) ews" 


as a discrete analogue to the Laplace transformation 
Qi = [#0 en df (1.1) 
0 


of the function /(). The computational techniques based on the Laplace 
transformation suggest then approximation methods employing its discrete 
counterpart. It is more common to introduce a new variable 


z= es" (1.2) 


and to drop the factor h. The resulting operation, defined by 
3f= Di nh) a (1.3) 
n=0 


is frequently called the z-transformation of f. We shall affix the subscript s to 3, 
whenever 2 is replaced by s, so that 


h3,f=h 3? (nh) (1.4) 


is the Riemann sum corresponding to  f. | 
The function 3 / is the generating function for the sequence /(m h), 7=0,1,... 
which results from the discretization of f(t). Mathematically there may not be 


1) This paper represents research supported in part by the Office of Naval Research and in 


part by the Mathematics Research Center, US Army, Contract No. DA-11-022-ORD-2059. 
2) Mathematics Research Center, US Army. 
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much difference between dealing with 3 f or dealing directly with the sequence 
[#(0), f(2), ...]. But the two viewpoints suggest the invention of different 
computational schemes. For instance, M. CuéNop (1]*) and R. BoxER and 
S. THALER [2] have each proposed algorithms for the approximate solution 
of differential equations by such finite difference schemes. The former author 
operates directly with the sequence /(7 h) while the authors of [2] are moti- 
vated by the theory of Laplace transformation. The resulting finite difference 
methods have the same accuracy, as far as the order of magnitude of the error 
is concerned. However, the numerical answers as well as the computational 
steps of the respective methods differ in a non-trivial way. 

The papers quoted do not attempt a general mathematical error analysis. 
The purpose of the present investigation is to provide such an analysis for 
some of the techniques based on the z-transformation. 


2. Elementary Properties of the z-Transformation 


The discussion will be limited throughout to functions /f(¢) that are entire 
and of exponential type. This means that f(#) is a holomorphic (i. e., regular 
analytic) function of ¢ in the whole finite complex ¢#plane and that an in- 
equality of the form 

| < Cer!" (2.1) 


where C and x are constants, holds everywhere. While it is possible to extend - 

many of the arguments below to wider classes of functions these generalizations 

do not greatly increase the usefulness of the results to the usual applications. 
With these assumptions on f(t) the series in (1.3) is dominated by the series 


co 


“nh ~—n 
Ce es « 


n=0 


and it follows by the ratio test for infinite series that (1.1) converges for 
DM ibe Ae (202M 


| The z-transform is obviously a linear operation, i.e., f(t) = c, fit) + ce f(t) 
implies 3 f = cy 3h + C23 fy. This property is frequently exploited by applying 
the z-transformation termwise to a series expansion of the function in question. . 
Thus, if 


1) = D7 ee eat ee 


3) Numbers in brackets refer to References, page 417. 


Vol. VIII, 1957 Discrete Approximations to the Laplace Transformation 403 


with ’ bene 
r aarar Hi |, “: (2.4) 
one would write 
3t= 2) 4.3 Gar (2.5) 
r=] 


Since /(#) is an entire function the series in (2.3) converges for all t. Never- 
theless, the relation (2.5) requires a proof, since a sum of infinitely many 
terms is involved. In order to prove (2.5) observe first that 3 f, as defined in 
(1.3), is, for |z| 22+ 6, |h| Shy (6 > 0, hy > 0, arbitrary), a uniformly 
convergent series of holomorphic functions of / and is therefore itself a holo- 
morphic function of 4. Hence, it possesses a convergent series expansion in 
powers of h. The coefficient of #”-1 in this series is 


1 q’-1 co 1 oO = qt-1 ra 
(y —1)! ae = 2 Tinka ag eat ees locas ae ® 


provided the series 


1 . r—1 4(r—1) —n 
reap ye fo (nh) z 
converges uniformly. This it does, since the differentiation of a uniformly 
convergent series of analytic functions may be performed termwise. By virtue 
of (2.4) the second member of (2.6) equals 


hence, 
1 + r—1 y—n pr—1 __ 1 Ass . r—1 ,—n 
3f= : em hee be aR = by Pe poi z 
ies) t7-1 
a 


which proves (2.5). 

The expressions for 3 ¢’-1/(r — 1)! in closed form occur frequently in the 
applications. They will therefore be derived here. Let R,(z) be defined as the 
function with the power series expansion 


R,(2) ork mP a? . (2.2) 
then 
prt Ar} 


2 7 R18) - (2.8) 
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The series in (2. 7) converges for | z | > 1. The functions R,(z) can be written 
as finite sums with the help of the so-called ‘Stirling numbers of the second 
kind’ Sf. (Several tables of these positive integers are available. See [3], page 
168, footnote, and page 170.) One has, in fact, 


Ry(z) = —, (2.9) 
p 
R,,(2) ae) kh! St ae ee . (2.10) 


The identity (2.9) is obvious. The proof of (2.10) is based on the identity 


WP = Sh tra: S30 (ni — A) Aer +S tlt Delp +) (P<auli epee 
(2.11) 


(see [3], page 168). It follows from (2.11) that, for | z| > 1 und p > 0, 


R,(2) = by nb 27” 


n=0 


= {Spm + Spm (n— Vet Sh (1) vo (np he 


LW n=0 
= + Sta? 3? n(n —1) (n—p +127"? 
n=0 
eck ka i = gk 
a See Fa wae ee es Sr! Ter 
b ; 5 
ee Si (z — 1)k+1 


An alternate form of (2.8), obtained by simple algebraic manipulation, is 


pew 
R,(2) = sneer Da eae (2,127) 
with sad 


pu 
Te es (—1)P RR Chases 


If a function /(¢) is given by its power series (2.3) it can, of course, be 
calculated numerically for any value of ¢ by means of a suitably long partial ania 
of the series. An alternate method, which yields the values of {(n h) for all n 
is based on the z-transformation and is recommended, e. g., in [2]. It can be 
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expressed in the following rule, which is an immediate consequence of the 
preceding analysis. 


Rule 1. [} 


replace every 


and expand R,_,(z), as given by (2.9) or (2.10), in a series of powers of z-}. 
If the resulting double series is re-arranged according to powers of 2-1 the coeffi- 
cent of 2-” is equal to f(n h), because of (1.3), (2.5) and (2.8). 


The only error inherent in this numerical procedure is the one caused by 
the truncation of the infinite series involved. 


3. The z-Transformation and the Laplace Transformation 


In view of formulas (1.1) and (1.4) the expression (1.4) can be expected to 
approximate the Laplace transform £7. The difference between the two 
transforms can be expressed by means of the Euler-Maclaurin formula, a per- 
tinent form of which is given in [4], page 50. With the notation of the present 
paper that formula becomes 


Te hy qv-1 Se 
2f—43,f= DB [Gar Oe] + Re (3.1) 
y= ‘ al 
where the remainder R,, is 
mt. ae qm a 
a St af B, (1) a Ca (3.2) 
0 


Here m is an arbitrary positive integer, B, is the y-th Bernoulli number, and 
B,, (2) is the periodic function with period unity that coincides in the interval 
0 <z< 1 with the m-th Bernoulli polynomial. Since the definition of BER- 
NOULLI’s polynomials and numbers that are found in the literature are not 
always equivalent, it must be emphasized that the one meant here is that 
given in [4], page 19. It is easy to verify that the derivation of (3.1) and (3.2) 
given in [4], pages 49-50, is valid under the assumption (2.1), provided s > x. 

The remainder term R,,, is of the order of magnitude O(h™*+1). More precisely, 
there exists for every interval x+ O0Ss5S55, <605 012 0,.a constantyG.. 
independent of s and /, such that 


IR, | C,H? (wt 0SsSs, 0<bSh). (3.3) 
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In order to prove (3.3) set t = Wt in (3.2). It then takes on the form 


pier = T 7 

dk "a ml Elo v6 (es s) é at , 
0 

where y(t, s) is a polynomial of degree m in s whose coefficients are linear in 

the m first derivatives of f(t). The latter satisfy inequalities of the form (234), 


possibly with larger constants C. Hence, 


Pore ess nM YP st [ e- 9 de (M = const) , 
v=0 


0 


and the last integral converges for s > x. Thus, | R,, | SCA nh”, wheresee, 
is a constant, for x + 6 <s <s, < co. However, since (3.1) and (3.2) are valid 


for all m 
? Bs Amrit qm JG 
i, (m alte 1)! Be, E [7 (2) é ee - Rgces , 


which proves the stronger inequality (3. 3). 

A first conclusion that can be drawn from (3.1), (3.2) and (3.3) is that 
2f—h 3, f is of the order of magnitude O(h). More precisely, there exists a 
positive number A such that 


|27-h3,f| SAM (+ d6S5s85, OSX ASH). (3.4) 
Formula (3.1) suggests a refinement of the z-transformation. If 3*f is 
defined by { 
3*f=31-> 10). ae 
then (3.1), (3.2), (3.3) and the fact that B, = —1/2 show that 
Aer ret last ee (3.6) 


This transformation amounts to approximating the Laplace transformation by 
means of the trapezoidal rule. 

Boxer and THALER [2] propose a different modification, which can also be 
motivated by means of formulas (3.1) and (3.2). They define the transformation — 
3** f for f = 7#-4/(r — 1)! by 


(oak Fil 


rt Br 
om (y — i)! me (vy — 1)! +h ee (3.7) 


The equations (3.1), (3.2) and (3.3) show then that 


js | Aone 
Ee TAOS Gay] SA et OSs Sy ee 
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where A, is a constant. Hence 3** leads to a still better approximation to the 
Laplace transformation for these special functions, provided 7 > 1. 


The extension of the definition of 3** to general f(t) of exponential type 
meets with the difficulty that the series 


x KLE z Te (3.9) 


which seems to be the — way of defining 3** /, is in general divergent. 
In a purely formal way one finds that 


foe) yr-l co H=1 So ees 
Dp 43" Goa = 243 @aayr + 3a, SEW 4610) 
aris ea) 


r= 


at 


The first series in the right member converges and represents 3 /, by formula 
(2.5), but the second series is generally divergent (cf. [4], chapter 3, § 6). 

In order to arrive nevertheless at a definition of 3** f/f one may take 
recourse to the theory of asymptotic series. A power series 


not necessarily convergent, is called an asymptotic representation of a function 
g(x), in symbols 


x)~ - o, 2°) for 4 a+ 0:, 
“u=0 
if 


im a" | p(s) — 3, a =) eae aes KOie 
0 u=0 


Fo 


for all non-negative integers m. An important theorem in the theory of asymp- 
totic series [5] says that corresponding to any power series — convergent or not — 
there exists a function g(x), holomorphic for 0 < x S %p, of which it is an 
asymptotic expansion. The function g(x) is not uniquely determined, in fact it 
can be chosen in infinitely many ways, but all these functions differ by quan- 
- tities that are asymptotic to zero, i.e., by functions that are of order O(x™) 
for arbitrarily large m. Now let y(f, 2) be an otherwise unspecified holomorphic 


function of such that 

Wm Soa, EW ee (3.11) 
and define 3** / by 

3**f=3f+7(64), _ (3.12) 
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then (3.10) shows that 
Ah 


wn =1 
Bf 3, aed ap apm he gee (3.13) 
(Observe that the series in the right member is a power series in h.) 

Unless (0) happens to be zero, 3** f leads to an approximation of the 
Laplace transformation that is no better, as far as the order of magnitude is 
concerned, than 3* /. The former transformation is mentioned here, because 
it is used in [2]. Whether it leads to better numerical results than 3* in the 
applications made in [2] remains to be investigated. 


4. Approximate Solution of Linear Differential Equations with 
Constant Coefficients 


An ordinary initial value problem for an ordinary linear differential equation 
with constant coefficients and with a nonhomogeneous term w(t) that is entire 
and of exponential type can be solved, in a well-known way, by means of the 
Laplace transformation. This method leads to an expression for the Laplace 
transform Ly of the solution y of the problem which has the form 

P(s) 
1+ Q(s) ’ 
where P(s) and Q(s) are holomorphic at s = coand have a zero there. Actually, 
Q(s) is a polynomial in s~1, but this is immaterial here. In order to find y from 
(4.1) one may expand the right member into a series in powers of s~1, say 


Sy) == (4.1) 


Oye Sa oe (4.2) 
r=1 


from which one finds, by termwise inversion, 


lee) 


DS 2 eae aa (4.3) 


r=1 
This series may be evaluated directly or by means of Rule 1 of Section 2. 


An alternate method, based more directly on the theory of the z-transfor- 
mation, is sometimes used in electrical engineering. For its description let 


P(s)"= bsg Gs, “O(s\ 5s OES (4.4) 9 
Then ee af om 
pt) = £1 P(s) = Ja, (4.5) 
a) = £7.09) = 374, (4.6) 
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The functions £(¢), g() are of exponential type (cf. [6], page 62). Formula (4.1) 
can be written 


__*? 
AY asses Te eae (4.7) 
In view of the inequality (3.4) the right member of (4.7) can be approximated 
ag h 3s igs 
1+h3.9° oe 


This is an analytic function of e~** = z-!, which can be expanded, by long 
division, into a power series in e-’”, say 


h3 e snk. 
1 a q =h S5,¢ ; (4.9) 


Since the expression (4.8) is an approximation to  y, it is plausible that y,, is 
approximately equal to y( h). Because of the inequality (3.6) and the ana- 
logous inequality for 3** one may expect that the same procedure carried 
out with 3* or 3** will lead to still better approximatious. If this plausibility 
consideration is valid one has the following rule for the approximate calculation 
of y from (4.1). 


Rule 2. Expand P(s) and Q(s) in (4.1) in powers of s-1. Replace every s—" by 


hT 
W—Dr ty (3) 


divide by h and expand the resulting function of z in powers of 2-1. Then the 
coefficient of z-" is approximately equal to y(n h). Still better approximations are 


obtained 1f s-1 is replaced by h 
h (2) — 2? 
or if every s-" is replaced by 
hr of B; 
@ = ie pear | ) yl 


The second variant of this rule corresponds to the use of the operator 3*, 
the third to that of 3**. 

A proof of the validity of Rule 2 will be given here for the version based on 
the operator 3*. The proof for the ordinary z-transformation is analogous but 
simpler. The reasoning does not follow the preceding heuristic argument, but 
employs directly the fact that 


OE ee ot Da an! cy ek 


where f, = p(w h), In = (0 hh). 
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Let BP co 
T+ h3*¢ dy Vi Fae (4.11) 


Inserting the series (4.10) into (4.11) and multiplying by the denominator one 
obtains, after an interchange of summations, 


oo co n—1 
1 os —n __ 1 ~ on os gon 
L pet Steere (14 5 ha) Beet 4S (Tinton) 


Comparison of like powers of z~! leads to the recursion formulas 


2 ~ py = nm (1 ay +h qo) 5 Yo» (4.12) 
n—1 
by = (Lt ZAG) Fn th 3 Im tnm (> 0), (4.13) 


from which the y, can be determined successively. If yp is eliminated from 
(4.13) by means of (4.12), the formula (4.13) takes the form 


n—1 
i ~ i) a 
2 (14+ 5 ha) Moc (1 ba zh G0) Pn Poth D Vatnem Lee oe 
m=1 
(4.14) 
where the last summation is to be replaced by zero for m = 1. 
In order to compare ¥,, with y( h) = y, the equation (4.7) for y is written 
in the form Lf= Ly+ Ly Lq. Taking inverse Laplace transforms and — 
making use of the convolution theorem the integral equation . 


10+ | 0 att —7) (4.15) 


is obtained. Now set ¢= mh and replace the last integral by the trapezoidal — 
approximation formula and its remainder term, which is of the form h? w,, — 
where w, is a function of ¢= nh and h bounded for0 <h<h,0<t<T 
(h,, T finite but arbitrary). Then (4.15) becomes | 


1 es i | 

Po Va Pe RY, 1D) Ve, Geen tate eo (n > 0). ele 
m=1 

For n = 0 formula (4.15) becomes ) | 

Po = Mo - Gay | 


Eliminating y, from (4.16) by means of (4.17) one sees that (4.16) differs from — 
(4.14) by terms of the order O(h?). Hence, if one sets, for abbreviation, 


| 

' 

8, =n — n> (4.18) 
a 

{ 
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subtraction of (4.16) from (4.14) yields the equations 
1 n—1 
(+52 a) On +h S) Sy Gum =H dy, (mn =1,2,..), (4.19) 
m=1 


where ,, has the same properties as w,. 
In order to appraise 4,, let 


Then (4.19) implies 


n—1 


1 | 4 
[|1+ zh — tS tanl | 4, SEK 
I _ m=1 


For sufficiently small values of T and h there exists a positive number e inde- 
pendent of / such that the expression in brackets is not less than e, form < T/h, 
because the second term is the Riemann sum of an indefinite integral. Hence, 


in such an interval 
i3,,5 4,558 —. 

The last inequality proves that Rule 2, when carried out with the operator 
3*, leads to an approximation of order O(/?) for y(~h), » > 0 [but not for 
y(0)]. If the z-transformation itself is used the approximation is O(/) for ali 
mn = 0. 

The analysis of Rule 2 when the operator 3** is employed can be performed 
by the same method, but the details are tedious and will be omitted here. It 
turns out that the result is in general no better than if 3* is used, as far as the 
order of magnitude of the error is concerned. However, if the first k (k > 0), 
derivatives of #(¢) and g(#) vanish at f= 0, the operator 3** does lead to a 
better result. The error is then O(/?**). For the proof of this fact the integral in 
(4.15) must be approximated to the order O(h***) by means of the MacLaurin 
formula. Even in this case, however, it suffices to replace the first nonzero 
terms in the series (4.4) by the corresponding expression h 3** ¢"/(r — 1)! in 
applying Rule 2. For the later terms h 3 ¢"/(r — 1)! is an adequate substitution, 
if the error is judged by its order of magnitude. 


5. More General Differential Equations 


Boxer and THALER [2] propose an interesting procedure for the integration 
of linear differential equations with variable coefficients and for certain non- 
linear differential equations. The present section contains a brief discussion 
_ of the validity of this part of [2]. 
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The authors of [2] describe their algorithm for differential equations of 


the form 
m—-1 


ym + Se, 9” = vl), (5.1) J 
v=0 
where the g, may be any analytic functions of ¢, y, y’,..., y?~”. (vy is an 
abbreviation for d’y/dt.) The function y(t) must be entire and of exponential 


type. 
It is doubtful that anything is gained by using the operators 3* or 3** 


for such equations, rather than the simple z-transformation, since the procedure 
contains other approximate steps which are of order O(h), unless a considerable 
amount of additional labor is expended. The arguments will therefore be 
carried out for the z-transformation only. Furthermore, the discussion will 
be limited to the special case 


vy +py=ylt), p=olhy,’), (5.2) 


in order to simplify the manipulations. The treatment of the more general 
differential equation (5.1) requires no new ideas. 

In order to find an approximation to the solution of (5.2) with the initial 
values 


y(0) =%, (0) = By, (5.3) 9 


Boxer and THALER suggest the following algorithm. Let = (0, %, Bo), 
then the differential equation problem with constant coefficients 


Yo + Yo Vo= PE), Yo) =%» Yo(0) = Bo (5.4) 


can be approximately solved by Rule 2 of Section 2. To this end a certain 
fraction of the form 
3 Pl) 


aeons (5.5) 


where p(t) is an entire function of exponential type, has to be expanded accord- 
ing to powers of z~?. This expansion is best performed by long division, which 


yields successively approximate values for the numbers y(n h). In order to- 
obtain an approximation to the solution of (5.2) itself BoxER and THALER 
propose that, after having found y,(0) in this manner, g») be replaced by 


p(h, y(n), y'(A)) and that the next step of the division be carried out with 
the modified denominator. Then one continues the division, replacing gy at 
every step by (nh, y(nh), y'(h n)). Of course, y(n h) and y'(wh) are not 
known. In computational practice the substitution must therefore be performed 


with the nearest already available approximation to y(n h) and with some 


approximation to y'(m h) based on the previous calculations. 
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For the mathematical analysis of this method the following two simple 
facts are needed. 


1. Let /(¢) be of exponential type, and let /(0) = 0. Then 
23/H=3HE+h). (5.6) 


This is proved by the series of equalities 


2310) =F) Minh) A= YS f((n +1) "= BHM. 


n= 


_ 


2. A simple application of the convolution theorem shows that 
; 
s2 2 y(t) = a} (t — 2) w(t) de. (5.7) 
0 


Application of Laplace transformation to the problem (5.4) yields 
__ 2 Pl) 
ed caer + Okt 


where 
t 


Ppartevthrttaet= el (ta y0 dob By ttt): 
0 


Let 


t 
ne) = 36-3 (¢ — 1) p(t) trata, (5.8) 
0 
The result of the first step of the long division of (5.5) is a relation of the form 
ll saa wt") 5.9 


which can be regarded as the definition of the analytic function ti(27*). Since 

(0) = a by (5. 8), it follows from the definition of the z-transformation that 
Xo(O) = HX - 

Hence, (5.9) and (5.8) lead to 


t 


ne) = 23 le 1) y(t) dt + Bo! 


0 


+ 3 —2ahg3t- 


414 Wo.rcanc Wasow ZAMP 


Application of (5.6) changes this into 


mlz) = 3 fe — t) p(t +h) dv + Bolt + h) + % — dy h Po(t + 7 
ie (5.10) 


t 


-a{fe-ayermere pea, 


where 
0 


a4 = 9 + By h — ao Yoh? — |r ple + h) dt (5.11) 


Br = Bo— a0 poh + | ple+ H) de. (5.12) 


Let now @,, be some numbers such that, for all 7 in theintervalO SxnuhXT 
under consideration 


lpn— (nh, y (nh), y'(nh))| SKh (K= const) . (.13) 


If y(0), y(h,..., y(m — 1) h) are already known to within an accuracy of 
order h, such a number g, can always be found. This point will be discussed 
more fully below. 


The fraction = 
xi (2—*) 
1+h,p3t 


is the basis of the second step in the long division. It is a power series in 2~+ 
whose leading term is «,. It will now be shown that «, is an approximation to 
y(h), and that the repetition of the divisions leads to a sequence of numbers «, 
which approximates (7h). The precise statement to be proved is as follows. 


Theorem. Se @ = (0, %, Bo) and assume that the p,, n > 0, satisfy (5.13). 
With y(2-1) given by (5.8) let the functions y,(2-1) be defined recursively by 


Laila) Coe =1__ Xn+i(2—") 


Then there exist constants T, hy, M, independent of h and n but dependent on , 
yp, K such that ‘ 


lxn(0) —y(nh)| SMh (0S nh<T,0ShSh,), (5.15) 


where y(t) is the solution of (5.2), (5.3). 
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Proof: The proof consists in showing inductively that 


t 
ete) 7 KS fF (¢ a T) y(t +n h) at + Bn t+ ml (5.16) 
| 
where «, = x,,(0) and £, is a certain constant [cf, (5.19)| such that 
An=|%n— y(nh)|+|Ba—y'(nh)| <Nnk® (N=const, 0Snh< Te 
(5.02) 


For n = 0, the relation (5.16) has already been proved and (5.17) is obvious. 
Assume that these relations are true for n = k. By (5.14) 


Xnsr(2-1) = 2 An(2-1) — 2 ay — 2 Oy Qrh Bt. 


If y,(z~*) is replaced by the right member of (5.16) with n = k; a calculation 
analogous to that in the ei from n = 0 to m = 1 in (5.10) leads to 


Anvil? -3f/ (¢— 1) + (R+1)h) dt + Bryrt+ usp, (5.18) 
where m 
tensa = oe + Bult — on Gh? — [rp(t+ (b+ 1) h) dr, (5.19) 
—h 
0 
Bess = Be- On eh + [p(t + (R41) A) dt, (5.20) 
ou 
—¥) 


in generalization of (5.11) and (5.12). Equation (5.18) proves (5.16) for 
m=k-+1. 

In order to show that (5.17) is true forv=k+1,i.e., that p44, Bory are 
approximation to yz4, = y((k +1) h) and y’z4, = y/((k + 1) A), respectively, 
the given differential equation (5.2) for y must be used. One has, by TayLor’s 
formula, 


Vert = Vp thy tho, (R+1IAST, [a2 
with : 
eet (e) |, 5.22 
ete pees PD 1) (5.22) 
and from (5.2), 
(k+1)h (k+1)h . 
Vent = Va 7 / p(t, y(t), v’(t)) v(t) dt + A w(t) dt | 
kh kh (5.23) 


0 
= yi — hyp OUI, ny i) +f wE+ (+A) h) dt + Io, 
=F 
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with 
|», | S max 
=ts a 


Thanks to assumption (5.13) the relation (5.23) can be written 


0 


Yar = %— Bh eet [vlc + +1) h) det Wop, (5.25) 
—h 


where y¥# is uniformly bounded for k + 1 < T/h. 

Now subtract (5.21) from (5.19) and (5.25) from (5.20), form the absolute 
values of the two members in each of the two resulting equations and add 
them. By virtue of (5.13) the quantities | p, | are bounded for m S T/h, the 
bound depending on 9, py, K and T only. Furthermore, 


|a,| <NRAF+ max | y()], 
OstsT 
because of the assumed truth of (5.17) for n = k. In consequence of these two 


facts, of the inequality (5.22) and of the boundedness of v¥ there results, after 
a short calculation, the inequality 


Ina S In + AMht+t BNEM+CH, (5. 26) 


in which A, B and C are positive constants independent of N, k and h. They 
may depend on JT but they need not be increased, if T is decreased. The 


inductive assumption that (5.17) is valid for 1 = k shows that (5.26) implies — 


its validity for n = k + 1 provided 
ANH + BNRK+ Ch? aN fe. 
Since k h < T, the last inequality is certainly satisfied, if 
ANT+BNTA+CEN. 


This will be the case if N and T are chosen so that the inequalities 


N—C 


Nhe ee en ee 


a [y@) (4 9), 9’) ] | (5.24) | 


are true. This may require a reduction of the original size of T, but, as was - 


mentioned before, this can be done without changing A, B or C. This completes 
the induction from k to k + 1 and hence the proof of the theorem. 


In the light of the inequality (5.17) it is clear that an acceptable choice of 
the numbers @,, is 
Pa (n h, &n, Br) : 


The numbers 8, could be found successively from (5.20). Since the exact 


integration involved in the use of this formula may be inconvenient, a sequence ~ 
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Bx obtained by approximating the integral by a quadrature formula with an 
error that is O(h?) could be used instead. A still-simpler procedure consists in 
replacing y, in y by (Vn_1 — Vn—9)/h, which, because of (5.21) and the con- 
tinuous differentiability of y’(¢), changes m by an amount of the order O(h). 

The proof of the foregoing theorem shows that the method of Boxer and 
THALER is closely related to the simple stepwise integration procedure in which 


approximate values y*, y,* for y(mh) and y'(n h) are successively calculated 
from the equations 


Inti Vat Verh, yh. =" — oF @ (nh, y*, yi) h+ yp (whyh. 


These approximation formulas can easily be justified and they differ from 
(5.19), (5.20) by terms of order O(h?). Like the method of [2] they yield, in gen- 
eral, approximations of order O(), and they seem to be at least as simple. 
It must also be said that numerical solutions with this order of approximation 
are not considered very accurate in numerical analysis. Several fairly simple 
methods that yield higher precision are known. However, the order of magni- 
tude of the error term is not the ultimate criterion for the usefulness of a 
numerical method. Since a finer appraisal of the error is probably very difficult, 
the question is best decided by numerical experiments. 
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Zusammenfassung 


Einige numerische Anwendungen der Transformation 


3t= DY tah) 2 
n=0 


einer Funktion /(¢) (sogenannte z-Transformation) werden in dieser Arbeit auf 
ihre mathematische Giiltigkeit hin untersucht. Es zeigt sich, dass manche dieser 
Techniken bekannteren Differenzenmethoden sehr ahnlich sind. Die Arbeit 
bespricht ferner die Beziehung zwischen der Laplace-Transformation und der 
z-Transformation sowie zwei verbesserte diskrete Annaéherungen zur Laplace- 


_ Transformation. 


(Received: April 16, 1957.) 
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First Integrals of Dynamical Systems 


By OrENeE Bottema, Delft?) 


1. We consider a mechanical system with m degrees of freedom and the coor- 
dinates 41, 7s, --» » Ini it is holonomic, conservative and nongyroscopic; the kinetic 
energy T and the potential energy V do not contain ¢ explicitely. A set S of 
equations of motion is given by the m Lagrangian equations 

ah OF oT OV 
c= _— - == 
dt 04x 99% 09% 


If one knows m SX m mutually independent first integrals of the system 


FQ, 9) = & (iS 2 eee) (2) 


0. (1) 


it is a customary procedure to omit m equations of (1) and to solve the set S’ of x 
equations consisting of (2) and £,— 0-(k=1,2,...,7—m). The set Sis 
considered equivalent to S. 

Obviously the first integrals are consequences of the equations of motion and 
as such they are without doubt necessary conditions. It does not mean, however, 
that S’ is a set of sufficient conditions. A trivial example is the following. In 
the above case a first integral is given by the energy equation 


DEN (3) 
It is well known that it can be derived from (1) by the identity 


: d 
=r + Vs 


making use of the fact that T is a homogeneous quadratic function of g;. If we 
take the set S’ consisting of L, = 0 (k = 1,...,— 1) and (3), then a solution 
of it satisfies g, L, = 0; therefore it is not necessarily a solution of (1) in the 
case q, = 0. 


2. We give a somewhat more complicated example where the assumed 
equivalency of S and S’ could easily give rise to an erroneous result. 

Consider a double pendulum moving in a fixed plane under no forces. It 
consists of two homogeneous rods 4B (length J,, mass m,) and BC (length /,,” 


mass m,), A being a fixed point. If we take as coordinates the angles m, and @, 
of AB and BC with a fixed direction, we have 


1 . SARE 1 
P= 3 4 i+ B cos(oy — 92) o1 2 + = COR, (4) 
where 
ASSEN Gah Bae 1, 1 G 1 2 
3 1 » Me) '4 , ms ain = 3 Me 'e. (5) 


1) Technische Hogeschool. 
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For the moment of momentum of the system about A we have 


T= A 9, + Bcos(py — 2) (y+ 2) + C He. (6) 
Hence two first integrals of the system are 
1 : oe ‘| ; 
FL=qZ A vit Beos(g— o2) digst 7 CH= 4, (7) 
F, = A 9, + B cos (py — 2) (1 + 2) + C by = Cy. (8) 


Obviously gy, = 2/2 + wt, y= wt, where o is a constant, is a solution of the 
set (7, 8), with c,= (A + C) w?/2, c,= (A+ C) w. The conclusion would be 
that if for ¢ = 0 4B and BC are at right angles and have equal angular velocity, 
this situation is maintained during the motion. 


It is easily seen that the conclusion is wrong. The Lagrangian equations of 
motion are 


L, =A 9, + Bcos(9; — $2) $2 + Bsin(y, — v2) = 0, (9) 


L, =C G+ Bcos(y; — pe) $, — B sin(p, — 2) pj = 0. (10) 
Obviously 


Ed 
Sie yee P2=— wt 


is mot a solution of this set. If for ¢ = 0 


TU : ¥ 
Pa aF8 2? i= Po — OO, 
then 


3. From (7) and (8) it follows that 


G,= bapa A 9; 9 + Bcos(p — $2) (G1 G2 + %1 G2) — B sin (~, — 2) (11) 
x (p, — 2) Pr G2 + C 2 = Of 


Gus = A G, + Bcos(y, — %) (1 + G2) — Bsin(y, — $2) (Pt — 93) 


(12) 
+CG,=0. 
Obviously we have the following identities 
Gy= lit Pots, G,=1,+L,. 


Therefore G, = 0 and G, = 0 are consequences of L; = 0 and L, = 0. On the other 
hand ; ; : 
. (1 — G2) Ly = Gi — G2 Ge, (G1 — G2) La = —Gi t+ 1G. 

From G,= 0 and G, = 0 it follows (p, — Gs) Li = (G1 — 92) Ly = 0. Hence a 
solution of the set (7), (8) is not always one of (9), (10), a sufficient subsidiary 
condition being 9, + Pz. 
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Zusammenfassung 


Wenn man eine oder mehrere der Lagrangeschen Gleichungen durch erste 
Integrale ersetzt, dann sind die neuen Bedingungen zwar notwendig, aber nicht 
immer hinreichend. 


(Received: February 18, 1957.) 
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9° Assemblée du Comité suisse d’optique (CSO.) 
du 22 mai-1957 a Yverdon et Ste-Croix 


Le 22 mai 1957 le Comité suisse d’optique, section photonique, a tenu son 
assemblée annuelle A Yverdon et Ste-Croix. Au cours de cette séance six com- 
munications ont été présentées. Nous donnons ci-aprés le résumé de ces com- 
munications. , 


Reziprozitatsfehler bei photographischen Schichten. Von J. EcceErt, Lis 
rich). 

Nach der Reziprozitatsregel von BunseN und Roscoe soll die photoche- 
mische Wirkung nur von der Anzahl absorbierter Quanten, also von dem Produkt 
It aus Bestrahlungsstarke J und Bestrahlungszeit #, nicht aber von den Fak- 
toren selbst abhangen. Fiir die photographische Wirkung trifft diese Regel nur 
beschrinkt zu, wie besonders SCHWARZSCHILD feststellte. Diese Reziprozitats- 
fehler werden in Ziirich mit einem automatisch arbeitenden Gerat fiir den 
Bereich. von 7 Zehnerpotenzen und fiir verschiedene Entwicklungsarten ermit- 
telt, um das Verhalten der Entwicklungskeime topographisch und kinetisch zu 
untersuchen. 


Etudes de certaines propriétés optiques du corps solide. Par R. MERCIER 
et B. V. Vitroz, Lausanne’). 


Ces travaux vont prochainement étre publiés. 


Zur Physik und Technik des Kalten Spiegels. Von R. KrAut, Balzers?). 
Unter Kalten Spiegeln wird eine Anordnung einer Vielzahl von diinnen, 
interferenzfahigen Schichten verstanden, die das sichtbare Spektrum im wesent- 
lichen reflektiert, wahrend das infrarote ungehindert passieren kann. Der Aufbau 
erfolgt durch (A,)/4 ”)-Schichten mit méglichst grosser »-Differenz benachbarter 
Teilschichten. Infolge hdéherer Ordnungen treten im Ultravioletten weitere 
Reflexionsmaxima auf, die aber durch Absorption iiberlagert sind, wahrend im 
Sichtbaren und Langwelligen die Schichtsubstanzen absorptionsfrei sein miissen. 
Die Absorptionsfreiheit wird nach einem patentierten Verfahren dadurch 
erreicht, dass wahrend des Aufdampfprozesses eine Aufoxydation der zur Ver- 
dampfung gelangenden Suboxyde in die am héchsten abgesdttigte Sauerstoff- 


1) Photographisches Institut der ETH. 
2) BPUL. ! 
3) Geradtebauanstalt Balzers. 
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verbindung erfolgt, indem dem Molekiilstrom der zu verdampfenden Substanz 
ein Sauerstoffstrom mit einem Partialdruck tiberlagert wird, welcher den O,-Dis- 
soziationsdruck der abgesattigten Verbindung kompensiert. 

Da die verwendeten Schichtsubstanzen hohe chemische Bestandigkeit, 
Warmefestigkeit und grosse mechanische Harte besitzen, eignet sich der Kalte 
Spiegel auch als Oberflachenspiegel, so dass sowohl eine Schlifflache eingespart 
wie auch schlieren- und blasenhaltige, warmebestandige Glaser verwendet werden 
k6énnen. Das Reflexionsvermégen im Sichtbaren ist hdher als bei Ag- und Al- 
Spiegeln. Bei Verwendung Kalter Oberflachenspiegel in Verbindung mit Bogen- 
lampen tritt zudem bei Beschddigung der Oberflaiche durch Abbrand nur ein- 
facher Lichtverlust ein, wahrend bei riickseitig belegten Spiegeln infolge des 
Strahlenganges der Verlust doppelt so gross werden kann. Der von einem Kohle- 
bogen iiber einen Kalten Spiegel auf ein Vakuum-Thermo-Element abgebildete 
Krater ergab nur ein Viertel der Thermospannung verglichen mit dem iiber einen 
Al- oder Ag-Spiegel abgebildeten Krater. 


Vollkommene Objektive. Von R. Sterrier, Aarau‘). 

Im wesentlichen deckt sich der Inhalt des Referates mit den nachstehend 
aufgefiihrten Veréffentlichungen des Referenten: 
Uber die optische Abbildung von Flachen und Réumen, Optik 12, 529-543 (1955). 
Uber vadialsymmetrische optische Medien, Optica Acta 3, 101-103 (1956). 


Gesichtspunkte beim Bau eines universellen Spektral-Photometers. Von 
E. LtscHer, Lausanne’). 

Die Hauptanforderung, die an einen Monochromator gestellt wird, ist die, dass 
die Strahlungsleistung so wenig wie méglich verringert werden soll. Man wird 
also den Lichtleitwert®) so gross als méglich wahlen: 


A A. spt ‘ 
f 
Dabei bedeuten: AA Bandbreite; A Auflésungsvermégen; / Spalthohe; d Durch- 
messer; f Brennweite des Kollimators. 

Da ein méglichst grosser Spektralbereich erwiinscht ist, wurden als Kollima- 
toren spharische Spiegel anstelle von Linsen verwendet. Die Spiegel sind mit einer 
aufgedampften Aluminiumschicht versehen, damit auch im UV gearbeitet werden 
kann, da Silber bei 3600 A einen steilen Abfall des Reflexionsvermégens’) zeigt. 

Als Monochromator-Montierung eignet sich die Anordnung nach CZERNY- 
TurRNER®’) gut, da sowohl mit Prismen als auch mit Reflexionsgitter gearbeitet 
werden kann. Da eine grosse Spalthéhe ausgeniitzt wird (25 mm), muss der Se- 
kundarspalt bei Verwendung von Gittern wegen der Kriimmung der Spektral- 
linien bei Prismenmonochromatoren ausgewechselt werden kénnen. WIHLM’®) 
konnte zeigen, dass bei dieser Montierung Asymmetriefehler (Koma) der spha- 
rischen Spiegel vermieden werden und ferner bei nicht zu grossem Offnungs- 
verhaltnis (bei unserer Montierung 1:5) die spharische Aberration gegeniiber 
dem Astigmatismus vernachlassigt werden darf. 


4) Kern AG. 

5) EPUL und optisches Laboratorium Metrohm. 

6) G. Hansen, Optik 1, 227 (1946). 

7) G. B. Sasrne, Phys. Rev. 55, 1064 (1939). 

8) M. Czerny und A. F. Turner, Z. Phys. 61, 792 (1930). 
®) H. Wiuum, Rev. Opt. 33, 425 (1954). 
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Bei Monochromatoren fiir das Schumann-Gebiet muss die Anzahl der 
reflektierenden Flachen wegen dem geringen Reflexionsvermogen aller bisher 
untersuchten Materialien méglichst klein gehalten werden. Wir wahlten deshalb 
fiir dieses Gebiet eine Paschen-Runge-Gittermontierung, bei der das Gitter die 
einzige reflektierende Flache bildet. 

Der grosse Spektralbereich bedingt die Verwendung verschiedener Detek- 
toren, die sukzessive eingeschaltet werden kénnen. Im UV und Sichtbaren eignen 
sich Photomultiplier am besten. Thermoelemente werden fiir das IR montiert. 
Damit die Schwierigkeiten der Nullpunktstabilitat vom Gleichstromverstarker 
vermieden werden, wurde der Lichtstrahl vor dem Eintrittspalt zerhackt und 
mit einem Wechselstromverstarker gearbeitet. Die Zerhackerfrequenz muss 
wegen der Tragheit der Thermoelemente niedrig (11 Hz) gewahlt werden. 

Damit Gaszellen (1,5 m Absorptionslange), Absorptionsgefasse im Schmelz- 
ofen, verwendet werden kénnen, wurde der eigentliche Absorptionsraum moglichst 
gross ausgebildet. Dies bedingt allerdings eine erneute Kollimation nach dem 
Messgefass vor den Empfangern. 

Die universelle Verwendung eines Spektralphotometers bedingt auch eine 
Umschaltung auf verschiedene Lichtquellen (UV-Hg-Lampe, W-Lampe, Nernst- 
Brenner, Flamme, Lyman-Lampe). Die Anordnung dieser verschiedenen Licht- 
quellen ist so getroffen, dass diese auf einem Halbkreis mit einem spharischen 
Spiegel im Zentrum stehen. Der Hohlspiegel bildet die jeweilige Lichtquelle auf 
den Eintrittsspalt ab. 

Eines der Hauptprobleme beim Bau von Spektralphotometern bildet die 
méglichst weitgehende Unterdriickung von Streulicht. Dies kann nur durch 
geeignete Anordnung von Blenden erreicht werden. Fiir nichtmoduliertes Licht 
ist der Empfanger ohnehin unempfindlich. 


Ein neuartiges Gerat fiir die Fernseh-Abtastung von Filmen. Von TH. 
STUTZ, Ziirich?®), 

Im Auftrage der Schweizerischen PTT baute die Abteilung fiir industrielle 
Forschung des Institutes fiir technische Physik der ETH ein weiteres Modell 
eines Fernseh-Filmabtasters fiir 16-mm-Schmalfilm, nachdem der erste Prototyp 
seit dem Dezember 1953 im Fernsehstudio Bellerive dauernd im Betrieb steht. 
Es handelt sich um einen Leuchtpunktabtaster mit einer neuartigen Konstruk- 
tion zum optischen Ausgleich der Bewegung des mit beliebiger Geschwindigkeit 
durchlaufenden Filmes. 

Grundsatzlich ist in jedem Filmabtaster ein Film-Wiedergabegerat mit einer 
Einrichtung zur Fernsehaufnahme zusammengebaut. Leider kann man aber von 
den verschiedenen Méglichkeiten der Filmwiedergabe und der Fernsehaufnahme 
nicht die einfachsten und bewahrtesten kombinieren. Wollte man zum Beispiel 


mit dem seit friihesten Zeiten beliebten Leuchtpunktabtaster einen auf einem - 


normalen Kinoprojektor vorgefiihrten Film iibertragen, so liesse sich die Dunkel- 
pause nicht auf den durch die Fernsehnorm vorgeschriebenen Wert von weniger 
als 2 ms reduzieren. 

Eine Betrachtung aller Méglichkeiten zeigt, dass konventionelle Filmprojek- 
toren eine speichernde Fernseheinrichtung oder dann der einfachere Leucht- 
punktabtaster einen der heute nicht mehr gebrauchlichen Filmprojektoren ohne 
Dunkelpause bedingen. Bei diesen nichtintermittierenden Projektoren mit 
gleichférmig durchlaufendem Film ist es bei einem einfachen Verhaltnis der 


10) Abteilung fiir industrielle Forschung des Institutes fiir technische Physik der ETH. 
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Filmbild- und der Fernseh-Ablenkfrequenz moglich, die Filmbewegung zur Bild- 
zerlegung mitzubeniitzen. Die einzige Lésung, die gleichzeitig alle Film- und 
Fernsehnormen bewaltigt, ist der optische Ausgleich der Filmbewegung. Hier 
wird das Bild des gleichférmig bewegten Filmes mit optischen Mitteln stillgesetzt, 
so dass die Fernsehiibertragung vom Filmtransport unabhangig ist. 

Die konstruktive Verwirklichung von Projektoren ohne Dunkelpause bietet 
etwas andere Schwierigkeiten als diejenige von intermittierenden. So muss zum 
Beispiel die Optik dem Abstand zweier aufeinanderfolgender Bilder angepasst 
werden, und auch die Filmtransporteinrichtung muss den Schrumpfungszustand 
des Filmes beriicksichtigen. Mangel in der Bildfeldausleuchtung oder eventuelle 
Unscharfen sind besonders unangenehm, weil sie sich gegeniiber dem Bildinhalt 
bewegen. 

Fiir das hier zu beschreibende Gerat kam nur ein Leuchtpunktabtaster in 
Frage, da dieser eine besonders hohe Bildscharfe aufweist sowie keine heiklen 
Einstellungen und beim Farbfernsehen keine Deckungsprobleme kennt. Die 
Filmwiedergabe geschieht mit optischem Ausgleich der Filmbewegung, weil das 
Gerat auch imstande sein muss, verschieden normierte Amateurfilme zu iiber- 
tragen. Damit auch die Filmwiedergabe keine besonderen Anspriiche an die 
Bedienung stellt, soll die Beriicksichtigung der Filmschrumpfung vollautoma- 
tisch erfolgen. 

Das Prinzip der neuen Lésung ist in Figur 1 dargestellt. Darin sind von links 
nach rechts abgebildet: der Leuchtschirm einer Kathodenstrahlréhre, eine den 
Leuchtschirm ins Unendliche abbildende Linse, ein horizontalachsiges Blenden- 
rad, zwei durch eine Vertikalebene getrennte Halften eines Objektives, die den 
auf dem Leuchtschirm einen Fernsehraster schreibenden Leuchtpunkt auf je 
ein Bild des Filmes abbilden, sodann folgt hinter dem Film eine Kondensorlinse, 
die den gesamten Lichtstrom auf die dahinter folgende Photokathode sammelt. 
Zwei Greifer sind mit je einer der beiden Objektivhalften starr zusammengebaut, 
und diese beiden Gebilde sind langs einer Parallelen zur Filmftihrung beweglich. 
Zwei Kurvenscheiben bewegen jeweils das eine der beiden Halbobjektive gleich- 
férmig abwarts, wahrend das andere sich etwas rascher vom vorhergehenden zum 
nachstfolgenden Filmbild aufwartsbewegt, so dass sich wahrend einer gewissen 
Zeitspanne, die in Figur 1 festgehalten ist, beide Objektive im Abstand eines 
Filmbildes mit dem Film bewegen. Wahrend dieser Zeit wird durch das erwahnte 
Blendenrad die Abbildung vom jeweils unteren auf das obere Halbobjektiv tiber- 
blendet. Dabei ist das Blendenrad so ausgebildet, dass es von der einen Objektiv- 
pupille genau die gleiche Flache freigibt, wie es von der anderen zudeckt, so dass 
sich die Helligkeit wahrend der Uberblendung nicht andert. In Figur 2 ist die 
Objektivbewegung schematisch dargestellt und die Uberblendung durch die 
schraffierte Flache angedeutet. Die Greiferspitzen, die jeweils bei der Abwarts- 
bewegung in ein Perforationsloch eingreifen, sind so ausgebildet, dass sie den 
Film nur abwartsstossen kénnen. Bei der Aufwartsbewegung gegen den Film 
federn die Greiferspitzen zuriick. 

Da der Film nur von den Greifern transportiert wird und daher genau die 
gleiche Bewegung macht wie die jeweils aktive Objektivhalfte, wird in jedem 
Zeitmoment der auf dem Leuchtschirm zeilenweise geschriebene Fernsehraster 

-relativ stillstehend auf ein Filmbild abgebildet. Unabhangig von der Absolut- 
bewegung des Films durchleuchtet der Leuchtpunkt das Filmbild Punkt fiir 
Punkt. Der die Photozelle erreichende Lichtstrom entspricht jeweils der Hellig- 

 keit des betreffenden Bildpunktes, und der elektrische Strom in der Photozelle 

stellt das zu iibertragende Fernsehsignal dar. 
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Nur wahrend der in Figur 1 dargestellten Uberblendungsphase muss der 
Objektivantrieb den Schrumpfungszustand des Filmes beriicksichtigen. Zu 
diesem Zweck wird durch eine nicht dargestellte Einrichtung standig die Film- 
lange pro Bildteilung gemessen und nach dem Resultat die wirksame Lange der 


Figur 1 


' Prinzipschema. 


beiden in Figur 1 unten dargestellten Hebel verandert, die mit ihren vorderen 
Enden auf den Kurvenscheiben laufen und die Objektive vertikal bewegen. Zu 
dieser automatischen Langenmessung wird wegen unvermeidlicher zufalliger St6- 
rungen eine integrierende Anordnung nach folgendem Prinzip verwendet: Vor 
dem Bildfenster wird der Film durch eine ungezdhnte Rolle transportiert, die 
von der Kurvenscheiben-Welle aus tiber ein Regelgetriebe angetrieben wird, Nur 
wenn die Einstellung des Regelgetriebes der Filmlange entspricht, bleibt die Lange 
der Filmschlaufe iiber dem Bildfenster konstant. Eine Veradnderung dieser Schlaufe 


verstellt das Regelgetrieb ; . 5 : of 
Sih, AG gelgetriebe und damit das auf den Objektivantrieb uibertragene 


Vol. VIII, 1957 Varia — Miscellaneous = Divers 425 


Figur 2 
Die Bewegung der beiden aus je einer Objektivhalfte und einem Greifer bestehenden Gebilde. 


Neben den erwahnten Vorteilen aller Leuchtpunktabtaster weist das neue 

Gerat folgende Vorteile auf: 

1. Die wirksame Offnung der Optik ist durch keine konstruktiven Gegeben- 

heiten begrenzt. 

Es wird nur die minimale Zahl optischer Elemente verwendet. 

. Die Anordnung erzeugt beim Uberblenden keine Helligkeitsschwankungen. 

. Es entsteht relativ wenig Streulicht. 

Die Beriicksichtigung der Filmschrumpfung geschieht automatisch und ohne 

optische Kompromisse. 

. Durch den Greiferantrieb wurde eine Filmtransportgenauigkeit méglich, die 
derjenigen intermittierender Projektoren mindestens ebenbiirtig ist und keine 
libertriebenen Prazisionsanforderungen an die Antriebsmechanik stellt. 
Zum Schluss sei allen jenen gedankt, ohne die das Gerat nicht hatte gebaut 

werden k6nnen, allen voran dem Leiter unseres Institutes, Herrn Prof. E. Bau- 

MANN, sodann Herrn D. DIEBOLD, dipl. El.-Ing., und der Firma Kern in Aarau, 

als der Herstellerin der speziellen Optik. 


ab wh 
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TH. Stutz, Ein neues Verfahren zur Fernseh-Abtastung von Filmen, Techn. Mitt. 
Pro 2Nt. tL (1954). 


Rechenautomat 704IBM 


Die International Business Machines Corporation hat kiirzlich in Paris einen 
elektronischen Rechenautomaten vom Typus 704 in Betrieb genommen. Es 
handelt sich um den gegenwartig leistungsfahigsten Rechenautomaten in Europa. 

Die IBM stellt in grossziigiger Weise die Summe von ungefahr einer Million 
Schweizer Franken in der Form von Benutzungszeit der 704 der europdischen 
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Wissenschaft zur Verfiigung. Fiir die Schweiz macht dies bis Jahresende etwa 
51% Stunden Rechenzeit auf diesem Automaten aus. : 

Forschungsinstitute der Hochschulen oder eventuell der Industrie konnen 
von dieser Méglichkeit Gebrauch machen. Diesbeziigliche Antrage sind zu richten 
an den Unterzeichneten (Prof. Dr. E. St1EFEL, Eidgenéssische Technische Hoch- 
schule, Ziirich) oder an die IBM-France (5, Place Vendome, Paris) direkt. Sie 
miissen enthalten: 

1. Eine Schilderung des vorgeschlagenen wissenschaftlichen Problems. 

2. Eine Beschreibung der anzuwendenden numerischen Methode fiir das auto- 
matische Rechnen. 

3. Eine Schatzung fiir die beanspruchte Rechenzeit auf der 704. ~ 

Die Lésung des Problems auf der 704 muss von dem Antragsteller oder von 
einem durch ihn bezeichneten Bevollmachtigten vorbereitet, programmiert und 
in Paris durchgefiihrt werden. Auskunft iiber die Maschine 704 und die zuge- 
hérige Programmierungstechnik erteilt die Extension Suisse der IBM (Talacker 
30, Ziirich). E. Stiefel 


Arbeitstagung iiber die Ausbildung von Ingenieuren und Mathemati- 
kern in numerischer Mathematik unter Beriicksichtigung der elek- 
tronischen Rechenanlagen 


In der Zeit vom 27. bis 29. Mai 1957 fand in Miinchen am Mathematischen 
Institut der Technischen Hochschule eine Arbeitstagung tiber die Ausbildung 
von Ingenieuren und Mathematikern in numerischer Mathematik unter Bertick- 
sichtigung der elektronischen Rechenanlagen statt, an der Deutschland, Frank- 
reich, Osterreich und die Schweiz durch Fachexperten vertreten waren. 

In Referaten und Diskussionen kam tibereinstimmend die Ansicht zum 
Ausdruck, dass die modernen Grossrechenanlagen einen massgebenden Einfluss 
auf die Mathematik als Ganzes, ihren Anwendungsbereich in anderen Wissen- 
schaften und ihre Auswirkungen in Technik und Wirtschaft ausiiben. Daraus 
ergeben sich insbesondere weitgehende Riickwirkungen auf die Gestaltung des 
mathematischen Unterrichts in allen Stufen. 

Die Arbeitstagung empfiehlt daher ihren Teilnehmern, sich in ihren Landern 
fiir die Erfiillung folgender Forderungen einzusetzen: 

1. Staérkere Beriicksichtigung des numerischen Rechnens und seiner mit den 
modernen Rechenanlagen verkniipften Weiterentwicklung im mathemati- 
schen Unterricht der Technischen Hochschulen und Universitaten, wobei 
weniger an eine Vermehrung der Stunden als an eine der neuen Entwicklung 
angepasste Auswahl des Stoffes gedacht wird. 

2. Es wird empfohlen, das Fach «Numerische Mathematik» in die Abschluss- 
priifung fiir Mathematiker und Physiker aufzunehmen und auch den Inge- 
nieuren in der Abschlusspriifung Gelegenheit zu geben, Kenntnisse auf diesem 
Gebiet nachzuweisen. 

3. Auch an den mittleren und hdheren Schulen ist eine Reform und Vertiefung 
des Unterrichts im numerischen Rechnen erforderlich, welche hauptsachlich 
durch eine richtige Ausbildung der Lehramtskandidaten und eine entspre- 
chende Orientierung der bereits tatigen Lehrer erreicht werden muss. 

Es sei auch darauf hingewiesen, dass sich durch den zunehmenden Einsatz 
der modernen Rechenanlagen in Wissenschaft, Technik und Wirtschaft erfolg- 
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versprechende Berufsaussichten fiir die Absolventen mittlerer und hoherer 
Schulen ergeben. 

Die Teilnehmer der Tagung sind iibereingekommen, zu gegebener Zeit neuer- 
dings zusammenzutreten, um iiber die in den einzelnen Lindern inzwischen 
getroffenen Massnahmen zu berichten-und weitere Erfahrungen auszutauschen. 

CREMER, SAUER, STRUBECKER, WALTHER (Deutschland) 
DuRAND, KuNTZMANN (Frankreich) 

INZINGER (Osterreich) 

STIEFEL (Schweiz) 


ee ee eee OO eee 
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Plane Waves and Spherical Means. Applied to Partial Differential 
Equations. Von Fritz JOHN (Interscience Publishers, New York 1955). 172 S., 
10 Fig.; $4.50. 

Das Gebiet der partiellen Differentialgleichungen ist trotz seiner grossen 
praktischen Bedeutung auch heute noch voll ungeléster Probleme, und die Zahl 
der sich in den verschiedenen Zeitschriften damit befassenden Arbeiten ist beacht- 
lich. Fritz JOHN gibt in seinem Buch eine Zusammenfassung einiger Resultate, 
die sich unter gewissen einheitlichen Gesichtspunkten, namlich mittels elemen- 
tarer Identitaten zwischen Integralmitteln, gewinnen lassen. In einem ersten Teil 
wird die Zerlegung beliebiger Funktionen in Funktionen vom Typus ebener 
Wellen (Radon-Transformation) betrachtet und auf das Anfangswertproblem 
homogener hyperbolischer Gleichungen mit konstanten Koeffizienten sowie auf 
die Konstruktion von Fundamentall6sungen im Kleinen von linearen elliptischen 
Gleichungen und Systemen mit analytischen Koeffizienten angewendet. In 
einem gewissen Sinne eine Verallgemeinerung bedeutet die Heranziehung von 
spharischen Mitteln, die im anschliessenden Teil erfolgt und beniitzt wird, um die 
Differenzierbarkeit der Lésungen von linearen und nichtlinearen elliptischen 
Gleichungen und Systemen mit gentigend regularen Koeffizienten sowie die 
Analytizitat der Lésung im Falle analytischer Koeffizienten zu beweisen. 

Die Darstellung ist sehr iibersichtlich und so elementar als méglich gehalten. 
Das Buch diirfte alle, die mit partiellen Differentialgleichungen zu tun haben, 
interessieren. E, Roth-Desmeules 


Continuous Transformations in Analysis. Von T. Rapo und P.V. 
REICHELDERFER ([(Grundlehren der mathematischen Wissenschaften, Bd. 75] 
Springer-Verlag, Berlin, Géttingen, Heidelberg 1955). 442 S., 53 Abb.; DM 59.60. 

Nach Uberlegungen von Banacu, ViTaALt und andern lassen die fundamenta- 
len Begriffe der beschrankten Variation, der absoluten Stetigkeit und der Ablei- 
tung einer reellen Funktion mit mehreren Variablen eine geometrische Deutung 
zu, und zwar im Bereich der Transformationen in einem n-dimensionalen eukli- 
dischen Raum, die durch die reellwertigen stetigen Funktionen gestiftet werden. 
Das vorliegende Werk setzt sich zum Ziel, die topologischen und masstheore- 
tischen Grundlagen der Theorie der reellen Funktionen von diesem Standpunkt 
aus systematisch darzustellen. Die erforderlichen Hilfsmittel aus der Topologie 
und der Gruppentheorie werden kurz zusammengestellt, ebenso, jedoch bedeutend 
ausfiihrlicher, die Begriffe der algebraischen Topologie. Obwohl die Verfasser in 
erster Linie die allgemeine Theorie im u-dimensionalen Raum bieten, werden 
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noch die Falle n = 1 und ” = 2 gesondert behandelt. — Die grundlegende Idee der 
stetigen Transformationen findet beispielsweise in der Variationsrechnung An- 
wendung. Das Buch wendet sich — im Hinblick auf die Darstellung — vorwiegend 
an Mathematiker und Studierende der Mathematik, doch spielen die stetigen 
Transformationen natiirlich in der angewandten Mathematik ebenfalls eine Rolle, 
so dass sich auch etwa die Physiker theoretischer Richtung dafiir interessieren 
werden. E. Roth-Desmeules 


Photoconductivity Conference. Herausgegeben von R. G. BRECKENRIDGE, 
B. R. RussEtt, und E. E. Haun (John Wiley & Sons, New York 1956). 653 S., 
210 Abb.; $ 13.50. 

Im November 1954 trafen sich auf dem Gebiete der Photoleitung fiihrende 
Wissenschafter aus verschiedenen Landern in Atlantic City, um die standig an 
Wichtigkeit zunehmenden Probleme der Photoleitung zu diskutieren. Als Resul- 
tat dieser Diskussionen liegt heute ein dreissig Vortrage umfassendes Buch vor 
uns. Trotzdem ein solcher Konferenzbericht nicht die Einheitlichkeit eines Lehr- 
buches besitzen kann, gelang es den Herausgebern durch geschickte Gliederung 
der behandelten Themen, iiber die blosse Berichterstattung hinaus eine auch dem 
Nichtfachmann zugdngliche Einfiihrung in das vielschichtige Gebiet der Photo- 
leitung zu geben, ohne dass dabei die Aktualitat der Vortrage gelitten hat. 

Im 1. Teil wird in zwei Arbeiten die phanomenologische Theorie der Photo- 
leitung behandelt. Der 2. Teil befasst sich in sechs Vortraégen mit den kompli- 
zierten Lichtabsorptionsprozessen in Kristallen. Die sechs Arbeiten des 3. Teils 
setzen sich mit dem Verhalten der durch Lichtabsorption befreiten Ladungstrager 
in Festk6rpern auseinander. Im 4. Teil werden die Eigenschaften verschiedener 
Photoleiter diskutiert, und schliesslich folgen im letzten Teil neun Referate iiber 
neueste experimentelle Ergebnisse an verschiedenen Photoleitern. 

Das Buch erganzt in hervorragender Weise die recht sparlich vorhandene 
Literatur tiber Photoleiter. E. Mooser 


Science and Information Theory. Von LEon BrILLouin (Academic Press 
Inc., New York 1956). 320 S., 74 Fig.; $ 6.80. 

In zwanzig Kapiteln werden in sehr konzentrierter Form die wesentlichen 
Aspekte der Informationstheorie und ihrer Anwendung in verschiedenen Gebieten 
der Wissenschaft dargestellt. Der formalen Analogie des Begriffes Information 
mit demjenigen der negativen Entropie der Thermodynamik wird eine eingehende 
Betrachtung gewidmet und werden diejenigen Falle betrachtet, bei denen sich 
eine tatsdchliche Ubereinstimmung ergibt. Der Begriff Negentropie wird aber 
auch dort angewandt, wo kein Zusammenhang mit der Thermodynamik besteht, 
namlich in der Ubertragung von Sprache in Laut und Schrift. 

Ausgehend von der Definition der Information nach SHANNON wird die An- . 
wendung generell diskutiert und am Beispiel der englischen Sprache erldutert. 
Das Prinzip der Kodierung, die Kapazitat eines Ubertragungskanals, allgemeine 
Kodierungsprobleme, fehleranzeigende und -korrigierende Kode folgen. Mit 
einem Kapitel iiber Signalanalyse mittelst Fourier-Methoden und Zeitauslese 
und einem weitern die thermodynamischen Prinzipien zusammenfassend, sind 
damit alle Grundlagen zusammengestellt. Thermische Agitation, Brownsche 
Bewegung, Thermisches Rauschen in elektrischen Kreisen (Nyquist-Formel) 
finden eine kurze Darstellung. Es folgen Betrachtungen, welche dem Autor be- 
sonders zusagen, namlich das Negentropie-Prinzip der Information, das Problem 
von Maxwells Damon, die Anwendung des Negentropie-Prinzips der Information 
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in der Physik und Beobachtung und Information. Im letzten Viertel des Buches 
werden noch folgende interessante Themen behandelt: Informationstheorie in 
Zusammenhang mit der Unbestimmtheitsrelation und physikalische Grenzen der 
Beobachtung; das Negentropie-Prinzip der Information in der Ubertragungs- 
technik; Schreiben, Drucken und Lesen; das Problem des Rechnens (Rechen- 
maschinen) und zum Schluss etwas iiber Information und Organisation und iiber 
Probleme, welche ausserhalb des Bereichs der Informationstheorie liegen. 

Das vorliegende Werk ist heute das einzige, welches eine umfassende und 
konzentrierte Darstellung der Informationstheorie mit ihren vielfachen Anwen- 
dungen gibt. Auf die Originalaufsatze wird mit Fussnoten hingewiesen. Das Buch 
entspricht einem wirklichen Bediirfnis, und wir kénnen dem Verfasser fiir seine 
Arbeit sehr dankbar sein. H. Weber 


Handbuch der Physik. Herausgegeben von S. Friccr. Bd. 14; Kéilte- 
phystk I (Springer-Verlag, Berlin 1956). 349 S., 215 Fig.; DM 72.—. 

Dieser erste der zwei Bande iiber Kaltephysik enthalt fiinf englisch geschriebene 
Artikel. 

J. G. Daunts Bericht tiber die Herstellung tiefer Temperaturen — bis Wasser- 
stofftemperaturen — bildet eine fiir Physiker dusserst niitzliche Darstellung der 
Gasverifliissigungsprobleme. Die wichtigsten Methoden sind so beschrieben, dass 
der Tatsache, dass Physiker und Kaltephysiker keine Kalteingenieure sind, 
Rechnung getragen wird. S.C. CoLztins Diskussion tiber Heliumverfliissiger, 
deren industrielle Herstellung er massgebend gefdrdert hat, enthalt eine inter- 
essante Ubersicht iiber moderne Methoden. D. K. C. MacDona tp berichtet in 
einfacher Form iiber die Phanomene der elektrischen Leitfahigkeit bei tiefsten 
-Temperaturen. Die Beschreibung experimenteller Methoden seines Laborato- 

-riums diirfte fiir viele niitzlich sein. K. G. KtEmMENsS bearbeitet die Warmeleit- 
fahigkeiten. Angesichts der vielen experimentellen Untersuchungen der letzten 
Jahre ist seine Zusammenfassung der Resultate vom theoretischen Standpunkt 
aus sehr begriissenswert. P. H. KEEsom und N. PEARLMAN berichten iiber spe- 
zifische Warmen. 

Die Hauptschwierigkeit bei der Herausgabe lag aber offenbar in der Platzbe- 
schrankung. Der Leser hatte es begriisst, wenn die Aufsdtze noch ausfihrlicher 
waren. J. L. Olsen 


Verformung und Fliessen des Festkérpers. Kolloquium Madrid, 
26.-30. September 1955. Herausgegeben von R.GRAMMEL (Springer-Verlag 
Berlin, Géttingen und Heidelberg 1956). 324 S., 188 Abb.; DM 37.50. 

Das Buch enthalt die Vortrage und Diskussionsbeitrage vom Madrider Kol- 
~ Joquium 1955, das von der Internationalen Union fiir Theoretische und Ange- 
wandte Mechanik unter dem Vorsitz von R. GRAMMEL mit dem Zweck durch- 
gefiihrt worden ist, eine engere Zusammenarbeit zwischen Mechanikern und 
Festkérper-Physikern anzubahnen. Es gliedert sich in die folgenden Themata: 
1. Versetzungen und Plastizitat; 2. Nichtlineare Elastizitat und Vermischtes; 


3. Viskoelastizitat und Relaxation. as , 
Von den insgesamt 33 Beitragen seien im folgenden nur die allgemeinen 


Vortrage genannt: 

1. G.I. Taytor, Strains in Crystalline Aggregate. A. H. COTTRELL, Disloca- 
cations in Crystals. A. SEEGER, Neuere mathematische Methoden und phystkalische 
Ergebnisse zur Kristallplastizitaét. P. G. HopGE jr., The Theory of Piecewise Linear 


Isotropic Plasticity. 
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2. H. KaupERER, Ein nichtlineares Elastizitdtsgesetz ; Aufbau und Anwendungs- 
modglichkeiten. L. R. G. TRELoAR, Large Elastic Deformations in Rubberlike 
Materials. Ses 

3. M.A. Brot, Variational and Lagrangian Methods in Viscoelasticity. 


F. Scuuttz-Grunow, Viskositét und Zeitwirkungen im nichtlinearen Bereich. 
H. Ziegler 


Prazisionsmessungen von Kapazitaten, Induktivitaten und Zeit- 
konstanten. Von E. BLecuscumiptT. Band I (Friedr. Vieweg & Sohn, Braun- 
schweig 1956). 166 S. 84 Abb.; DM 11.80. 

In der Sammlung «Verfahren und Messkunde der Naturwissenschaft» ist 
soeben die zweite neubearbeitete Auflage der Monographie von ERICH BLECH- 
SCHMIDT erschienen, welche die vielfach zerstreuten Mitteilungen tiber die Verfah- 
ren zur Prazisionsmessung von Kapazitaten und dielektrischen Verlusten in tber- 
sichtlicher und knapper Form zusammenfasst. Der véllig neue Teil iber Induk- 
tivitatsmessungen wird als Heft 14 II in Kiirze erscheinen. Das kleine Werk 
enthalt als Einleitung eine Ubersicht iiber die Begriffe der Kapazitat der 
Dielektrizitats-Konstante sowie der Einheiten. Nach einer Diskussion der 
Kapazitaétsnormale werden die Messmethoden und die verschiedenen Bricken- 
anordnungen eingehend behandelt. 

Die vorliegende kritische Zusammenstellung diirfte vor allem fiir Physiker 
und Ingenieure von Nutzen sein, erlaubt sie doch eine rasche und zuverlassige 
Orientierung tiber die zahlreichen Faktoren, welche die Genauigkeit solcher 
Messungen beeinflussen k6énnen. E. Baldinger 


Advances in Electronics and Electron Physics. Von L. Marton. Band 6 | 
(Academic Press Inc., New York 1954). 538 S., 203 Fig.; $11.80. 

Die von L. Marton herausgegebene, jahrlich erscheinende Sammlung von 
zusammenfassenden Berichten tragt neuerdings den erweiterten Titel «Advances 
in Electronics and Electron Physics». Damit soll deutlich gesagt sein, dass auch 
rein wissenschaftliche Belange aus der Physik des Elektrons zur Sprache kom- 
men, Fragen also, die ausserhalb des Rahmens dessen stehen, was heute land- 
laufig unter Elektronik verstanden wird. Dies geht auch fiir Band 6 aus den 
Titeln der Beitrage hervor (Beitrag 2 und 8): 

A.B. PrpparD, Metallic Conduction at High Frequencies and Low Temperatures ; 
E. ABRAHAMS, Felaxation Processes in Ferromagnetism; J.Smit and H.P. J. Wijn, 
Physical Properties, of Ferrites; H. F.Ivey, Space Charge Limited Currents; 
W.M. WEBSTER, A Comparison of Analogous Semiconductor and Gaseous Electro- 
nics Devices; M.E.Hatne, The Electron Microscope — A Review; R.G.E.HUutTTER, 
Traveling-Wave Tubes; J. VAN DEN HANDEL, Paramagnetism. 

Im Grunde genommen hat sich jedoch an der Art der Auswahl nichts geandert, 
zeigen doch auch die friiheren Bande die gleiche wertvolle Tendenz, den nur - 
fiir angewandte Elektronik interessierten Leser aus seinem engeren Gebiet heraus- 
zufiihren und damit der iiberhandnehmenden Spezialisierung entgegenzuwirken. 
Schon aus diesem Grunde ist auch dem vorliegenden Bande ein weiter Leserkreis 
zu wiinschen. A.A. Rusterholz 


Elements of Statistics. Von H.C. Fryer (John Wiley & Sons, New York 
1954). 262 S., 34 Fig.; $4.75. 

Das Buch stellt einen modernen Lehrgang der mathematischen Statistik fiir _ 
Nichtmathematiker, die keine Vorkenntnisse auf diesem Spezialgebiet besitzen, 
dar. Vorerst werden die Begriffe Wahrscheinlichkeit und mathematische Erwar- 
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tung erlautert und die wichtigsten statistischen Masszahlen eingefiihrt. Nach der 
Besprechung der Binomial- und Normalverteilung werden verschiedene Schat- 
zungs- und Testverfahren, insbesondere x°- und ¢-Test, behandelt, dann werden 
die Fragen der linearen Regression und Korrelation besprochen. Es wird ver- 
sucht, dem Leser den zentralen Grenzwertsatz der Wahrscheinlichkeitsrechnung 
und seine Auswirkungen auf die Anwendbarkeit der statistischen Methoden klar- 
zumachen. Der Autor hat sich bemiiht, mit méglichst wenig Mathematik aus- 
zukommen, und hat um so grésseres Gewicht auf grundsatzliche Fragen gelegt. 
Vor allem wird der Leser immer wieder auf die Voraussetzungen aufmerksam 
gemacht, auf denen die behandelten Verfahren beruhen. Der Text wird durch 
sehr viele Beispiele erlautert, die aus den verschiedensten Anwendungsgebieten 
der Statistik ausgewahlt worden sind. Diesem Lehrbuch ist eine méglichst grosse 
Leserschaft vor allem in den Kreisen der Naturwissenschafter zu wiinschen. 

P. Schmid 


Progress in Low Temperature Physics, Vol. I. Herausgegeben von 
C. J. GorTER (North Holland Publishing Company, Amsterdam 1955). 418 S., 
ian bigs: htt. 30250, (61 sh,-$ 8: —). 

In 18 Artikeln schreiben 21 fiihrende Kaltephysiker iiber ihre Spezialgebiete. 
Etwa die Halfte des Buches ist den Erscheinungen der Superfluiditat und Supra- 
leitung gewidmet. Der Rest befasst sich mit Para-, Ferro- und Antiferromagne- 
tismus und einigen andern, etwas willkiirlich ausgewahlten Themen. 

Inhalt und Form der Artikel sind recht verschiedenartig. Die Themen sind 
zum Teil allgemeiner, zum Teil sehr spezieller Natur. Unabhangig davon vermit- 
teln die einen Autoren unter vollen Literaturangaben eine umfassende kritische 
Ubersicht tiber den gegenwartigen experimentellen und theoretischen Stand ihres 
Gebietes, wahrend andere sich mit knappen und skizzenhaften Darstellungen 
begniigen, deren Wert zum Teil etwas fragwiirdig ist. Fast alle Artikel betreffen 
bereits publiziertes Material. R. P. FEyNMAN’s Kapitel tiber fliissiges Helium 
jedoch ist eine neue und originelle Arbeit. Dafiir ist sie einseitig und weit davon 
entfernt, eine wirkliche Ubersicht zu geben. 

Im ganzen lasst sich sagen, dass es sich um ein heterogenes, aber interessantes 
und in mancher Hinsicht niitzliches Buch handelt. A. Thellung 


Handbuch der Physik - Encyclopedia of Physics. Bd. 15: Kaltephysik II. 
Herausgegeben von S. Friicce (Springer-Verlag, Berlin 1956). 477 S., 318 Fig.; 
DM 112.—. 

Nach einer Einfiihrung in die magnetischen Eigenschaften bei tiefen Tempe- 
raturen (J. VAN DEN HANDEL) gibt D. pr KLERK einen umfassenden Uberblick 
iiber Theorie und Methode der adiabatischen Entmagnetisierung, in welchem auf 
die Eigenschaften der verschiedenen hierfiir geeigneten Salze ausfiihrlich einge- 
gangen wird. Den experimentellen Teil der Supraleitung behandelt B. SERIN, die 
Theorie J. BARDEEN, wobei besonderes Gewicht natiirlich auf die Londonsche 
Theorie und den Zwischenzustand gelegt wird. Den Abschluss bildet eine Dar- 
stellung der Eigenschaften des fliissigen Heliums mit besonderer Beriicksichtigung 
von Helium II (K. MENDELssouN). Dem vollstandig in englischer Sprache abge- 
fassten Text liegt im allgemeinen die bis 1955 bekannte Literatur (insgesamt 
rund 1200 Zitate) zugrunde, jedoch enthalten Nachtrage auch Arbeiten von 1956. 
Leider konnte nicht mehr beriicksichtigt werden, dass es inzwischen gelungen ist, 
durch Kernentmagnetisierung Temperaturen von etwas tiber Oma Ke Zumerner= 
chen. Zusammen mit dem ersten Teilband der Kaltephysik wird so ein umfassen- 
der Uberblick der Physik der tiefsten Temperaturen gegeben, wobei der vorlie- 


asi 
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gende Band fast ausschliesslich Erscheinungen gewidmet ist, die im Gebiet des 
fliissigen Heliums oder der adiabatischen Entmagnetisierung zu beobachten sind.. 


- Im Hinblick auf die recht gediegene Arbeit der Autoren ist der hohe Preis des : 


Bandes, der manchen Interessenten vom Kauf abhalten diirfte, zu bedauern. 
P. Grassmann 


Quantum Chemistry. Von Wa LtER Kauzmann (Academic Press Inc., — 


» New York 1957), 744 S., 160 Fig., $12.00. 
Das vorliegende Buch wendet sich an Chemiker, die nur massige mathema- 
tische Kenntnisse haben und sich tiber die Anwendungen der Quantenmechanik 


in der Chemie orientieren méchten. Vorausgesetzt werden nur die elementarsten | 


Grundlagen der Analysis; die in diesem Buche ben6tigten mathematischen Hilfs- 
mittel werden in einem einleitenden Kapitel in einfachster Weise entwickelt und 
an Beispielen der Physik ausserhalb der Quantenmechanik ausfiihrlich erértert 
(138 Seiten). Nach einer nicht sehr tiefschiirfenden Einfiihrung in die Quanten- 


_.. mechanik (118 Seiten) werden atomare Systeme (98 Seiten) und die Theorie der — 
’ chemischen Bindung (142 Seiten) behandelt. Ein Exkurs iiber nichtstationare ~ 


Systeme (203 Seiten), in dem im wesentlichen die Wechselwirkung der Materie 
mit der elektromagnetischen Strahlung behandelt wird, beschliesst das Buch. Die 
Darstellung ist durchweg sehr breit. Die Anwendung gruppentheoretischer Me- 
thoden, die das Versténdnis der meisten hier behandelten Fragen ausserordent- 


lich erleichtert und vertieft hatten, wurde durchweg vermieden. Der Autor 


schreibt im Vorwort: ‘The primary aim of this book is to help chemists under- 
stand the more important concepts that quantum mechanics has introduced 


- into chemistry. A secondary aim is to examine critically the limitations of many — 


_ of these concepts’. Im Gegensatz zum ersten Ziel diirfte nun leider das zweite 
“nicht in dem wiinschbaren Masse erreicht worden sein. Es ist zwar sicherlich nicht 
- einfach, die Giiltigkeitsgrenzen der in der Quantenchemie angewandten Approxi- 


mationen aufzuweisen; eine sorgfaltigere und kritischere Behandlung dieser — 


Fragen, insbesondere bei den Problemen der chemischen Bindung, ware aber 
gerade fiir diejenigen Leser, an die sich dieses Buch wendet, ausserordentlich 
wichtig gewesen. A. Primas 


Théorie des circuits de télécommunication. Von Vitotp BELEVITCH 


-(Ch. Uystpruyst, Librairie Universitaire, Louvain 1957). 384 S., 339 Fig. ; 


SbEr. 450.22: 
Das vorliegende Buch behandelt in sehr gedrangter Form den heutigen Stand 
' der Theorie der linearen elektrischen Netzwerke. W. CauvErR, H. W. BopE und 


E. A. GUILLEMIN stehen der Arbeit zu Gevatter, doch sind viele Entwicklungen ~ 


_in origineller Weise abgewandelt, meist zum Vorteil des Lesers. Die Auswahl des 


ee: Stoffes, seine Anordnung und Behandlung verraten die Vertrautheit des Verfas- _ 
‘sers mit der Materie. Leider haben sich eine ganze Anzahl kleiner Druckfehler 


ausser den im Errata aufgefiihrten eingeschlichen. 


Neben den notwendigen einfiihrenden Definitionen werden die Wellenpara- 


. meter, die allgemeine Theorie der passiven Netzwerke, die klassische Filtertheorie, 

die moderne mit gegebenem Betriebsverhalten, die Spulen und Transformatoren 
‘und der gegengekoppelte Verstarker behandelt. Erganzende Betrachtungen zur 
Synthese von. passiven Netzwerken, Einschwingvorgangen und eine kommen- 


~. tierte Ubersicht der Literatur vervollstandigen in wertvoller Weise den Inhalt. 
_ Das Buch fiillt eine Liicke im franzésischen Sprachgebiet und kann dem theore- 


tisch interessierten Fachmann sehr niitzlich sein. H, Weber 
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Neben den in der einschlagigen Buchliteratur iiberwiegend behandelten Angelegen- 
heiten des Geld- und Zahlungsverkehrs werden hier auch Devisenausgleichsfragen, 
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nisse angeht, so gehen sie im allgemeinen nicht iiber das hinaus, was beim Absolven- 
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